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;?The  maxima/minima  of  one  response  surface  (measure  of 
effectiveness,  HOE)  is  calculated  when  the  domain  is  con¬ 
strained  by  consecutive  levels  '(aspirations)  of  a  second 
response  surface.  Thus,  the  optimal  trade-off  between  the 
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two  responses  (MOE’s)  is  displayed  over  an  entire  operating 
domain,  providing  new  and  unique  capabilities  to  evaluate 


options. 
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The  study  applies  this  development  to  a  strategic 
force  structure  modeled  in  a  linear  programming  model.  The 
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methodology  determines  the  optimal  force  structure  to  maxi¬ 
mize  one  MOE  given  the  aspiration  levels  of  the  other  MOE. 
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Therefore,  the  decision  maker  has  the  optimal  force  struc¬ 
tures  to  achieve  a  range  of  demands/aspirations  and  also 
maximize  the  second  goal. 

The  methodology  has  applications  to  multicriteria 
optimization  when  relatively  large  numbers  of  domain  con¬ 
trol  variable  combinations  are  involved.  This  technique 
may  be  a  tool  to  design  the  optimal  weapon  system  charac¬ 
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teristics  or  choose  the  optimal  mix  of  weapon  types  to 
fulfill  dual  mission  requirements. 

This  work  is  an  extension  of  previous  work  which  used 
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response  surfaces,  experimental  design,  and  mathematical 
programming  to  conduct  multi-dimensional,  sensitivity 
analysis  and  weapon  system  comparisons. 
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CHAPTER  ONE 


The  Research  Problem 

Introduction 

This  study  is  concerned  with  a  system's  optimal 
structure  of  operation  when  determined  by  multiple  measures 
of  merit.  Its  purpose  is  to  offer  a  methodology  to  aid  the 
multicriteria  decision  process  when  involved  in  a  relatively 
large  number  of  alternatives.  In  particular,  it  considers 
a  system's  optimization  of  one  measure  of  merit  (response) 
in  comparison  to  a  second  measure  of  merit  (response).  The 
research  is  an  extension  of  Palmer  W.  Smith's  work  (29,30) 
in  application  of  response  surface  methodology  to  analytical 
math  modeling. 

This  first  chapter  introduces  the  problem  and  outlines 
the  study  to  follow.  A  brief  background  and  literature  re¬ 
view  lends  further  insight  to  the  problem  and  offers  some 
recommended  references  for  the  various  methodologies  to  be 
studied.  With  this  foundation,  the  chapter  then  delineates 
the  problem  statement  and  research  objectives.  Finally,  the 
scope  of  the  research  example  is  defined.  Following  the 
study,  the  appendices  list  the  model  being  considered  for 
generating  the  surfaces,  the  regression  design,  the  regres¬ 
sion  analysis,  and  various  samples  of  computer  code  employed. 

This  study  involves  many  disciplines,  i.e.,  research 
design,  response  surface  methodology,  xe  ression  analysis, 
Lagrange  multipliers,  and  non-linear  n.  th  programming.  In 
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view  of  this,  it  is  beyond  the  scope  of  this  paper  to 
attempt  an  in-depth  review  of  all  these  disciplines.  There¬ 
fore  ,  the  literature  review  is  primarily  concerned  with  an 
overview  of  the  current  methodologies  for  multicriteria 
decision  making.  The  remainder  of  the  literature  review  is 
intended  more  as  a  descriptive  definition  of  the  methodol¬ 
ogies  that  are  employed  in  the  study.  References  to  these 
different  methodologies  will  be  made  throughout  the  text  and 
are  used  during  the  research. 


Background 

A  system's  effectiveness  is  often  determined  by  more 
than  one  measure  of  merit.  These  measures  of  merit  are 
usually  determined  by  the  objectives  of  the  system's  de¬ 
signer  and  are  not  necessarily  of  the  same  numeric  scale  or 
units  of  measure  (non-commensurate ) .  For  instance,  a 
system's  commensurate  goals  may  be  to  increase  retail  sales 
and  decrease  overhead  costs.  Since  both  these  measures  are 
monetary,  the  two  goals  can  be  combined  to  a  single  measure 
of  merit,  profit.  If,  however,  the  two  goals  were  increased 
retail  sales  and  decreased  production  time,  then  the  measure 
of  merit  would  be  money  and  time  (non-commensurate'. 
Obviously,  system  policies  which  increase  the  prior  and  de¬ 
crease  the  latter  are  sought  by  the  decision  maker.  "ovever, 


the  actual  relati  hips  between  the  various  goals 


their 


respective  optimi  ion  can,  at  times,  be  vague  ana  elusive 
for  the  decision  .  er.  Often,  the  analyst  attempts  to 
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transform  the  different  measures  of  merit  (MOE's)  into  a 


common  unit  such  as  dollars.  The  transformation,  although 
appealing,  may  raise  some  serious  questions.  Ignizio  (15: 
375)  exemplifies  this  problem  by  considering  the  dilemma  of 
minimizing  highway  funding  and  also  minimizing  highway  in¬ 
juries  and  deaths.  Few  decision  makers  will  want  the  task 
of  trading  off  dollars  for  lives.  Since  a  common  unit  is 
not  always  attainable,  then  the  system's  operation  is 
characterized  by  more  than  one  measure  of  effectiveness 
(HOE).  A  vector  of  multiple  MOE's  is  applied.  Each  vector 
element  is  the  level  of  attainment  for  one  of  the  decision 
maker's  (system's)  objectives,  i.e.,  an  MOE.  Once  the  com¬ 
bination  of  possible  MOE's  is  determined  over  the  operating 
range,  the  decision  maker  must  still  make  a  trade-off  or  a 
decision  as  to  which  is  the  optimal  solution. 

Literature  Review 

In  light  of  the  decision  maker's  problem,  several 
methodologies  have  been  developed  to  aid  in  his  thought 
process.  Some  of  the  current  methodologies  for  multiple 
comparisons  are  multiple  goal  programming  (MGP),  multiple 
attribute  utility  theory  (MAUT),  and  multiple  objectives 
optimization  theory  (MOOT). 

The  first  technique,  multiple  goal  programming,  is  an 
extension  of  linear  programming  and  is  presented  in 
Ignizio 's  text  (15).  Essentially,  multiple  goal  programming 
sequentially  solves  a  series  of  linear  programs.  However, 


M 1 V  13.  1..  ■  1 11  ■>  y.j»  5^e  p.  pi  yj  |^.  p.  pi  pj  HI 


instead  of  a  single  objective  function  and  setting  hard 
values  on  the  right-hand-side  of  the  constraint  set,  MGP 
uses  deviational  variables  in  the  objective  function  and 
also  the  constraint  set.  Therefore,  the  sequential  linear 
program  attempts  to  minimize  the  distance  or  range  of  the 
deviational  variables  from  their  right -hand-side  values 
(aspiration  levels)  (15:376).  With  this  technique,  the  MGP 
allows  the  analyst  to  sequentially  go  through  a  series  of 
lexographically  ordered  goals  and  either  minimize  or 
maximize  the  deviation  from  the  aspiration  level  set  in 
the  right-hand-side  value.  This  technique  avoids  many  of 
the  linear  programming  cases  of  infeasibility.  Since  each 
MGP  objective  is  made  of  deviational  variables,  the  algo¬ 
rithm  often  finds  how  "close*'  one  can  come  to  the  tradi¬ 
tional  linear  programming  optimal  solution  even  though  the 
actual  aspiration  level  (right-hand-side)  may  not  be 
attainable. 

However,  a  disadvantage  in  the  use  of  the  MGP  algo¬ 
rithm  may  occur  if  the  goals  require  a  mix  of  deviational 
variables  from  different  right-hand-side  values.  In  this 
case,  the  goal  in  question  requires  the  deviational  vari¬ 
ables  to  be  commensurate  in  both  common  units  of  measure 
and  equal  numeric  scale.  If  a  single  goal  contains  non- 
commensurate  deviational  variables,  the  algorithm  will 
treat  the  unit  change  in  one  as  equal  to  a  unit  change  in 
the  other.  Thus,  if  a  right-hand-side  range  for  one  de- 
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viational  variable  were  in  terms  of  tenths,  whereas  a 
second  were  in  terms  of  hundreds,  the  MGP  algorithm  would 
be  unable  to  distinguish  the  deviation  of  tenths  versus 
that  of  hundreds.  Likewise,  it  would  equally  weight  the 
value  between  unalike  measures  of  effectiveness  such  as 
time  and  money. 

Thus,  the  MGP  solution  is  a  lexographic  ordering  of  the 
goals  based  on  the  deviations  from  the  goals'  aspiration 
levels.  It  does,  as  stated  earlier,  yield  a  solution  even 
if  all  goals  are  not  attained.  This  MGP  solution  is  the 
"best/optimal"  given  the  order  of  goals  and  the  relative 
deviation  from  the  unattained  goals.  See  Ignizio  (15:16) 
for  further  multiple  goal  programming  discussion. 

The  next  two  techniques,  MAUT  and  MOOT,  were  developed 
to  handle  multicriteria  decision  theory  problems.  The 
first,  MAUT,  is  used  to  elicit  from  the  decision  maker  his 
preference  function  in  the  trade-off  between  the  system's 
measures  of  effectiveness  over  the  range  of  possible  opera¬ 
tion.  This  preference  function  applies  its  criteria  weight 
to  each  alternative  policy  and  then  searches  for  the  al¬ 
ternative  with  the  greatest  weighted  value.  This  approach 
requires  a  good  deal  of  elicitation  time  with  the  decision 
maker  and  also  a  calculation  of  each  alternative. 

The  second  approach,  MOOT,  is  a  multiple  optimization 
method  which  attempts  to  identify  a  non-dominated  solution 
set  (NDSS).  The  NDSS  is  the  set  of  solutions  or  alterna- 


tives  that  cannot  be  bettered  by  any  other  alternatives. 

This  dominance  is  based  on  an  alternative's  set  of  measures 
of  effectiveness  (MOE's).  Each  alternative's  performance 
is  evaluated  by  its  MOE's.  Thus,  associated  with  each  al¬ 
ternative,  ,  is  a  set  or  vector  of  MOE's,  v^  (11:1).  The 
non-dominated  solution  set  is  determined,  then,  by  a  com¬ 
parison  between  each  alternative's  MOE  vector  (v^).  When 
compared  to  another  MOE  vector,  each  element  of  the  non- 
dominated  vector  must  be  greater  than  or  equal  to  the  com¬ 
pared  vector  and  at  least  one  of  the  non-dominated  vector 
elements  must  be  greater  than  its  respective  counterpart  in 
the  compared  vector.  Unlike  MAUT,  this  technique  of  non- 
dominance  does  not  require  the  development  of  an  exacting 
trade-off  between  the  individual  measures  of  merit  which 
characterize  the  system's  operation.  Therefore,  the  NDSS  is 
an  optimal  set  of  alternatives  (A^)  with  an  MOE  vector 
better  than  any  other  alternative  MOE  vector.  This  solution 
set  is  also  referred  to  as  a  pareto-optimal  frontier.  Both 
Keeney  and  Raiffa  (18)  and  Zeleny  (34)  offer  texts  for 
studying  multicriteria  decision  theory.  Dr.  Dewispelare ' s 
paper  (11)  offers  an  excellent  review  of  multicriteria 
decision  theory  and  also  presents  a  combined  MOOT  non- 
dominated  solution  set  and  MAUT  preference  function  approach 
to  give  the  optimal  recommendation  to  the  decision  maker. 

In  either  case,  MOOT  or  MAUT,  the  techniques  require 
an  appreciable  time  with  the  decision  maker  and  also  require 
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his  arbitrary  choice  or  trade-off  between  the  measures  of 
merit.  Both  approaches  require  a  calculation  be  made  upon 
each  possible  alternative  which,  in  the  case  of  a  very 
complex  or  large  problem,  may  take  a  considerable  amount  of 
time. 

Besides  these  three  major  areas  (MGP,  MAUT,  and  MOOT), 
the  literature  review  also  covers  a  technique  employed  by 
the  chemical  engineering  industry.  In  1960,  Lind,  Golden, 
and  Hickman  (21)  demonstrated  the  use  of  a  research  design 
to  develop  a  response  surface  for  yield  (MOE  1)  based  on  the 
mix  of  three  possible  chemicals  (independent  domain  vari¬ 
ables).  Their  second  MOE,  cost,  was  then  calculated  for  the 
possible  combinations  of  their  independent  variables.  Thus, 
two  response  surfaces  representing  their  two  measures  of 
effectiveness  were  created.  Each  surface  was  then  reduced 
into  iso-value  contour  lines,  similar  to  those  used  in  typo¬ 
graphy  to  define  elevations.  To  then  optimize  the  system, 
the  two  sets  of  contours  were  then  overlaid  to  find  the 
maximum  ratio  of  yield  to  costs. 

In  1981,  Khuri  and  Conlon  (19)  showed  the  optimization 
of  several  responses  (MOE's)  by  use  of  a  regression  poly¬ 
nomial.  In  this  case,  the  authors  show  that  they  can 
generate  a  global  "best”  solution  over  multiple  response 
surfaces  of  the  same  degree  over  the  same  operating  domain. 
This  "best"  solution  is  generated  from  the  individual 
maxima  of  each  response.  A  distance  measure,  r,  is  devised 


to  move  from  the  original  starting  point  of  the  individual 
maxima  to  a  common  point  of  the  highest  estimated  value. 

This  distance  function  is  essentially  the  distance  moved 
from  a  maximum  divided  by  the  variance  of  said  estimate. 

Given  this  distance  measure,  the  analyst  may  then  find  the 
global  "best"  solution  for  the  multiple  MOE's  of  equal 
degree  response  surfaces. 

Thus  far,  the  literature  review  has  focused  on  the 
area  of  global  optimization.  Now  the  attention  is  shifted 
to  offering  the  decision  maker  a  range  of  optimal  solutions 
for  a  system's  operating  domain. 

P.  W.  Smith  (29,30)  has  shown,  through  the  use  of 
orthogonal  design  and  multivariate  regression,  the  ability 
to  accurately  predict  the  increase  in  a  measure  of  effective¬ 
ness  due  to  each  of  the  operating  domain  elements.  He 
demonstrates  the  optimal  rate  of  MOE  increase  per  domain 
element  in  the  orthogonal/independent  regression  coefficients. 
Likewise,  a  pairwise  trade-off  between  domain  elements  is 
determined  by  the  ratio  of  the  respective  independent  re¬ 
gression  coefficients.  Thus,  Smith's  regression  equation 
may  be  used  to  predict  which  domain  element (s)  should  be 
used  in  the  achievement  of  some  aspiration  level  set  by  the 
decision  maker.  This  solves  the  problem  of  having  to 
totally  rerun  a  set  of  linear  programs  and  thus  gives  the 
analyst  a  very  powerful  tool  if  the  aspiration  level 
should  change  for  the  decision  maker.  This  regression 
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prediction  will  be  good  as  long  as  it  is  used  within  the 
original  domain  limit  of  the  research  design.  Further  use 
of  Smith's  study  and  the  extension  to  this  research  will  be 
defined  in  Chapter  Two  for  the  proposed  methodology. 

To  adequately  review  the  proposed  methodology  requires 
preliminary  knowledge  of  several  disciplines  such  as  re¬ 
sponse  surface  methodology,  multiple  regression,  and  re¬ 
search  design.  Although  each  is  a  subject  in  itself,  a 
brief  definition  and  overview  are  offered  as  follows. 

Response  Surface  Methodology.  Response  surface  method¬ 
ology  was  originally  proposed  by  G.E.P.  Box  and  K.B.  Wilson 
in  their  article  in  the  Journal  of  the  Royal  Statistical 
Society  in  1951.  Since  then,  it  has  been  well  covered  in 
many  texts  such  as  Cochran  and  Cox  (6)  and  Myers  (22).  The 
concept  of  the  response  surface  is  as  the  name  says,  a  sur¬ 
face  such  as  that  of  the  earth.  We  may  think  of  a  surface 
as  responses  or  measures  of  effectiveness  that  are  a 
function  of  their  domain  elements.  In  the  analogy  to  the 
earth’s  surface,  we  can  think  of  elevation  as  being  the 
coordinates,  latitude  and  longitude.  Thus,  the  elevation 
changes  as  we  move  across  the  domain.  This  surface,  then, 
may  be  expressed  in  a  functional  form  where  the  dependent 
variable,  r,  is  the  measure  of  effectiveness  and  is  a 
function  of  the  independent  variables,  d^,  where  k  runs 
from  1  through  n,  n  being  the  dimension  of  the  domain  and 
n+1  being  the  dimension  of  the  entire  surface.  Eq  (1)  is 


representative  of  this  functional  relationship 


rij  Ri(dl'  d2'  ••*dn) 


(1) 


where  i  =  the  type  of  response  considered  and  j  =  the  domain 
location  (i.e./  j  =  (d^,  . ..dn)  )•  Eq  (!)  also  demon¬ 

strates  the  response  surface's  functional  relationship  for 
"i"  multiple  response  surfaces  over  the  same  domain.  Thus, 
as  in  our  original  example,  the  earth  elevation  is  a  measure 
of  merit  that  would  be  response  one  over  the  domain  of  lati¬ 
tude  and  longitude.  Likewise,  a  second  response  such  as  the 
average  annual  rainfall  could  also  be  mapped  over  the 
earth's  longitude  and  latitude.  This  would  create  a  vector 
of  two  measures  of  effectiveness,  elevation  and  average 
rainfall,  as  a  function  of  the  same  domain.  This  is  the 
same  as  MOOT's  vector  of  measures  of  effectiveness  used  for 
the  non-dominated  solution  set  determination.  Figure  la  is 
representative  of  a  typical  response  surface  where  we  have 
two  domain  elements,  d^  and  d2*  Figure  la  is  the  geometric 
representation  of  a  three-dimensional  space,  and  Figure  lb 
collapses  that  figure  into  iso-value,  two-dimensional  con¬ 
tour  lines  similar  to  those  used  in  typography.  Each  con¬ 
tour  line  is  representative  of  Eq  ( 1 )  where  the  value  of  r 
is  set  at  a  constant.  Thus,  we  can  think  of  each  contour 
as  r^,  r2#  r^f  and  so  on.  The  examples  shown  are  no  greater 
than  three-dimensional  to  allow  for  graphic  illustration. 
Obviously,  however,  from  Eq  (1),  we  can  see  that  a  response 
surface  can  be  of  any  dimension  depending  on  the  number  of 
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Fig  1.  Example  Response  Surface 


variables  in  our  domain.  Thus,  the  purpose  of  the  response 
surface  is  to  help  us  visualize,  as  it  were,  the  functional 
relationship  between  our  measure  of  effectiveness  and  the 
independent  variables  in  our  operating  domain.  There  are 
also  numerous  techniques,  such  as  gradient  search,  for  ex¬ 
ploring  the  surface  in  search  of  stationary  (critical) 
points.  These  techniques  and  further  discussion  of  re¬ 
sponse  surfaces  are  best  covered  in  Myers'  text  (22).  Next 
we  consider  the  research  design. 


11 


Research  Design.  The  research  design  is  a  powerful 
sampling  technique  used  by  the  experimenter  to  try  to  ap¬ 
proximate  the  response  surface  mentioned  above.  The  number 
of  possible  permutations  of  our  domain  set  can  grow  very 
quickly  as  the  range  of  each  domain  element  increases.  Thus 
an  attempt  to  evaluate  each  point  may  prove  tedious  if  not 
impossible  in  some  cases.  For  example,  if  a  three-element 
domain  were  to  have  respective  ranges  of  10,  15,  and  20 
units  each,  then  the  number  of  permutations  for  all  com¬ 
binations  would  exceed  some  3,600  samples.  Although  this 
number  seems  large,  consider  if  the  ranges  were  increased  to 
20,  40,  and  50  respectively.  Rather  than  3,696  samples, 
43,911  samples  would  be  required.  Obviously,  total  enumera¬ 
tion  of  the  surface  becomes  impossible  as  the  problem's 
domain  range  increases.  Thus,  the  analyst  must  approximate 
the  surface  by  sampling  different  points  across  it.  If  the 
surface  is  known  to  be  flat,  the  analyst  would  simply  sample 
each  of  the  corner  points  and  then,  connecting  these  corners 
be  able  to  completely  define  the  surface  in  question.  If, 
on  the  other  hand,  the  surface  is  not  flat  or  unknown,  this 
sampling  of  strictly  the  corner  points  is  insufficient  to 
describe  the  changing  curvature  over  the  surfaces. 

The  theory  of  research  design  is  well  documented  in 
numerous  texts  (6,14,22,32)  and  also  in  numerous  articles 
such  as  Box,  et  al  (3,4).  In  general,  the  design  is 
chosen  by  the  experimenter  based  on  some  prior  knowledge 


of  the  response  surface  or  some  assumption  about  it.  Thus, 
a  2  design  is  chosen  if  the  experimenter  believes  there  are 
no  second  order  terms  describing  the  surface.  If,  upon 
testing,  his  data  lacks  goodness  of  fit,  the  researcher  must 
re-evaluate  his  original  assumption  about  the  curvature  of 
the  surface  and  choose  a  smaller  region  to  examine  or  per- 
haps  a  higher  order  model.  The  3  design  allows  for  second 
order  terms  in  the  characteristic  equation  of  the  surface 
and  thus  would  be  the  next  choice.  This  characteristic 
equation  is  the  regression  of  the  surface  and  is  discussed 
in  the  next  definition. 

Thus,  the  research  design  is  nothing  more  than  a  choice 
in  the  number  and  location  (domain  elements'  levels)  of 
samples  to  be  taken.  The  power  of  the  design  is  its  ability 
to  represent  the  experimental  region  in  as  few  samples  as 
possible.  This  power  of  design  is  developed  by  choosing 
the  domain  levels  (sample  locations)  such  that,  when  put 
into  an  array,  the  resultant  matrix  is  orthogonal.  In  this 
design  matrix,  each  row  is  a  sample  and  each  column  is  one 
of  the  domain  elements  at  its  particular  level  for  that 
sample.  Thus,  if  we  had  a  three-element  domain,  and  ten 
samples,  we  would  have  a  ten-row  by  three-column  matrix  for 
our  design.  Appendix  E  shows  the  design  matrix  applied  in 
this  study.  The  true  power  of  the  orthogonal  design  matrix 
is  in  its  ability  to  generate  independent  coefficients  in 
the  general  linear  regression  model  mentioned  below.  To 


accomplish  this,  our  design  matrix  (D)  is  orthogonal;  thus, 

the  transpose  times  the  original  matrix  equals  a  diagonal 

matrix,  D'D  =  A*I.  There  are  many  published  designs  which 

accomplish  this  (2,3,4,  and  8).  This  study  is  primarily 
,  > 

concerned  with  a  3  fractional  factorial,  rotatable  design 
introduced  by  Box  and  Behnken  (2).  This  design  will  be 
further  discussed  in  Chapter  Two  of  our  methodology. 

Multiple  Regression.  As  mentioned  earlier,  response 
surface  is  merely  a  geometric  means  to  help  interpret  the 
functional  relationship  shown  in  Eq  ( 1 ) .  The  researcher 
chooses  the  experimental  design  in  hopes  of  achieving  a 
good  characteristic  sampling  of  the  surface  in  question. 

This  characteristic  sampling  is  then  expressed  in  a  multiple 
regression  equation.  The  multiple  regression  equation  is  a 
simplified  polynomial  expression  of  a  functional  relation¬ 
ship  depicted  in  Eq  (1).  The  accuracy  of  this  polynomial 
lies  in  its  ability  to  predict  points  on  the  surface. 

Thus,  the  goodness  of  fit  is  determined  by  the  residual 
error  between  the  polynomial  prediction  and  that  of  the 
actual  experimental  value.  This  study  is  primarily  con¬ 
cerned  with  linear  regression  models  rather  than  the 
curvilinear  regression.  DeGroot  (9:5-10)  and  Devore 
(  10:422)  both  give  excellent  theoretical  descriptions  of 
least  square  estimates  and  their  use  in  regression 
estimates.  Likewise,  Kleinbaum  (20:136)  shows  regression 
application  to  fitting  response  surfaces.  Eq  (2)  is  a 


polynomial  approximation  for  the  functional  relationship 
in  Eq  (1).  This  equation  can  also  be  expressed  as  in 
Eq  (3)  in  matrix  form,  where  r  is  a  column  vector  of  re¬ 
sponses,  B  is  the  vector  of  coefficients,  and  D  is  the  de¬ 
sign  matrix  as  described  previously.  In  minimizing  the 
error  term,  e,  in  Eq  (3),  it  has  been  shown  in  many  texts 
(9:513,  22:29)  that  the  maximum  likelihood  estimator  which 
minimizes  the  variance  is  B,  as  shown  in  Eq  (4).  The  power 

r  =  B0  +  B1d1  +  B2d2  ...Bndn  +  Bnd!2  +  B22d2  ’  ** 

2 

+  B  d  +  B, -d. d~  . . .  +B ,  w  lNdd  , 
nnn  1212  (n)(n-l)  n  n  -  1 

r  =  D  •  B  +  e 


B  =  (D'D)"1  D*r 


of  the  orthogonal  design  is  again  shown  in  the  solution  of 
B  estimate.  As  stated  earlier,  the  power  of  the  D'D  term, 
which  collapses  into  a  diagonal  identity  matrix,  is  to  create 
the  orthogonal  polynomials  which,  in  turn,  cause  each  of  the 
coefficients  in  the  B  vector  to  be  independent  from  one 
another  (31:72).  The  power  and  use  of  this  independence 
between  the  regression  coefficients  is  brought  up  and  used 
in  the  methodology  and  referenced  in  Smith's  work  also. 
Finally,  the  error  term  in  the  regression  equation  is  a 
measure  of  our  lack  of  ability  to  predict  the  true  value  on 
the  response  surface.  This  error  term  is  composed  of  two 
components,  response  variance  and  model  bias.  Wonnacott 


s 
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is  no  stochastic  variance  in  our  error  term.  Thus,  designs 
which  minimize  the  bias  induced  by  the  model  would  be  of 
optimal  choice  for  this  study.  Karson  (16,17)  shows  the 
criteria  for  a  minimum  bias,  two-dimensional  design  to  be 
orthogonality,  rotatability ,  and  equal  spacing.  However, 
there  is  no  three-level  design  which  meets  the  criteria  for 
minimizing  only  bias.  The  Box  and  Behnken  design  (2) 
chosen  for  this  study  is  designed  to  reduce  both  measure¬ 
ment  variance  and  model  bias. 

The  last  three  topics,  response  surfaces,  research  de¬ 
sign,  and  regression  analysis,  have  each  been  briefly  re¬ 
viewed  in  comparative  isolation  from  one  another.  In 
actuality,  response  surface  methodology  incorporates  all 
three  topics  simultaneously.  This  is  evidenced  by  the  in¬ 
ter-relationship  between  the  design  level  and  the  order  of 
the  regression  equation.  The  three-way  interdependence  is 
completed  by  the  regression  equation's  error  term  which 
measures  the  adequacy  of  the  researcher's  estimate  of  the 
surface. 

The  reader  should  see  Myers'  "Response  Surface 
Methodology"  (22)  for  a  combined,  in-depth  coverage  of 
these  subjects. 


Problem  Statement 


There  is  a  conceptual  problem  for  a  system's  decision 
maker  between  multiple  objectives,  their  respective  measures 
of  merit  (MOE's),  and  the  optimal  system  structure/operation 
to  achieve  these  respective  goals.  Methodologies  mentioned 
in  the  literature  review  aid  the  decision  maker  in  this 
process  but  still  depend  heavily  upon: 

(i)  The  decision  maker's  MOE  trade-off  weight? 

(ii)  The  use  of  commensurate  goals;  and 

(iii)  The  calculation  of  each  alternative  to  de¬ 
termine  the  non-dominated  solution  set . 

Although  these  problems  may  not  all  arise  in  a  particular 
study,  they  become  more  prevalent  as  the  range  and  number 
of  the  domain  elements  increase.  This  is  exemplified  by 
optimal  mix  problems  mentioned  earlier. 

Research  Objectives 

Find  the  optimal  relationship  between  a  system's  set 
of  measures  of  effectiveness  over  a  relatively  large 
operating  domain. 

Subobjectives : 

(i)  Generate  the  surfaces  for  comparison  over  the 
same  domain; 

(ii)  Apply  a  research  design  to  estimate  the 
surfaces  in  regression  equations; 

(iii)  Choose  response  surface,  r^,  to  be  constrained 


(iv) 


(vi) 


by  the  second  response  surface,  r 2 .  Represent 
the  constraining  surface,  in  iso-value  con¬ 

tour  lines; 

Determine  the  constrained  maxima  of  r^  for 
each  of  the  ^  contours; 

Map  the  constrained  maxima  function  for  r , ;  and 
Reverse  the  roles  of  r^  and  ^  and  show  the 
constrained  maxima  function  for 


Scope 

The  purpose  of  this  study  is  to  optimally  compare  mul¬ 
tiple  response  surfaces  over  the  same  operating  domain. 

Thus,  the  development  of  a  highly  sophisticated  model  to 
generate  such  surfaces  is  not  deemed  necessary.  This  study 
uses  a  deterministic,  aggregated,  arsenal  exchange  model  to 
look  at  the  problem  of  optimal  force  mix  in  strategic 
planning.  This  model  is  discussed  further  in  the  method¬ 
ology,  Chapter  Two,  and  again  in  Appendix  A. 

Since  the  true  functional  relationship,  Eq  (1),  is  not 
known,  this  study  assumes  a  second  order  relationship 
between  the  MOE  and  the  domain  variables.  This  decision  was 
based  on  the  prior  work  of  Smith  (30)  on  similar  systems. 
Smith  found  the  damage  expectancy  measure  of  effectiveness 
to  be  monotonically  increasing  and  a  second  order  regression 
polynomial  gave  excellent  results.  In  view  of  this  prior 
knowledge,  this  research  will  start  by  using  a  second 
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order  regression  and  three-level  design. 


General  Methodoloc 


Represent  one  surface  in  a  set  of  iso-value 


constraint  lines  and  determine  the  interval 


for  those  lines? 


(ii)  Find  the  constrained  maxima  of  r^  subject  to 


the  iso-constraint  of  Tc^i  and 


(iii)  Map  the  constrained  maximum  function  in  terms 


of  a  constraining  response. 


The  next  chapter  expands  the  methodology  and  discusses 


the  specific  theory  involved.  Chapter  Three  will  then  apply 


this  proposed  methodology  to  an  example.  The  analysis 


section  of  the  methodology  is  covered  in  Chapter  Four. 


CHAPTER  TWO 


Methodology 


Introduction 

This  chapter  details  the  methodology  for  the  comparison 
and  optimization  for  one  response  surface  to  another.  It 
builds  on  the  general  methodology  described  in  Chapter  One 
and  expands  on  those  areas  of  the  literature  review  which 
are  pertinent  to  the  proposed  methodology. 

This  methodology  consists  of  12  steps.  Each  step  is 
presented  both  in  general  format  and  specific  application 
for  Chapter  Three.  Prior  to  the  actual  stepwise  presenta¬ 
tion,  however,  an  overview  of  the  theoretical  basis  is 
covered.  This  overview  serves  as  the  "big  picture"  for  the 
step-by-step  optimization  process. 

Theory 

The  theoretical  basis  for  this  study  is  best  described 
by  Lagrange  multipliers.  Protter  (27:176)  demonstrates  the 
use  of  Lagrange  multipliers  to  determine  the  constrained 
maxima  of  one  function  (Rg)  over  the  domain,  d^,  subject  to 
the  constraints  of  other  functions  (R^,  R2»«.»)  over  the 
same  domain.  Eqs  (5)  through  (9)  demonstrate  the  general 
methodology  for  Lagrange  multipliers  where  RQ  ( d is  the 
objective  to  be  optimized  subject  to  the  requirements  of 


optimize  =  Rg(d^) 


(5) 
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subject  to  R^(d^)  =  0 


The  Lagrange  multipliers  are  introduced  to  form  a  new  func¬ 


tion  R^(d^/  X  j )  where  d^  is  the  original  domain  and  X  j  is 


the  multiplier  for  the  j —  constraint. 


optimize  R£(di#  =  R0(di}  +  X  j  *  Rj^di 


Thus,  the  problem  is  now  to  optimize  the  new  function 


R^(d^,  ^  j ) •  This  is  accomplished  by  simultaneously  solving 


the  set  of  partial  derivatives  of  Ri  set  equal  to  0. 


solve 


R q/0  d^  =  0  for  all  i's 


R^/0  X  j  =  0  for  all  j 's 


The  solution  to  this  set  of  equations  is  the  set  of  critical 


(stationary)  points  for  RQ  subject  to  Rj.  Protter  (27:178) 


gives  several  numeric  examples  of  this  technique. 


Chapter  One's  literature  review  describes  the  response 


surface  graphically  in  Figure  1/  and  functionally  in  Eq  (1). 


The  review  also  describes  the  researcher's  approximation 


for  the  response  surface  functional  relationship  in  a  multi- 


variable  regression  polynomial ,  Eq  (2).  The  combination  of 


these  techniques  with  Lagrange  multiplier  optimization  lays 


the  foundation  for  this  study. 


Figure  2a  depicts  the  surface  to  be  optimized,  R^djO. 
Figure  2b  depicts  the  constraint  surface  R^(d^)  in  both  the 
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n+1  problem  dimension  and  collapsed  into  contour  lines  in 
the  n  dimensional  domain.  Thus,  each  constraint  contour 


is  a  subset  of  the  total  surface  R^(d^)  with  a  constant 

value  r  .  Therefore,  the  surface  R.  (d.)  is  represented 
1  j  1 

as  a  set  of  contour  lines  in  the  operating  domain.  Eq  (10) 

depicts  this  set  of  j  constraint  contours.  The  number 

rl 

of  contours  required  to  characterize  the  surface  is  de¬ 
pendent  on  the  surface  complexity. 
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The  question  of  contour  interval  criteria  is  determined  in 
step  9  of  the  Methodology  presentation.  Figures  2a  and  2b 
pictorially  define  the  two  measures  of  merit  (MOE)  which  are 
to  be  optimized  one  to  the  other.  Figure  2c  combines 
Figure  2a ' s  surface,  RQ(d^),  and  Figure  2b' s  contours, 
r^  =  R1(d^).  The  final  step  is  to  find  a  maximum  value  of 
Rg(d^)  when  constrained  to  each  contour,  r^  =  R^d^).  Thus, 
there  are  j  sets  of  Lagrange  multiplier  constrained  maxima 
problems  to  be  solved,  Eq  (11). 


optimize  each  rQ  =  R£(dA,  =  Ro^di^  +  ^  j  ’  (Ri(dj.) 

•  V 


The  solution  set  (rQ  ,  r  )  depicts  the  optimal  RQ(d^)  for 

j  J 

a  given  R^(d^)  value  over  the  entire  operating  domain  (d^). 

This  theoretical  presentation  must  be  further  constrained 


(10) 


(11) 


for  most  applications.  The  theory  allows  the  range  of 
each  domain  element,  d^,  to  be  unlimited.  In  practical 
application,  the  domain  is  limited  to  some  experimental/ 
operating  range.  In  fact,  it  is  this  range  which  defines  the 
research  design  and,  therefore,  limits  the  area  to  which  the 
regression  polynomial  may  be  applied  (6:336).  This  problem 
is  exemplified  in  the  case  of  the  polynomials  approximating 
each  surface  being  second  order.  In  this  case,  the  Lagrange 
optimization  will  seek  to  reach  the  global  critical  point. 
This  point,  although  mathematically  correct,  may  lie  out¬ 
side  the  design  region  under  consideration.  Therefore,  the 
upper  and  lower  bounds  for  each  domain  element  must  be  in¬ 
cluded  in  the  problem,  Eq  (12).  These  additional  con¬ 
straints  require  a  Lagrange  multiplier  for  each  bound. 

Ui  ^  di  <  1. 

Thus,  for  n  domain  elements,  an  additional  2*n  X  's  must  be 
included  in  the  optimization  function,  R^(d^,  X j,  Xu /  X ^  ) . 

The  methodology  to  solve  this  new  inequality  con¬ 
strained  optimization  function  is  shown  in  Hadley  &  Whitin' s 
text  (13:436).  Essentially,  the  process  involves  a  pairwise 
comparison  between  inequality  constraints  and  the  original 
problem  to  determine  which  constraints  are  active.  An 
active  constraint  is  one  which  actively  restricts  the  optimal 
domain  of  the  original  problem.  These  active  constraints  are 
set  to  equalities  and  the  problem  is  resolved.  The  inactive 


constraints  are  less  restrictive  than  the  active  ones  and 


therefore  can  be  ignored.  "The  effort  required  to  solve  the 
problem  increases  rapidly  with  the  number  of  [inequality] 
constraints/"  (13:437).  Due  to  this  increased  effort,  an 
alternative  approach  is  also  offered  in  step  10  of  the 
methodology. 

In  conclusion,  this  section  has  shown  the  theoretical 
basis  for  the  proposed  methodology.  Figure  2  has  shown  the 
two  MOE's  (rQ,  r^)  in  question  as  response  surfaces  (RQ,  R^ ) 
and  collapsed  the  constraining  surface  into  iso-value  con¬ 
straint  lines  (r^  =  R1(d^)  )  in  the  operating  domain. 

Finally,  the  constrained  maxima  of  the  optimization  surface 
subject  to  each  contour  is  found  using  Lagrange  multiplier 
techniques.  It  is  the  solution  set  of  j  maximized,  ordered 
pairs  which  depict  the  optimal  trade-off  between  the  two 
surfaces.  The  remainder  of  this  chapter  is  dedicated  to 
the  actual  step-by-step  process  to  achieve  this  optimization. 

Stepwise  Methodology 

Figure  3  is  a  flowchart  of  the  proposed  methodology. 

It  serves  as  an  outline  to  this  chapter  and  the  entire  study. 
The  first  three  steps  are  systems  analysis  techniques  and  the 
remaining  steps  are  optimization  specific.  Obviously, 
these  first  three  steps  cannot  summarize  the  many  systems 
analysis  texts.  However,  this  study  is  primarily  concerned 
with  an  optimization  technique;  these  three  steps  are  in- 
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eluded  to  re-emphasize  their  importance  in  a  complete  re¬ 
search  project. 

Step  1  -  Problem  Definition.  The  problem  statement 
acts  as  the  keystone  and  focus  for  the  research.  In  de¬ 
fining  the  research  problem,  Quade  (28:69)  offers  a  series 
of  questions  to  help  the  analyst  interrogate  the  system's 
designers/decision  makers.  The  purpose  is  to  construct  as 
lucid  and  specific  a  task  definition  as  possible.  One 
question  bears  reiteration  throughout  the  study: 

"Is  it  the  right  problem  anyway?  Might  it 
not  be  just  a  manifestation  or  a  symptom  of  a 
much  larger  or  deeper  problem?  Would  it  be 
better  to  tackle  this  larger  problem  if  there 
is  one?"  (Quade  28:69) 

As  the  research  progresses,  this  question  must  be  reiterated 
since  greater  problem  insight  may  require  a  redefinition  of 
the  research  task.  It  is  this  iterative  process  which  the 
flowchart.  Figure  3,  depicts  with  arrows  returning  to  step 
one. 

For  Chapter  Three's  example,  the  problem  is  to  find 
the  optimal  mix  of  a  given  set  of  strategic  force  types  to 
fulfill  two  missions.  The  first  mission  is  to  maximize  the 
counterforce  (CF)  capability  given  the  force  must  meet  the 
countervalue  (CV)  requirements.  The  second  mission  reverses 
the  roles  of  the  first  and  considers  the  maximization  of  CV 
given  the  required  level  of  CF. 

Step  2  -  Variable  Selection.  The  first  step  focuses 
the  research  on  either  the  system's  inability  to  reach  an 
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aspiration  level  (goal)  or  a  new  goal  of  that  system.  Thus, 
step  two  begins  by  determining  the  system's  objectives  and 
their  respective  goals/aspiration  levels  set  by  the  decision 
maker.  Each  goal's  level  of  achievement  is  a  measure  of 
effectiveness  and  becomes  the  problem's  dependent  variable. 

The  independent/control  variables  are  defined  by  the  domain 
in  which  the  system  is  operating.  One  of  the  major  tasks  of 
the  analyst  will  be  to  determine  which  of  the  domain  elements 
have  a  dynamic  effect  upon  the  measure  of  effectiveness  under 
consideration.  Those  that  are  correlated  to  the  measures  of 
effectiveness  will  become  the  control  variables,  while 
those  that  are  not  will  be  treated  as  constants.  If,  later 
in  the  research,  data  does  not  appear  to  be  reacting  properly, 
it  may  be  caused  by  an  actual  control  variable  that  was 
thought  to  be  a  constant  and  left  out  of  the  domain.  An 
alternative  perspective  is  to  redefine  the  total  system  in 
terms  of  the  model  being  assumed.  This  model  consists  of 
strictly  the  MOE's  under  consideration  and  those  factors 
believed  to  be  the  primary  cause  of  MOE  changes.  Thus,  all 
other  elements  of  the  domain  now  become  part  of  the  environ¬ 
ment  in  which  this  model  operates. 

The  next  section  discusses  the  problems  associated  with 
scenario  dependence  in  the  selection  of  variables  and  their 
respective  ranges. 

Chapter  Three's  measures  of  effectiveness  are  levels 
of  counterforce  damage  and  countervalue  damage  as  determined 


by  the  aggregate  arsenal  exchange  model.  The  control 
variables  are  the  amount  of  each  type  of  weapon  available. 
For  this  example,  four  out  of  five  weapons  are  ranged  and 
one  weapon  is  held  constant.  The  number  of  variables  not 
included  as  control  variables  for  this  study  is  large  when 
considering  the  variables  of  the  linear  programming  model, 
each  of  the  right-hand-side  values,  and  each  of  the  equa¬ 
tions  in  the  matrix.  Those  variables  not  included  in  the 
study  are  discussed  more  thoroughly  in  steps  four  and  five. 

Step  3  -  Operational  Limits  of  Variables.  This  step 
defines  the  range  of  each  control  variable.  Step  two’s 
variable  selection  combined  with  these  ranges  define  the 
system  for  the  research  study  and  also  the  environment  or 
scenario  under  which  it  will  operate.  Quade  (28:193)  dis¬ 
cusses  the  importance  of  this  scenario/environment  to  the 
overall  analysis.  He  discusses  how  an  analysis  may  be 
biased  by  the  choice  of  scenario.  Likewise,  a  well  chosen 
scenario  may  aid  in  the  robustness  of  the  study  to  changes 
from  the  initial  problem  statement  in  step  one.  Therefore, 
the  research  should  iterate  back  to  step  one  to  insure  that 
the  problem  definition  and  research  objective  can  be  de¬ 
fined  using  the  choice  of  variables  in  step  two  and  the 
range  of  the  control  variables  in  step  three. 

For  the  example  in  Chapter  Three,  the  range  of  the 
control  variables  is  first  selected  to  start  at  zero  and 
run  to  the  upper  limits  of  450,  750,  1040,  and  424  for  each 


of  the  four  weapons  systems  respectively.  As  we  will  see 
from  the  example,  the  lower  limits  of  this  problem  were 
changed  to  300,  450,  720,  and  200  respectively,  after 
initial  analysis. 

Step  4  -  Data  Source/Generator.  In  any  research,  the 
question  of  what  resources  will  be  available  to  the  analyst 
must  be  addressed  early  in  the  study.  These  resources  will 
include  the  personnel  available,  computer-time ,  the  ex¬ 
perimental  or  previous  studies'  data,  and  existing  model 
availability.  It  will  be  from  these  resources  that  the  data 
for  analysis  will  be  generated.  The  choice  between  re¬ 
sources  will  be  determined  by  both  the  cost  and  time  limita¬ 
tion  upon  the  study  and  also  the  level  of  detail  and  variance 
acceptable  for  the  solution. 

For  this  research,  the  model  is  a  deterministic, 
aggregated  arsenal  exchange  model.  The  model  uses  linear 
programming  to  solve  for  the  optimum  use  of  weapons  using 
damage  expectancy.  The  model  considers  ten  types  of  targets 
and  five  types  of  weapons.  The  objective  function  seeks  to 
maximize  the  damage  expectancy  given  the  effectiveness  of 
each  type  of  weapon  on  each  type  of  target  and  the  capability 
to  strike  a  target  twice  with  the  same  weapon  type.  Further 
details  of  the  model  used  to  generate  the  data  are  available 
in  Appendix  A. 

Since  a  simplex  algorithm  is  generating  each  data 
point,  the  response  surface  is  actually  a  surface  of  all 


optimal  points.  Thus,  for  each  n-tuple  of  the  domain 
variables,  there  is  a  vide  range  of  possible  damage  ex¬ 
pectancy  which  could  be  considered,  but  the  response 
surface  only  contains  the  optimal  DE  from  each  n-tuple  of 
weapons . 

As  previously  mentioned,  the  choice  of  control  variable 
(weapons  availability)  is  not  inclusive  of  all  possible 
variables  in  the  linear  programming  (LP)  algorithm. 
Therefore,  the  true  functional  relationship  between  the 
damage  expectancy  and  the  domain  is  not  known  due  to  the 
complexity  of  the  LP  algorithm.  The  LP  domain  consists  of 
the  study's  control  variables,  the  other  right-hand-side 
values,  all  tableau  variables,  the  number  and  type  of  equa¬ 
tions  in  the  matrix,  and  finally  the  objective  function. 
Chapter  One,  Eq  (1),  depicts  this  actual  functional  rela¬ 
tionship  of  the  response  over  its  LP  domain.  This  actual 
relationship  is  unknown  to  the  researcher  and,  therefore, 
the  assumed  regression  model  is  chosen  to  estimate  the 
relationship.  In  the  next  chapter's  research  example,  only 
the  control  variables  (right-hand-side  weapon  availability) 
are  varied.  Thus,  the  example  regression  approximates  the 
true  functional  relationship  of  optimal  DE  with  only  the 
four  weapons  availability  variables.  This  is,  however,  a 
good  estimation  of  the  surface,  since  changes  across  it  will 
be  a  function  of  only  those  variables  which  are  changing 
in  the  overall  LP  algorithm.  In  this  case,  the  four  out  of 


five  weapons  chosen  for  the  study  are  the  control  vari¬ 
ables.  The  effect  of  the  constant  fifth  weapon  is  displayed 


in  the  Bq  coefficient  of  the  regression  polynomial. 

Thus,  a  no-variance  surface  of  optimal  DE  points  is 
generated  from  the  LP  arsenal  model  and  estimated  in  the 
four-variable  regression  model. 

Step  5  -  Initial  Surface  Estimation.  The  researcher's 
initial  estimate  for  the  regression  model  of  the  response 
surface  should  be  determined  by  some  prior  knowledge  of  the 
actual  system  or  a  similar  system.  If  the  system's  change 
in  response  is  unknown,  then  either  a  simple,  first  order 
model,  which  may  have  a  poor  fit,  is  chosen,  or  a  higher 
order  model,  which  may  require  unnecessary  sampling,  is 
chosen. 

For  this  study,  prior  knowledge  of  similar  systems 
was  available  from  Smith  and  Mellichamp's  work  (30).  They 
found,  in  studying  a  similar  model,  that  a  second  order 
regression  polynomial  gave  an  excellent  prediction  of  the 
damage  expectancy  caused  in  a  five-weapon  arsenal.  Based 
on  this  prior  knowledge,  the  study  estimates  the  example 
arsenal  will  respond  in  a  similar  manner  and  the  damage 
expectancy  will  be  a  monotonically  increasing  function 
which  may  require  a  second  order  regression  polynomial. 

Step  6  -  Research  Design  Selection.  Chapter  One 
displays  the  power  of  the  research  design  as  a  sampling 
tool  to  gather  data  points  for  use  in  the  regression 
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estimate.  This  step  now  defines  the  criteria  to  choose  a 
design.  Smith  and  Mellichamp  (30)  displayed  the  analytic 
advantage  of  independent  co-efficients  for  the  comparative 
trade-off  between  control  variables  with  respect  to  the  HOE. 
Therefore,  one  of  the  selection  criteria  will  be  orthogonality. 

A  further  criterion  is  brought  about  by  the  use  of  a 
deterministic  model.  As  shown  in  Chapter  One,  the  de¬ 
terministic  model  has  no  measurement  variance  and,  thus,  a 
primary  criterion  now  becomes  that  of  minimizing  the  bias. 
Unfortunately,  Karson's  minimum  bias  criteria  is  only  shown 
for  a  two-level  design  (17)  and  no  minimum  bias,  only  3 
designs,  exist.  Hence,  from  the  previous  step  we  will  need 
a  design  which  will  allow  for  a  second  order  equation  and 
orthogonality. 

This  criterion  is  satisfied  by  Box  and  Behnken's 

1C  k 

fractional  factorial  3  design  (2).  This  3  design  samples 
points  from  each  corner  of  a  cube  and  each  point  of  a  star 
which  is  co-located  with  the  surface  of  the  cube  and  at  the 
center  of  the  cube.  Appendix  B  shows  this  research  design 
and  the  value  of  each  of  the  measures  of  effectiveness  for 
each  sample  point.  The  design  is  given  in  two  formats, 
coded  and  uncoded.  The  coded  format  has  transformed  the 
range  of  each  domain  variable  from  upper  and  lower  limits 
to  plus  and  minus  one  respectively,  and  zero  for  the  mid¬ 
point.  Eq  (13)  shows  this  transformation  where  x  is  the 
mid-point  in  the  range  between  the  upper  and  lower  limit 


and  d  is  the  actual  range.  Thus,  for  the  example  in 


Chapter  Three,  the  first  domain  variable,  v^,  ranges  from 
300  to  450.  For  this  case,  x  is  375  and  d  is  150. 

coded  x  =  (x  -  x)/d/2 

Using  Eq  (13),  we  see  that  the  lover  limit  of  300  is  trans¬ 
formed  to  a  -1  and  the  upper  limit  of  450  to  +1  and  the  mid¬ 
point  to  zero.  The  transformed  design  simplifies  checking 
for  the  orthogonality  of  the  research  design.  Likewise,  if 
said  regression  algorithm  does  not  have  an  automatic  co- 
variance  computing  tecnnique,  the  transformed  matrix  may  be 
used . 

For  example,  in  Chapter  Three,  the  power  of  this  design 
can  be  displayed  in  the  number  of  samples  which  will  be 
required  for  the  combinations  of  the  ranges  of  each  of  the 
domain  variables.  If  a  full  enumeration  of  the  surface  were 
done,  it  would  require  greater  than  three  billion  samples. 
However,  the  example  research  design  uses  25  points.  It 
should  be  noted  that  the  original  Box  and  Behnken  design 
required  an  additional  two  repetitions  on  the  center  point 
for  degrees  of  freedom  to  estimate  the  variance  error. 
However,  this  repetition  of  the  center  point  is  not  re¬ 
quired  since  this  example  uses  a  deterministic  model  with 
only  bias  error. 

Step  7  -  Global  Regression  Estimate.  Eq  ( 1 )  shows  the 
functional  relationship  betwen  the  measure  of  effectiveness 
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and  its  control  variables.  This  step  now  estimates  that 
functional  relationship  in  a  regression  polynomial  using 
the  sampled  points  from  the  design  in  step  six.  There  are 
many  computer  regression  packages  available  to  do  this. 

For  the  case  in  point,  the  Statistical  Package  for  the 
Social  Sciences  (SPSS)  (25:320)  is  used  to  calculate  the 
regression  equation.  A  stepwise  regression  option  is 
chosen,  thus,  allowing  the  researcher  to  follow  the  intro¬ 
duction  of  each  independent  coefficient  and  re-check  for 
this  independence.  A  sample  of  the  regression  output  is 
found  in  Appendix  C.  It  should  be  noted  that  since  we  are 
using  a  deterministic  model,  the  F  statistic  is  meaningless 
except  as  a  ratio  of  the  sum  of  squares  explained  to  the 
sum  of  squares  unexplained,  which  may  be  an  indicator  of 
the  bias  error. 

The  normal  statistical  inferences  drawn  from  the  F 
statistic  cannot  be  used  in  this  deterministic  regression 
model.  The  reason  lies  in  the  basic  assumption  that  there 
is  a  distribution  of  y  for  each  value  of  x  in  the  regression 
(33:15).  Thus,  the  assumption  can  be  thought  of  as  a  con¬ 
ditional  probability,  P(Y^  x^ ) .  However,  with  the  use  of  a 
deterministic  model,  this  conditional  probability  equals 
one.  Therefore,  the  error  term  in  our  generalized  regres¬ 
sion  equation  is  not  a  measure  of  both  the  measurement 
variance  and  the  model  bias  but,  rather,  only  the  bias  due 
to  the  assumed  model  and  its  inability  to  fit  the  actual 


response  surface.  Since  the  F  statistic  is  a  ratio  of  the 
variance  explained  by  the  regression  to  the  unexplained 
variance,  we  see  that  with  the  use  of  this  deterministic 
model,  the  unexplained  variance  is  actually  a  measure  of 
the  lack  of  fit  of  our  regression  polynomial  to  the  surface 
divided  by  the  degrees  of  freedom,  Eq  (14).  Since  the  model 
has  no  variance  in  measurement,  it  would  be  wrong  to  make 
statistical  inferences  as  to  the  sample  population  from  the 
F  statistic. 

F  =  Var,  Explained  _  (Yj  ~  Y)  /p _ 

Var,  Unexplained  ,a  .  2,,  s 

(Yl  ~  Yj_ )  /(n-p) 

where 

p  =  #  of  regression  variables 

Step  8  -  Regression  Error.  This  step  tests  the  re¬ 
gression  polynomial's  goodness  of  fit  for  each  surface.  As 
stated  earlier,  the  normal  F  statistical  inferences  cannot 
be  used.  Therefore,  the  primary  criteria  will  be  residual 
error  between  the  predicted  and  actual  values.  The  residual 
error  predicted  for  each  of  the  design  points  is  included  in 
the  SPSS  output.  To  increase  the  confidence  in  the  regres¬ 
sion  equation,  a  random  sampling  of  residual  error  is  also 
taken. 

Step  9  -  Contour  Interval  Determination.  As  mentioned 


in  Chapter  One,  we  must  represent  our  constraining  surface 
as  a  set  of  iso-value  contour  lines.  The  criteria  as  to  how 
many  lines  or  the  interval  between  these  lines  lies  pri- 


marily  in  the  last  two  steps.  If  the  error  is  small  in  our 
prediction  from  our  regression  polynomial,  then  the  in¬ 
terval  between  our  contour  lines  may,  likewise,  be  small. 

The  purpose  of  the  contour  is  to  show  the  change  in  slope 
of  the  surface.  Thus,  if  we  are  unable  to  distinguish  a 
change  in  slope,  there  is  no  point  in  having  a  contour  at 
that  point  in  the  projection  of  the  surface.  Using  this 
logic,  we  can  look  at  the  standard  deviation  of  the  error 
from  steps  seven  and  eight  above.  If  we  assume  that  we  will 
be  unable  to  distinguish  a  change  in  response  unless  it  is 
outside  of  the  error  of  the  estimate,  then  we  would  not 
want  a  contour  interval  any  closer  than  that  of  the  error. 
Thus,  an  interval  criteria  can  be  established  based  on  the 
standard  deviation  of  residual  error.  Since  the  regression 
has  only  bias  error  of  a  least  square  fit  of  the  design 
data,  the  bias  should  be  a  normal  distribution.  Using 
this  assumption,  a  range  of  two  standard  deviations  of 
residual  error  would  be  a  minimum  interval  for  selection. 
Outside  this  range,  one  would  be  assured  that  at  least  97% 
of  the  time,  the  model  will  be  able  to  estimate  a  change  in 
the  response  surface  that  was  not  caused  by  the  regression 
error.  This  criteria  then  sets  a  minimum  value  for  the 
contour  interval.  The  maximum  interval  must  be  determined 
by  the  complexity  of  the  surface  and  the  degree  of 
accuracy  required  in  the  solution. 

For  the  example  in  Chapter  Three,  the  standard  devia- 


tional  errors  were  so  small  that  the  interval  was  arbitrar¬ 
ily  chosen  at  a  value  larger  than  the  minimum  that  could  be 
allowed.  This  arbitrary  choice  was  based  on  an  estimate 
believed  to  be  a  sufficient  number  of  samples  to  adequately 
define  the  trade-off  between  the  two  measures  of  effective¬ 
ness  . 

Step  10  -  Constrained  Maxima  Determination.  As  de¬ 
scribed  earlier  in  the  chapter,  Lagrange  multipliers  can  now 
determine  the  maximum  of  response  Rg,  as  compared  to  each 
contour  line  found  in  the  preceding  step.  However,  it  was 
also  shown  that  the  Lagrange  multiplier  technique  could 
become  tedious  with  the  number  of  multipliers  required  for 
the  inequality  constraints.  Thus,  with  the  number  of  con¬ 
straints  and  no  computer  libraries  readily  available  to 
solve  the  Lagrange  technique,  an  alternative  optimization 
technique  was  found.  This  technique  is  based  on  the  work  of 
Fiacco  and  McCormick  (12)  in  nonlinear  programming.  The 
actual  computer  library  that  is  used  was  developed  by 
Mylander  (23).  This  algorithm  finds  the  minimum  of  either 
a  linear  or  nonlinear,  multivariate  objective  function, 
subject  to  linear  or  nonlinear  equality  and  inequality 
constraints.  The  algorithm  uses  Fiacco  and  McCormick's 
modified  objective  function  which  is  formulated  by  using 
the  original  objective  function  and  a  penalty  function. 

It  then  starts  an  iterative  process  in  choosing  new 
estimates  to  minimize  this  new  objective  function.  The 
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routine  estimates  each  new  iterative  answer  by  an  extra¬ 
polation  process.  It  then  tests  for  convergence.  In 
seeking  its  next  iterative  solution,  it  also  uses  a 
gradient  search  routine  to  speed  the  process.  This  is  a 
powerful  algorithm  which  solves  Fiacco  and  McCormick's 
function  and  also  simultaneously  solves  a  dual  problem 
based  on  Lagrange's  technique.  For  further  discussion  of 
the  theory,  see  Chapter  8  of  Fiacco  and  McCormick  (12)  and 
also  Appendix  D,  a  description  of  the  SUMT  library  as 
composed  by  Captain  R.  M.  Floyd  for  the  Air  Force 
Institute  of  Technology. 

Step  11  -  Map  Solution  Set.  The  solution  set  consists 
of  the  maximum  points  for  each  constraint  contour, 

(max  rQ,  )  j  -  The  trade-off  between  the  two  measures  of 
effectiveness  may  be  displayed  graphically  by  plotting  the 
solution  pairs.  The  slope  of  this  graph  yields  the  trade¬ 
off  of  marginal  return  for  the  maximum  of  rQ  given  a  value 
of  r i .  Likewise,  the  n- tuples  which  created  this  optimal 
set  and  which  also  come  from  the  SUMT  program  can  be  plotted 
allowing  a  visual  picture  of  the  correlation  between  each 
domain  variable  and  its  use  in  this  trade-off. 

Step  12  -  Analysis.  This  final  step  is  the  interpreta¬ 
tion  of  the  data  that  is  gathered  in  the  past  eleven  steps. 
The  data  mapping  of  step  eleven,  when  compared  to  the  in¬ 
dependent  regression  coefficients  of  each  weapon  system, 
will  give  insight  into  the  effectiveness  of  each  control 
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variable  upon  the  single  measure  of  merit,  the  multiple 
measures  of  merit,  and  the  optimal  trade-off  in  both. 
Chapter  Four  of  this  study  will  show  the  comparisons  and 
insights  through  this  technique. 


CHAPTER  THREE 


Application  Example 

Introduction 

This  chapter  lends  numeric  reality  to  Chapter  Two's 
proposed  methodology.  The  application  is  a  generic  mili¬ 
tary  force  structure  problem  similar  to  that  studied  by 
P.  W.  Smith  and  J.  Mellichamp  (30).  All  but  the  last  two 
steps  of  the  methodology  will  be  addressed  in  this  chapter. 
The  data  mapping  and  analysis  are  postponed  until  Chapter 
Four.  The  chapter  presentation  closely  follows  the  actual 
chronological  order  of  the  research. 

As  stated  earlier,  this  generic  force  structure  prob¬ 
lem  is  similar  to  that  studied  by  Smith  and  Mellichamp. 

Their  research  shows  the  individual  effectiveness  of  each 
weapon  system  when  used  in  concert  with  the  remaining 
arsenal  available.  Their  study  uses  a  single  measure  of 
effectiveness  (HOE),  damage  expectancy  (DE).  Their  use  of 
independent  regression  coefficients  shows  both  the  in¬ 
cremental  MOE  value  of  each  weapon  and  the  marginal  return 
of  one  weapon  to  another.  Building  on  Smith  and  Mellichamp' 
work,  this  study  considers  a  similar  force  and  then  compares 
the  optimal  trade-off  between  two  missions  (MOE's), 
counterforce  (CF)  and  countervalue  (CV).  The  optimization 
of  this  trade-off  is  shown  in  the  following  methodology 
application. 


Problem  Definit ion 


The  decision  maker  must  determine  the  best  mix  of 
forces  to  fulfill  two  missions.  The  first  is  a  counter¬ 
force  (CF)  mission  and  the  second  is  a  countervalue  (CV) 
mission.  The  force  structure  must  be  optimized  so  that  it 
has  the  flexibility  to  perform  one  mission  or  the  other. 
The  alternatives  for  the  structure  are  the  possible  com¬ 
binations  of  five  different  weapon  systems'  respective 
ranges  of  values.  The  fifth  weapon  system  is  set  at  a 
fixed  level  of  360  warheads  prior  to  this  analysis.  The 
decision  maker's  problem  is  to  determine  the  optimal  force 
structure  for  a  range  of  different  CV  demands  (mission 
requirements)  and  for  the  remaining  four  weapon  systems. 

The  generic  weapon  characteristics  are  shown  in 
Table  I.  The  weapons  (w^)  may  be  thought  of  as  representa 
tive  components  of  an  arsenal.  For  instance,  w^  and  w2 
may  be  thought  of  as  two  different  types  of  ICBM's  and  w^ 
may  be  an  SLBM  component.  For  this  problem,  w^  is  assumed 
to  be  introduced  at  a  constant  level  of  360  warheads.  The 
two  characteristics,  yield  and  circular  error  probable 
(CEP),  define  the  level  of  destructive  force  on  a  target. 
CEP  is  the  radius  within  which  the  guidance  of  the  system 
will  land  the  warhead  fifty  percent  of  the  time.  CEP  and 
weapon  yield  are  given  in  thousands  of  feet  and  megatons 
respectively. 


W! 

W 

2 

w3 

W4 

W5 

CEP 

.200 

.150 

.240 

.600 

.200 

Warheads/Weapon 

1 

3 

10 

4 

6 

Yield 

2.0 

.17 

.05 

.35 

.20 

Reliability 

.85 

.85 

.85 

.85 

.85 

Daily  Alert 

.85 

.85 

.85 

.33 

.33 

Max  Warheads/Type 

450 

750 

1040 

424 

360 

The  target  base  defines  the  remainder  of  the  system. 

Table  II  shows  the  characteristics  of  the  targets.  There 

are  six  targets  in  the  CF  class  and  four  in  the  CV  class. 

The  CF  type  targets  may  be  representative  of  such  things 

3  ...  . 

as  ICBM  silos,  military  bases,  and  C  I  facilities.  Likewise, 
the  CV  type  targets  may  be  representative  of  such  things 
as  storage  or  industrial  facilities.  The  vulnerability 
number  in  Table  II  is  the  connection  between  the  target  and 
the  weapon. 

The  vulnerability  number  (VNTK)  consists  of  five 
alpha  numeric  characters  (1:34).  The  first  two  numeric 
characters  represent  the  target's  level  of  hardness  to  its 
most  vulnerable  kill  mechanism.  The  alpha  character  identi¬ 
fies  the  kill  mechanism  most  effective  against  that  target. 


WTO 
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TABLE  II 

Target  Base  Characteristics 


Type 

Number 

Vulnerability 

Number 

Diameter 

CV1 

140 

24Q0 

.  56 

CV2 

215 

1 3P1 

.49 

CV3 

430 

3 1P6 

o 

o 

• 

CV4 

520 

23Q0 

.31 

CF1 

450 

35P7 

.00 

CF2 

1000 

52P8 

.00 

CF3 

200 

39P0 

• 

o 

o 

CF4 

20 

22P1 

.36 

CF5 

150 

11P0 

o 

o 

• 

CF6 

100 

10Q1 

.79 

Finally,  the  last  number  shows  the  target's  susceptibility 
to  increased  duration  of  the  particular  effect.  The  linear 
programming  arsenal  exchange  model  requires  these  VNTK 1 s  to 
be  converted  into  single-shot  probability  of  kill  (SSP^)  co¬ 
efficients  for  the  objective  functions.  The  arsenal  exchange 
work  of  Robert  Bunnell  and  Richard  Takacs  (5)  required  the 
computer  coding  of  this  conversion  from  VNTK  to  SSP^.  With 
their  kind  help,  this  program  has  been  employed  to  compute 
the  SSP^'s  for  this  study.  Further  explanation  of  the 


target  characteristics  and  weapons  characteristics  used 
in  this  aggregated  arsenal  exchange  model  are  given  in 
Appendix  A. 

Variable  Selection 

As  mentioned  in  the  previous  chapter,  there  are  many 
possible  variables  in  the  arsenal  exchange  linear  programming 
model.  For  this  study,  we  are  primarily  concerned  with  the 
objective  function  value  and  the  right-hand-side  value  for 
the  number  of  weapons  available.  The  ability  to  destroy  the 
target  base  in  question  is  measured  by  the  objective  function 
value.  Since  two  separate  missions  are  considered,  CF  and 
CV,  the  measure  of  effectiveness  is  the  level  of  destruction 
to  each  of  the  respective  target  bases.  This  study  will 
refer  to  each  MOE  as  CF  or  CV  respectively,  where  each  is  a 
measure  of  the  objective  function  for  that  particular  model. 
Likewise,  for  consideration  of  a  combined  mission  of  both 
CF  and  CV,  this  study  will  refer  to  this  as  all  forces  (AF). 

The  control  variables  will  be  the  five  weapon  systems 
which  compose  the  proposed  arsenal.  Of  these  five  systems, 
the  first  four  will  be  varied  in  the  study  and,  therefore, 
included  in  the  regression  analysis.  The  fifth  system  is 
constant  and  will  be  shown  in  the  regression  section  to  be 
in  the  Bq  term  of  the  polynomial.  Similarly,  all  the  other 
variables  in  the  LP  which  remain  constant  do  not  have  to 
be  included  in  the  control  variables.  This  selection  of 
control  variables  thus  assumes  that  the  target  base  will 
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remain  constant,  otherwise  the  change  in  target  base  would 
likewise  have  to  be  considered  in  the  regression  analysis 
and  also  in  the  sampling  from  the  LP  model.  This  variable 
selection,  then,  leads  to  the  four  weapons  as  control 
variables  and  the  two  measures  of  merit,  CF  and  CV. 

Control  Variable  Limits 

The  range  of  each  control  variable  is  determined  by 
the  scenario  or  the  study  under  consideration.  For  this 
example,  weapon  five  is  introduced  and  maintained  at  a  con¬ 
stant  level  of  60  weapons  or  360  warheads.  The  remaining 
four  weapons  are  assumed  to  have  previously  determined 
upper  limits  shown  in  Table  I.  Since  all  but  one  of  the 
weapon  systems  has  multiple  warheads,  the  w^  variables  are 
in  terms  of  warheads  rather  than  actual  weapon  numbers. 

The  lower  bounds  for  these  variables  are  dependent  on  the 
range/degree  of  conflict  to  be  considered.  Initially,  the 
lower  limits  are  set  at  zero.  However,  after  the  regression 
error  was  determined  in  step  eight,  the  lower  limits  were 
re-evaluated  and  set  as  shown  in  Table  III.  This  revision 
is  further  explained  in  step  eight.  Regression  Error. 

Data  Source/Generation 

The  data  for  this  generic  problem  are  generated  from 
the  aggregated  arsenal  exchange  model,  which  is  detailed  in 
Appendix  A.  The  choice  of  this  particular  model  was  not 


Weapon  Ranges  (Warheads) 

Initial 

Revised 

0  <  W  <  450 

300  <  W  <  450 

0  <  w2  -  750 

450  <  W2  <  750 

0  <  W_  <  1040 

720  <  W3  <  1040 

0  <  W.  <  424 

200  <  W4  <  424 

—  4  — 

W-  =  360 

Wc  =  360 
o 

critical  to  the  methodology  example;  therefore,  a  simplified 
version  was  chosen. 

The  model  attempts  to  maximize  the  DE  from  the  weapons 
available  on  the  given  target  data  base.  The  separate 
missions  are  modeled  by  using  either  the  CF  or  the  CV 
components  in  the  objective  function.  Thus,  each  mission 
has  the  same  LP  tableau  except  for  the  change  in  objective 
function.  The  two  measures  of  effectiveness,  CV  and  CF, 
are  determined  by  the  two  different  objective  functions 
respectively.  Since  all  other  LP  components  remain  con¬ 
stant,  the  changes  in  the  respective  MOE  are  due  solely  to 
changes  in  the  right-hand-side  control  variables.  Thus, 


the  two  MOE  variables  are  dependent  upon  the  combinations 
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of  independent  variables,  wi ,  levels.  The  function  of  the 
independent  variables,  w^,  is  the  simplex  algorithm,  the 
objective  value  of  which  is  the  MOE  in  question. 


Initial  Surface  Estimation 

As  mentioned  earlier,  the  initial  surface  estimation 
is  based  on  the  study  of  Smith  and  Mellichamp  (30).  In 
that  study,  they  found  the  DE  to  be  a  monotonically  in¬ 
creasing  function  of  the  force  structure.  Thus,  the  study 
assumes  that  the  contour  of  each  surface,  CF  and  CV,  will 
be  similar  to  that  of  the  combined  mission  which  Smith  and 
Mellichamp  studied.  In  light  of  this,  a  second  order  re¬ 
gression  model  is  assumed. 


Research  Design  Selection 

Chapter  Two  shows  the  criteria  for  selection  of  a  re¬ 
search  design.  In  particular,  Smith  and  Mellichamp 's 
analysis  (30)  relies  on  independent  coefficients  which 
require  an  orthogonal  design.  A  minimum  bias  design  is 
optimal  for  the  deterministic  model  but,  as  previously 
mentioned,  there  are  none  available  for  greater  than  two 
dimensions.  Thus,  the  criteria  of  rotatability ,  equispacing, 
and  orthogonality  were  chosen  to  minimize  the  variance  and 
bias.  To  fulfill  these  criteria,  the  Box  and  Behnken's 
fractional  factorial  3k  design  was  chosen  (2).  The  original 
design  calls  for  27  samples.  This  design  was  modified  to  25 
samples  because  the  model  is  deterministic.  The  difference 
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TABLE  IV 


Research  Design  and  Sampled  Values 


Design  Points  (Coded)  CF  CV  AF 


450  ( 

1) 

750  ( 

1) 

880  ( 

0) 

312( 

0) 

932.67 

1190.02 

1736.46 

450  ( 

1) 

450  ( 

-1) 

880  ( 

0) 

31 2  ( 

0) 

368.40 

1159.59 

1618.73 

300  ( 

-1) 

750  ( 

1) 

880( 

0) 

3 1 2  ( 

0) 

362.12 

1170.66 

1652.86 

300  ( 

-1) 

450  ( 

-1) 

88Q( 

0) 

3 1 2  ( 

0) 

797.85 

1139.65 

1504.04 

375  ( 

0) 

600  ( 

0) 

1040( 

1) 

424  ( 

1) 

873.84 

1165.56 

1677.36 

375  ( 

0) 

600  ( 

0) 

1040  ( 

i) 

200  ( 

-1) 

868.13 

1165.56 

1645.10 

375  ( 

0) 

600( 

0) 

720  ( 

-1) 

424  ( 

1) 

865.41 

1164.70 

1626.09 

375  ( 

0) 

6G0( 

0) 

720( 

-1) 

200  ( 

-1) 

859.61 

1164.31 

1585.47 

450  ( 

1) 

600( 

0) 

880  ( 

0) 

424  ( 

1) 

904.92 

1174.87 

1701.47 

450  ( 

1) 

600  ( 

0) 

880  ( 

0) 

200  ( 

-1) 

899.15 

1174.81 

1668.05 

300  ( 

-1) 

600  ( 

0) 

880( 

0) 

424  ( 

1) 

834.40 

1153.39 

1602.27 

300  ( 

-1) 

600  ( 

0) 

880  ( 

0) 

200  ( 

-1) 

828.60 

1155.32 

1567.22 

375  ( 

0) 

750  ( 

1) 

1040  ( 

1) 

312( 

0) 

901.59 

1180.33 

1717.52 

375  ( 

0) 

750  ( 

1) 

720  ( ■ 

-1) 

31  2  ( 

0) 

893.14 

1179.81 

1673.61 

375  ( 

0) 

450  ( 

-1) 

1040  ( 

1) 

312( 

0) 

337.39 

1150.28 

1594.91 

375( 

0) 

450  ( 

-1) 

720( 

-1) 

31  2  ( 

0) 

828.87 

1140.42 

1526.83 

450( 

1) 

600  ( 

0) 

1040  { 

1) 

3 1  2  ( 

0) 

906.24 

1175.30 

1708.49 

450( 

1) 

600  ( 

0) 

720(  • 

-1) 

3 1 2  < 

0) 

897.79 

1174.28 

1657.15 

300  ( 

-1) 

600  ( 

0) 

1040( 

i) 

312  ( 

0) 

835.76 

1155.32 

1612.06 

300(  ■ 

-1) 

600( 

0) 

720 ( ■ 

-1) 

3 1 2  ( 

0) 

827.24 

1153.96 

1554.73 

375  ( 

0) 

750( 

1) 

880( 

0) 

424( 

1) 

900.26 

1180  40 

1710.86 

375( 

0) 

750( 

1) 

880  ( 

0) 

200  ( • 

-i) 

894.50 

1180.2 

681.32 

375( 

0) 

450(- 

-1) 

880  ( 

0) 

424  ( 

1) 

836.03 

1149.85 

585.12 

375( 

0) 

450(  ■ 

-1) 

880( 

0) 

200  ( ■ 

-1 ) 

830.23 

1149.79 

i.  >40.05 

375( 

0) 

600( 

0) 

880( 

0) 

3 1 2  ( 

0) 

877.77 

1165.13 

1  35.88 

877.77  1165.13  1  35.88 


is  in  the  repetition  of  the  center  point  which  is  required 
for  degrees  of  freedom  in  an  experiment  which  has  no  vari¬ 
ance.  Table  IV  displays  this  design  for  the  revised  limits 
in  both  its  coded  and  uncoded  version,  as  described  in 
Chapter  Two.  The  coded  values  are  shown  in  parentheses 
next  to  their  uncoded  counterparts.  This  table  also  shows 
the  sample  values  of  each  MOE  (CV,  CF,  AF )  for  these  re¬ 
spective  design  points.  A  further  description  of  the  design 
is  in  Appendix  B. 


Global  Regression  Estimate 

The  above-mentioned  design  yields  25  data  points  for 
use  in  a  regression  analysis.  Each  regression  looks  at  the 
HOE  regressed  on  the  four  control  variables,  the  cross 
products  of  these  variables  and  the  square  terms  of  these 
variables,  Eq  (15). 


MOE  (wi)  =  B0  +  B1W1  +  B2W2  +  B3W3  +  B^  + 

+  B22W22  +  B33w32  +  B44W42  +  B12W1  '  W2 
+  B  W  •  W  ♦  Bm»,  •  W4  +  B23W2  •  W3 


+  B24W2  *  W4  +  B34W3 


W. 


The  computer  regression  library  chosen  was  the 
Statistical  Package  for  the  Social  Sciences  (SPSS)  with  a 
stepwise  option  (25:330).  An  example  of  a  full  regression 
is  shown  in  Appendix  C.  Likewise,  the  final  steps  of  each 


of  the  other  regressions  is  also  given  in  this  Appendix. 
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The  choice  of  the  stepwise  option  allows  the  researcher  to 
follow  the  introduction  of  each  coefficient  and  display  its 
independence  from  the  other  coefficients.  Eqs  (16),  (17), 
and  (18)  are  the  respective  regression  polynomials  using 
the  initial  limits  shown  in  Table  III.  When  applying  the 
revised  limits'  design  points,  the  regression  polynomials 
are  found  to  be  Eqs  (19),  (20),  and  (21)  respectively. 

AF  =  302.022  +  0.91918W1  +  0.94721W2  +  0.87583W3 
+  0.82253W4  -  0.00008W12  -  0.00014W22  -0.00027W3 

-  0.00018W42  -  0.00011W1W2  -  0.00002W1W3  -  0.00002W1W4 

-  0.00019W2W3  -  0.00022W2W4  -  0.00049W3W*4 

CF  =  370.603  +  0.7024W1  +  0.4473W2  +  0.1813W3  +  0.1604W4 

-  0.00016W12  -  0.000135W22  -  0.00007W32  -  0.00003W42 

-  0.00020W1W2  -  0.00002W1W3  -  0.00004W1W4  -  0.00003W2W3 

-  0.00004W2W4  -  0.00013W3W4 

CV  =  279.659  +  1.2307W1  +  1.0376W2  +  0.5759W3  +  0.7330W4 

-  0 . 00035W1 2  -  0.00035W22  -  0.00016W32  -  0.00036W42 

-  0.00070W1W2  -  0 . 00046W^W3  -  0.00037W1W4  -  0.00036W2W3 

-  0.00026W2W4  -  0.00040W3W4 
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(16) 


(17) 


(18) 


AF  =  323.747508  +  1.357119W1  +  1.335273W2  +  0.5225785W 

+  0.435G739W4  -  0.00020551^^  -  0.0002804W22  -  0.0000657W3 

-  0.0000239W42  -  0.0006909W1W2  -  0.0001248W1W 

-  0.0000485W1W4  -  0.0002518W2W3  -  0 . 000231 1W2W4 

-  0.0001166W3W4 

CF  =  504.542366  +  0.4715619W1  +  0.2955107W2  +  0.0288839W3 
+  0.0278609W4  +  0.00000031^  -  0.0000671W22  -  0.0000006W3 

-  0.0000005W42  -  0. 000001  5WjW3  -  0 .  0000009W.W 

-  0.0000007W2W3  -  0.0000006W2W4  -  0. 00000 13W3W 

CV  =  1027.039  +  0.1658222W1  +  0.1198222W2  +  0.0307574W3 

+  0.0019064W4  +  0000054W42  -  0.0000153W12  -  0.0000040W22 

-  0.0000075W32  -  0.0000129W1W2  -  0.0000175W1W3 

-  0.0000087W2W3  -  0.0000054W3W4 

The  first  SPSS  regression  of  the  revised  limits  rounded 
the  coefficients  to  the  fourth  decimal  place.  Thus,  it 
dropped  most  of  the  second  order  terms  from  the  CF  equation 
and  all  of  the  second  order  and  cross  terms  from  the  CV 
equation.  A  second  regression  on  the  revised  samples  was 
run  to  allow  the  seven  decimal  places  in  Eqs  (19),  (20),  and 
(21).  These  small  coefficients  are  probably  due  to  the 
smaller  region  and  slow  change  of  surface  contour  in  this 
region.  The  next  section  considers  the  goodness  of  fit  of 


these  polynomials  to  the  actual  surfaces  and  explains  the 
reason  for  revised  area  of  consideration. 


Regression  Error 

At  this  point  in  the  research,  it  was  determined  that 
the  initial  regression  polynomials  lacked  the  degree  of 
accuracy  sought  in  the  solution.  This  can  be  seen  by  look¬ 
ing  at  the  standard  deviation  of  the  residual  error  of  each 
regression.  The  standard  deviations  were  17.4,  7.9,  and 
20.3  for  all  force,  counterforce,  and  countervalue  re¬ 
spectively.  Since  this  error  was  thought  to  be  unacceptable 
further  error  random  sampling  was  not  accomplished. 

In  accordance  with  the  methodology,  the  study  returned 
to  steps  six  and  three  of  the  methodology.  In  step  six. 
Research  Design  Selection,  the  Box  and  Behnken  design  was 
reviewed  for  its  orthogonality  and  all  data  points  were  re¬ 
checked  as  to  accuracy.  This  review  showed  the  area  of 
consideration  to  be  too  large  for  the  second  order  response 
surface.  Thus,  a  redefinition  of  the  problem  scenario  or 
level  cf  conflict  was  considered.  The  redefinition  found 
the  more  realistic  ranging  of  variables  listed  in  Table  IV. 
Therefore,  the  revised  lower  limits  were  applied  in  step 
six,  a  new  sampling  taken,  and  the  methodology  proceeded 
as  designed. 

The  error  for  the  second  set  of  regression  equations, 
Eqs  (19),  (20),  and  (21),  was  found  to  be  much  smaller. 

Table  V  summarizes  the  regression  error  for  the  CF ,  CY,  and 


AF  surfaces  respectively.  The  table  shows  the  high  degree 
of  accuracy  with  the  revised  control  variable  limits.  The 
CF  and  CV  estimates  are  accurate  to  within  a  single  unit 
for  even  three  standard  deviations  of  error.  This  small 
error  is  significant  in  the  next  section. 

TABLE  V 

Regression  Error  Summary 


Surface 

Design  S.D. 

Of  Error 

Random  Mean 
Error 

Random  S.D. 
of  Error 

CF 

0.0131 

0.223 

0.204 

CV 

0.1569 

0.234 

0.222 

AF 

1.0066 

1.233 

1.154 

Contour  Interval  Determination 


Chapter  Two  determined  the  criteria  for  minimum  in¬ 
terval  to  be  no  closer  than  two  standard  deviations  of 
error.  From  the  above-mentioned  regression  error,  the 
minimum  interval  for  counterforce  and  countervalue  con¬ 
strained  surfaces  was  determined  to  be  no  less  than  one  unit 
In  the  case  of  AF,  it  was  determined  that  it  would  be  no 
less  than  three  units.  In  actual  fact,  it  was  not  neces¬ 
sary  to  utilize  such  small  steps  to  show  the  optimal 
trade-off  between  surfaces.  A  step  of  two  units  was  more 
than  sufficient  for  counterforce  and  countervalue  and  five 
units  for  all  force. 


Constrained  Maxima  Determination 

As  described  in  the  previous  chapter,  the  SUMT  program 
was  used  to  optimize  and  find  the  maximum  of  one  surface 
constrained  to  the  contours  of  another.  For  the  planned 
analysis  in  the  next  chapter,  three  cases  of  optimization 
are  looked  at.  The  first  two  cases  consider  the  optimiza¬ 
tion  of  CF  given  a  value  of  CV  and  CV  given  a  value  of  CF. 
Both  these  cases  use  simplified,  regression  polynomials 
described  in  the  next  chapter.  The  third  case  compares 
the  difference  between  the  use  of  simplified  and  detailed 
regression  polynomials.  The  data  for  these  cases  are  con¬ 
tained  in  Appendix  D. 

Conclusion 

This  chapter  has  applied  the  proposed  methodology  to 
a  generic,  military  force  structure  problem.  The  final 
results  are  compiled  in  Appendices  D  and  E.  The  next 
chapter  provides  an  analysis  of  these  results. 


CHAPTER  FOUR 


Analysis 


Introduction 

This  chapter  presents  the  study's  analysis  utilizing 
Chapter  Three's  example  data.  The  analysis  goals  are  twofold 
One  is  to  develop  insight  into  the  optimization  process  be¬ 
tween  the  measures  of  effectiveness  (MOE's)  and  their  re¬ 
spective  domain  elements  (n-tuples)  in  achieving  the  con¬ 
strained  maxima.  The  second  goal  is  to  offer  the  decision 
maker  this  analysis  in  a  usable  format.  Prior  to  the 
analysis,  redefinition  of  Chapter  Three's  works  and  some 
simplifying  assumptions  are  made. 

In  Chapter  Three,  the  MOE  response  surface  is  generated 
from  a  linear  programming  optimization  model.  Thus,  each 
point  on  the  surface  represents  the  optimal  use  of  its  re¬ 
spective  domain  elements  (n-tuple).  Since  each  point  is  an 
optimal,  the  surface  displays  the  optimal  MOE  for  the 
system's  entire  range  of  domain  elements.  Therefore,  the 
regression  polynomials  in  Chapter  Three  (Eqs  (19),  (20),  and 
(21)  )  represent  the  regression  estimates  for  the  optimal 
employment  of  the  domain  elements  to  achieve  their  re¬ 
spective  MOE.  In  terms  of  the  example,  each  equation  re¬ 
presents  the  optimal  force  structure  to  achieve  its  re¬ 
spective  mission  (MOE).  Using  the  above-mentioned  poly¬ 
nomial  definitions,  the  SUMT  program  selects  from  a  domain 
subset  (force),  already  optimally  designed  to  fulfill  one 


MOE  (mission),  the  best  mixture  of  domain  elements  (weapon 
systems, w^)  to  fulfill  a  second  MOE  (mission).  In  the 
analysis,  the  first  goal  is  to  determine  the  cause  for  the 
optimal  trade-off  between  these  MOE's  and  the  domain  n- 
tuples  which  generate  this  trade-off.  To  aid  this  discus¬ 
sion  of  constrained  maxima  (optimal  trade-off),  the  regres¬ 
sion  polynomials  are  simplified  by  rounding  the  coefficients 
to  four  decimals. 

These  simplified  polynomials,  Eqs  (19a),  (20a),  and 
(21a),  display  three  major  points  of  discussion  in  the 
chapter. 


=  323.7475  + 


i.357i*w1  -  0.0002*1^  -  o.ooo7*w1v;2 


-  0.0001*W1W3  +  1.3353*W2  -  0.0003*W2 

-  0.0003*W2W3  -  0.0002*W2W4  +  0.5226*W3 

-  0.0001*W  2  -  0.0001*W.,W .  +  0.4351*W„ 

3  3  4  4 


CF  =  504.542  +  0.4716*1^  +  0.2955*W2  -  0.301*W2 

+  0.0289*W3  +  0.0279*W4 


CV  =  1027.039  +  0.1658*W1  +  0.1198*W^  +  0.0308*W3 

+  0 . 0019*W„ 


First,  they  are  used  to  show  a  graphic  representation  of 
the  optimization  process  of  an  MOE  surface  constrained  to 
a  single  contour?  thus,  giving  some  pictorial  validity 


v/.v.  .v  .*  .*  .*  , 


by  demonstrating  graphically  the  results  of  SUMT.  Second, 
the  simplified  regression  equations  aid  in  the  analysis  of 
the  MOE's  optimal  trade-off  and  the  priority  process  for 
entering  domain  variables.  Finally,  it  is  hoped  that  while 
these  equations  serve  to  improve  the  quality  of  discussion, 
they  also  serve  to  exemplify  possible  pitfalls  in  rounding 
the  regression  coefficients. 

The  chapter  presentation  begins  with  the  above- 
mentioned  graphic  representation  of  the  optimizaton  process. 
This  is  accomplished  through  a  series  of  MOE  contour  plots, 
one  type  MOE  overlaid  upon  another  type  MOE.  Following  the 
graphic  optimization,  the  data  from  Chapter  Three  is  mapped 
and  analyzed  over  the  range  of  domain  elements.  The 
chapter  concludes  with  a  summary  of  the  methodology's  ad¬ 
vantages  to  the  decision  maker. 

Graphic  Optimization 

This  section  graphically  displays  the  SUMT  optimization 
of  CF  for  a  given  CV  contour.  The  simplified  regression 
polynomials,  Eqs  (20a)  and  (21a),  are  used  in  this  section 
of  the  analysis. 

In  Chapter  Three,  CF  was  optimized  for  CV  values  of  1154 
to  1224  in  two  unit  steps.  These  CV  lower  and  upper  limits 
are  determined  by  the  lower  and  upper  limits  of  the  domain 
elements.  The  standard  deviation  of  error  for  CV  was  found 
to  be  0.222.  Thus,  the  minimal  interval  criteria  was  evalu- 


ted  to  be  one  unit.  This  close  an  interval  was  not 


necessary  to  adequately  represent  the  trade-off  between  CF 
and  CV.  Thus,  the  contour  interval  of  two  units  was  chosen 
Appendix  E  contains  a  full  listing  of  the  example  MOE  com¬ 
parisons  and  their  respective  domain  elements.  A  review  of 
these  data  shows  is  constant  at  450  in  the  region  of 
CV  =  1176.  Therefore,  the  graphic  example  uses  CV  =  1176 
and  Wg  =  450. 

In  this  experimental  region  (CV  =  1176,  =  450),  the 

contour  plots  consist  of  variables  ,  Wg ,  W  ,  and  the  MOE. 
This  requires  representation  of  four  dimensions  on  a  two- 
dimensional  plot.  This  is  accomplished  through  multiple 
contours . 

The  comparison  data  plots  for  CF  and  CV  are  generated 
by  a  computer  program  which  iterates  each  feasible  combina¬ 
tion  of  control  variables  for  each  regression  polynomial. 
These  data  were  then  plotted  and  are  availabl  in  Appendix  E 
Several  of  the  plots  in  Appendix  E  have  been  reduced  in 
scale  and  are  used  in  combined  overlays  throughout  the  re¬ 
mainder  of  the  chapter. 

Figure  4  represents  the  counterforce  contour  value  of 
85C  and  also  860  (dashed  line)  with  ordinate  W  and  ab¬ 
scissa  Wg.  Each  diagonal  line  across  the  graph  is  a  con¬ 
tour  of  both  CF  and  weapon  one.  As  such,  each  is  identified 
by  the  pair  of  numbers  in  parentheses.  The  first  is  the 


MOE  and  the  second  is  the  W,  value. 
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Figure  5  represents  the  countervalue  of  1176  with  the 
various  diagonal  lines  due  to  the  stepwise  increase  in  W^. 
This  representation  of  five-space  in  two-dimensional  plots 
is  similar  to  that  used  by  Smith  and  Mellichamp  (30). 

If  Figures  4  and  5  are  combined,  the  intersection  of 
common  value  contours  reveals  the  area  of  common  domain 
between  surfaces.  Thus,  the  graphic  comparison  optimiza¬ 
tion  becomes  a  search  for  the  maximum  CF  value  with  a 
common  W-^  value  on  both  graphs.  Figure  6  is  the  overlay  of 
Figures  4  and  5. 


The  figure  displays  a  commonality  in  the  1176  CV  contour  and 
the  850  CF  contour;  however,  there  is  no  common  domain  point (s) 
for  the  860  CF  contours  and  the  contours  of  1176  CV.  There¬ 
fore,  Figure  6  has  bounded  the  graphic  optimization  search 
between  CF  values  of  850  and  860. 

Figures  7  and  8  have  similarly  overlaid  the  CV  plot  of 
1176  against  the  CF  plots  of  854  and  856  respectively. 

Figure  7  shows  the  common  points  for  values  of  at  435 
and  430.  Figure  8  has  no  common  points.  However,  Figure  8 
does  show  the  closest  point  to  a  common  W,  value  of  435. 

This  point  is  just  off  the  W.  axis  and  is  marked  "Z." 


Fig  8.  Overlay  Plot  of  CV  =  1176  and  CF  =  856 


Graphically,  the  optimal  point  lies  between  the  CF 
values  of  854  and  856,  and  the  value  lies  between  430 
and  435.  This  graphic  solution  compares  very  well  with 
that  predicted  by  the  SUMT  optimization  in  Appendix  E. 

The  SUMT  operating  domain  is  434.9,  450,  720,  and  400.3  for 
w:  through  respectively.  The  graphic  solutions  between 
Figures  7  and  8  show  this  to  be  the  operating  domain  upon 
which  the  plots  converge.  Likewise,  SUMT's  predicted 
maximum  CF  of  855.36  very  closely  approximates  the  "Z" 
point  just  off  the  axis  of  the  CF  856  plot. 


Therefore,  the  overlay  contour  plots  in  Appendix  E 
have  graphically  validated  the  optimization  process  used 
in  the  SUMT  algorithm.  This  plot  overlay  technique  will 
be  applied  again  in  the  next  sections. 

Data  Presentation 

At  the  conclusion  of  Chapter  Three,  the  data  had  yet 
to  be  presented.  This  section  displays  the  data  mapping 
that  will  be  used  in  the  following  section  for  the  analysis 
The  raw  data  are  presented  in  Appendix  E  and  are  presented 
here,  graphically,  in  Figures  9  and  10.  These  graphs  are 
shown  in  larger  scale  in  Appendix  E. 

These  graphs  and  most  of  the  analysis  are  based  on  the 
simplified  regression  polynomials,  Eqs  (20a)  and  (21a).  Th 
rounding  error  for  the  terms  less  than  four  decimal  places 
is  initially  assumed  small  and  simply  added  into  the  re¬ 
gression  error.  The  danger  of  this  assumption  is  discussed 
in  the  closing  of  the  analysis. 

Figure  9  represents  the  maximum  countervalue  available 
given  a  counterforce  value  requirement.  Counterforce  is 
displayed  on  the  left,  vertical  axis,  and  countervalue  on 
the  horizontal  axis.  This  graph  also  displays  the  optimal 
weapons  combination  to  achieve  this  maximum  CF  for  a  given 
CV.  The  number  of  weapons  is  given  along  the  right- 
vertical  axis.  For  example,  given  a  CV  of  1190,  the  maxi¬ 
mum  CF  is  approximately  880.  Likewise,  the  respective 
weapons  set  is  450,  536,  757,  and  393.  At  the  beginning  of 
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the  graph,  weapons  two  and  three  are  approximately  constant 
at  their  lower  limits  of  450  and  720  respectively.  Weapon 
one  is  steadily  increasing  to  its  maximum  limit  of  450. 
Weapon  one  reaches  this  upper  limit  at  approximately  1180 
CV .  Interpretation  of  this  figure  is  left  for  the  analysis 
in  the  next  section.  Figure  10  is  similar  to  Figure  9,  but 
the  roles  of  CV  and  CF  are  reversed.  Thus,  the  graph  shows 
the  maximum  CV  available  for  a  given  CF.  In  this  case, 
weapon  one  remains  constant  at  the  beginning  of  the  graph 
while  weapon  two  is  steadily  increasing  towards  its  upper 
limit.  At  approximately  850  CF,  weapon  one  begins  increas¬ 
ing  as  weapon  two  has  reached  its  limit.  Both  these  figures 
are  compared  and  discussed  in  the  following  section. 


Analysis 

This  section  is  concerned  with  three  areas.  First,  an 
explanation  of  the  trade-off  between  MOE's  and  what  causes 
the  changes  in  this  trade-off.  Second,  the  determination 
of  the  causes  for  the  force  structure  changes  and,  in 
particular,  the  apparent  randomness  of  W^  and  W^  values  on 
the  graphs.  Finally,  the  danger  involved  in  using  the 
rounded-off,  simplified  regression  equations  rather  than 
the  full  polynomials  is  discussed. 

Figure  9  shows  the  relative  trade-off  between  maxi¬ 
mum  CF  and  CV  as  approximately  two  regions.  Both  these 
regions  are  approximately  linear.  The  first  region 
runs  from  a  CV  value  of  approximately  1154  to  1178. 
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The  trade-off  between  MOE's  in  this  region  is  approximated 


by  a  straight-line  function  of  slope,  2.9.  Likewise,  the 
second  region  runs  approximately  from  1180  to  1206  CY  and 
is  closely  approximated  by  a  straight-line  of  slope,  1.52. 
These  approximate  slopes  may  be  accurately  calculated  by  the 
ratio  of  the  total  derivatives  of  each  of  the  regression 
polynomials,  Eqs  (20a)  and  (21a).  These  total  derivatives 
are  shown  in  Eqs  (22)  and  (23)  with  their  respective  examples 
in  Eqs  (22a)  and  (23a). 

dCF  =  («CF/aw1)*dw1  +  (aCF/aw2)*dw  +  (acF/aW-)*dW- 

+  (aCF/a  w4)*dw4 

dCF  =  (.4716)*dW1  +  ((. 2955 )-(. 002 ) *W2 ) *dW2  +  (.0289)*dW3 

+  (.0279)*dW4 


dev  =  (acv/aw.  )*dw. 

.  i  i 


dCV  =  (.1658)*dW1  +  (.1198)*dW2  +  (.0308)*dW3  +  (,0019)*dW4 


Thus,  the  trade-off  between  MOE's  is  equal  to  the  slope, 
dM0E1/dM0E2. 

From  Chapter  Three's  data  (Appendix  E),  the  ratio  of 
CF  to  CV  changes,  AC F/ACV,  may  be  compared  to  the  total 
derivative's  ratio.  For  example,  from  a  CV  of  1162  to  a 
CV  of  1168,  the  ACF/ACV  =  2.8817,  which  compares  to  the 


ratio  of  total  derivatives  of  CF  to  CV  of  2.8766, 

(dCF/dCV  =  2.8766).  Likewise,  the  second  region  from  CV 
1180  to  1206  has  a  slope,  ACF/ACV ,  equal  to  1.53,  which  com¬ 
pares  to  the  total  derivative  ratio  of  1.51. 

Since  the  change  in  weapons  (dVL  )  is  the  same  for  both 
the  CF  and  CV  total  derivatives,  the  ratio  of  each  co¬ 
efficient  may  be  compared  separately  if  the  other  dW ■ ' s  are 
relatively  small.  Thus,  in  this  case,  the  dW^  coefficient 
equals  .4716/. 1658  which  equals  2.84.  In  the  region  1162  CV 
to  1168  CV  (Figure  9),  only  is  actively  changing.  Thus, 
the  trade-off  between  MOE 1 s  in  this  region  is  closely  ap¬ 
proximated  by  the  dW^  coefficients'  ratio.  The  difference 
between  the  two  ratios  (2.8766-2.84  =  .0366)  is  due  to  the 
small  change  in  the  other  weapon  systems  (dW^). 

This  technique,  of  course,  becomes  more  difficult  when 
the  coefficients  have  cross  terms  or  square  terms  such  as 
the  W 2  term.  In  this  case,  an  average  value  of  W 2  for  the 
range  of  consideration  must  be  applied  in  the  calculations. 
Although  this  explains  the  relative  trade-off  or  slope  of 
the  maximum  CF  given  a  CV  line,  it  should  not  be  inferred 
as  the  effectiveness  of  a  weapon  system  to  increase  the 
maximizing  MOE.  Weapon  four,  for  example,  has  a  ratio  of 
14.68  for  its  MOE  trade-off.  However,  the  actual  change 
in  CF  per  weapon  four  is  very  small  with  its  CF  regression 
coefficient  of  .0279.  Consequently,  although  the  relative 
ratio  may  give  a  high  rate  of  return  in  terms  of  one  MOE  to 


another,  it  may  have  little  importance  as  to  the  actual 
change  in  total  value  of  the  maximum  MOE.  This  leads  to 
the  question  of  which  weapons  are  the  most  effective  in  the 
optimal  force  structure  to  maximize  a  secondary  MOE  given 
the  level  of  a  primary  constraining  MOE. 

Smith  and  Mellichamp's  work  (30)  shows  the  independent 
regression  polynomial  coefficients  used  to  determine  the 
incremental  increase  in  the  MOE  per  weapon.  Using  this 
rationale,  the  maximizing  MOE  prefers  the  weapon  system 
with  the  highest  coefficient  first.  Thus,  the  increase  in 
the  MOE  will  be  due  to  the  maximizing  MOE ' s  regression  co¬ 
efficients  and  the  amount  of  weapons  available.  The  number 
of  weapons  available,  however,  becomes  a  function  of  the 
level  of  the  constraining  MOE.  Thus,  given  a  CV  level, 
the  operating  domain  is  restricted  as  shown  graphically 
with  the  contour  lines.  Operating  within  this  restricted 
domain,  the  maximizing  function  attempts  to  trade-off  the 
highest  coefficient  weapons  system  with  the  number  available 
and  compare  its  relative  MOE  increase  to  that  of  introducing 
some  less  efficient  system  that  may  change  the  number  of 
weapons  available.  This  is  best  shown  in  Eq  (24). 

CF(W. )  =  eCF  *  [Number  of  W.  available] 

1  aw. 


CF ( W, )  =  ( .4716)* 


f(CV-(  .  1198)*W0-(  .0308)*W3-(  .0019  )  *W4  )  "j 


In  Eq  (24)  the  number  of  available  is  simply  the 
regression  polynomial  for  the  constraining  MOE  solved  for 
that  particular  weapon  system  (V?^).  Eq  (25)  shows  the 
example  of  maximizing  CF  in  terms  of  weapon  system  one.  The 
maximizing  algorithm  (SUMT)  must  trade-off  the  introducton 
of  a  lesser  efficient  weapon  system  against  the  value  of 
not  decreasing  the  number  of  weapons  available  for  con¬ 
tribution  in  Eq  (25).  This  is  numerically  shown  in  Eqs  (26) 

v  -.7  )  ,  and  (  28  )  . 

CF(W^)  gain  >  CF(W. )  gain 
.4716  ( ‘ X658  W3)  =  *188W3  >  .0289W3 

•471e  (TlMI  V  =  -0U6W4  ^  ,0279W4 

When  Eq  (26)  is  true,  then  the  algorithm  will  not  increase 
Wi.  Eq  (27)  is  the  gain  to  CF  as  a  function  of  weapon  one 
due  to  the  level  of  weapon  three  as  compared  to  the  value 
of  weapon  three,  as  deterined  by  its  CF  regression  co¬ 
efficient.  In  Eq  (27)  a  .188  gain  is  determined  in  the  CF 
equation  due  to  weapon  one,  whereas  only  .0289  is  due  to 
weapon  three.  This  case  is  not  true,  however,  in  Eq  (28) 
where  the  gain  for  weapon  one  is  only  .0116  as  compared  to 
the  actual  increase  of  .0279  if  weapon  four  is  introduced. 
Therefore,  the  algorithm  will  choose  to  enter  weapon  four 
rather  than  leave  it  out  in  favor  of  the  number  of  weapon 
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ones  available.  This  preference  for  W.  while  W,  is  in- 

4  i 

creasing  is  shown  in  Figure  9.  In  the  increasing  region, 
W4  is  at  or  near  its  upper  limits  while  is  at  its  lower 
limit. 

It  is  due  to  this  algorithm  trade-off  comparison  be¬ 
tween  each  domain  system's  relative  MOE  increase  that  the 
graph  of  the  data  appears  erratic.  This  is  a  especially 
true  for  the  smaller  coefficient  values  such  as  weapons 
three  and  four.  In  light  of  this  erratic  change,  no  attempt 
has  been  made  to  smooth  the  weapons'  curves. 

The  section,  thus  far,  has  explained  the  trade-off 
and  algorithm  optimization  of  the  force  structure  by  using 
the  simplified  regression  polynomials.  The  final  point  in 
this  analysis  section  is  to  determine  the  error  introduced 
by  this  simplification  and  its  implications. 

Figure  11  is  a  plot  of  max  CF  vs.  CV  using  the  full 
polynomials  given  in  Eqs  (20)  and  (21).  A  comparison  of 
Figures  11  and  9  shows  there  is,  indeed,  a  difference  in 
not  only  trade-off  value  but  also  in  the  weapon  systems' 
order  and  rate  of  entry. 

The  reason  for  the  change  in  slope  or  trade-off  is 
due  to  both  a  reduced  value  in  CV  and  the  inclusion  of  the 
square  and  cross  terms,  and  also  an  increased  value  in  CF 
due  to  this  more  accurate  use  of  the  square  and  cross 
terms.  At  first  glance,  this  does  not  appear  plausible  for 
the  case  of  CF.  However,  the  rounding  up  of  only  the  K-,  " 


CF  coefficient  causes  the  change  m  maximum  CF  available 
from  the  original  939  down  to  a  lesser  value  of  924  in  the 
simplified  equation.  Thus,  the  rate  of  trade-off  is  less  in 
a  simplified  case. 

The  change  in  the  force  structure  due  to  the  rate  at 
which  different  systems  are  entering  is  not  greatly  affected 
by  weapons  one  or  two.  However,  there  is  a  large  change  in 
weapon  three.  Again,  this  change  j s  due  to  the  trade-off 
done  by  the  optimization  technique  when  looking  at  weapon 
three's  introduction  in  comparison  to  reducing  the  numbers 
of  weapons  one  and  two  available. 

Throughout  the  analysis,  the  only  mention  of  the  bias 
error  listed  in  Table  V  has  been  in  the  minimum  contour  in¬ 
terval  determination.  Ignoring  such  a  small  error  has 
little  effect  on  the  MOE  trade-off  calculation;  however,  it 
may  affect  the  force  structure  numbers  of  the  lesser  effici¬ 
ent  domain  elements  (weapons).  Thus,  the  constraining 
numbers  of  weapons  available  should  be  recalculated  to  de¬ 
termine  if  two  standard  deviations  of  error  significantly 
change  the  optimal  n-tuples. 

In  conclusion,  the  analysis  has  found  and  been  able  to 
describe  the  relationship  between  measures  of  effectiveness 
and  the  optimal  force  structure.  Likewise,  it  has  shown 


that  the  original  assumption  of  rounding  to  the  fourth 
decimal  place  changes  the  overall  structure  of  the  problem. 


The  question  of  rounding  the  four  digits  originally 
arose  in  the  research  because  the  regression  package  used 
automatically  rounded  to  four  decimal  places  as  seen  in 
Appendix  C.  However,  when  the  package  did  compute  its 
residual  error  it  had,  in  fact,  maintained  the  values  of 
those  coefficients  to  the  degree  of  accuracy  that  the 
computing  device  used.  Thus,  rathe.;  than  using  the  SPSS 
automated  covariance  matrix,  the  problem  could  have  been 
avoided  by  using  the  coded  design. 

Decision  Tools 

The  above  methodology  should,  with  a  high  degree  of 
accuracy,  predict  an  optimal  force  to  achieve  a  maximum 
measure  of  effectiveness  given  a  constraint  level  of  a 
second  measure  of  effectiveness.  Thus,  the  decision  maker 
has,  in  Figure  11,  the  ability  to  see  not  only  his  optimal 
force  structure,  but  how  that  force  will  change  as  changes 
occur  in  his  MOE  requirements. 

So  far  this  methodology  has  given  the  decision  maker 
his  best  single  domain  set  for  a  given  constraint  MOE.  It 
is  foolish  to  -assume  that  in  any  analysis  all  the  factors 
would  be  known  or  that  some  which  were  assumed  to  be  constant 
in  the  operating  environment  may  not  have  changed.  In  light 
this, a  choice  of  optimal  solutions  would  probably  be 
'  *•  r  suited  as  a  tool  for  the  decision  maker.  The  SUMT 
•  ;•  finds  only  a  single  optimal  solution;  however. 


the  graphic  techniques  applied  in  the  opening  of  this 
chapter  can  also  show  a  near  optimal  set  of  solutions. 

Figure  12  shows  the  overlay  for  CV  1176  with  CF  plots  of 
850,  852,  and  854  CV.  Lines  connecting  the  common  domain 
points  for  each  CF  value  are  shown  in  Figure  12  as  iso-CF 
lines  (heavy  print  lines),  given  the  CV  level  of  1176. 

Thus,  any  of  the  force  structures  along  one  of  the  iso-CF 
lines  will  guarantee  that  CF  value  and  the  CV  of  1176.  For 
example,  the  line  AB  is  a  CF  value  of  854,  given  a  CV  value 
of  1176.  This  AB  line  contains  domain  ranges  of  lv.  from 
436  to  427,  from  720  to  765,  and  from  305  to  420.  Al¬ 
though  these  alternatives  have  given  up  approximately  1.4 
CF  from  the  SUMT  optimal  solution  of  355.4  CF,  it  may  be 
of  greater  overall  value  to  the  decision  maker  to  have  these 
multiple  options  available  to  him.  Thus,  these  lines  offer 
the  decision  maker  a  trade-off  from  a  single  optimal  solu¬ 
tion  to  a  set  of  multiple  solutions. 

Finally,  the  research  has  yet  to  consider  the  decision 
maker  weighting  the  value  of  one  MOE  to  another.  However, 
if  the  decision  maker  did  decide  to  weight  the  MOE's,  it 
would  have  little  change  on  the  methodology  as  it  has  been 
applied.  For  example,  if  the  decision  maker  decided  that 
CF  was  twice  as  valuable  as  CV,  the  only  change  would  be  to 
multiply  the  total  derivative  ratio  denominator  by  two. 

This  is  shown  in  Eqs  (29a)  and  (29b)  respectively. 


Thus,  the  slope  of  the  trade-off  curve  has  decreased; 
however,  the  optimal  force  to  create  that  curve  is  un¬ 
changed  . 

CF  =  2CV 

dCF  =  dCF 
d ( 2 • CV )  2 • dCV 

If,  on  the  other  hand,  the  decision  maker  were  to 
change  the  value  of  the  individual  targets  rather  than  the 
value  of  the  overall  HOE  class,  then  this  would  require  a 
complete  re-sampling  and  running  of  the  model. 

Conclusion 

This  chapter  contains  the  analysis  portion  of  the  re¬ 
search.  The  analysis  began  with  a  graphic  technique  to 
maximize  the  counterforce  MOE  given  a  countervalue  con¬ 
straint  contour  of  1176.  This  graphic  presentation  vali¬ 
dated  the  optimization  technique  and  the  SUHT  algorithm 
employed  to  calculate  the  constrained  maxima.  Following 
the  single  graphic  case,  the  maximum  CF  trade-off  given  a 
CV  level  was  analyzed  for  the  entire  experimental  domain. 

This  analysis  revealed  the  rate  of  trade-off  to  be  the 
MOE’s  ratio  of  total  derivatives  ( dMOEQ/ dMOEc ) .  This  maxi¬ 
mized  trade-off  was  determined  by  the  optimal  selection  of 
domain  elements  (n-tuples).  This  selection  was  shown  to  be 
a  complex  trade-off  between  each  element’s  ability  to  in¬ 
crease  the  maximized  MOE  (dMOEo/ddi )  and  the  availability 


of  each  element  as  determined  by  the  constraint  MOE 
(di  =  f(MOEc)  ). 

Finally,  the  decision  maker's  use  of  the  analysis  was 
greatly  enhanced  by  utilizing  "near  optimal"  solution  sets. 
These  solutions  were  shown  graphically  to  contain  multiple 
n-tuples  which  achieve  equal  value  MOE  (iso-CF  line,  Figure 
12).  Thus,  for  a  relatively  small  decrease  from  the  SUMT 
optimal  (less  than  one  percent),  the  decision  maker  was 
given  greater  flexibility  in  his  choice  of  operating  domain. 

The  advantages  of  this  analysis  are  further  discussed 
in  the  concluding  chapter. 


CHAPTER  FIVE 


Concluding  Remarks 

Research  Summary 

This  study  presents  the  optimization  methodology  for 
a  deterministic,  mathematically  modeled  system  as  defined 
by  the  aspiration  levels  of  one  measure  of  effectiveness 
(goal/mission)  and  optimization  of  a  second  measure  of  ef¬ 
fectiveness  (MOE).  The  methodology's  advantage  over  other 
multicriteria  decision  theories  is  its  ability  to  consider 
large  numbers  of  alternatives,  show  the  MOE's  comparisons 
throughout  the  range  of  alternatives,  and  use  non-commensur- 
ate  measures  of  effectiveness . 

The  research  example  demonstrates  the  methodology  ap¬ 
plied  to  a  generic,  aggregated  military  force  structure 
problem  defined  by  a  linear  programming  arsenal  exchange 
model.  The  example  has  two  measures  of  effectiveness, 
counterforce  (CF)  and  countervalue  (CV).  The  system  involves 
five  weapon  systems,  of  which  four  are  varied  over  their 
plausible  ranges,  and  ten  target  classes.  The  number  of 
feasible  force  structure  alternatives  over  these  ranges 
exceed  three  billion.  The  optimal  effectiveness  of  these 
alternatives  is  closely  approximated  (less  than  one  percent 
error)  by  use  of  response  surface  methodology. 

The  employment  of  a  central  composite  3  fractional 
factorial  research  design  enabled  the  generation  of  a  sur¬ 
face  for  each  MOE  from  a  minimal  number  of  model  runs.  The 


original  design  required  27  samples  per  surface  (2). 

However,  the  deterministic  model  has  no  variance;  thus,  the 
samples  are  reduced  to  25  per  surface.  With  these  two  sur¬ 
faces,  CF  and  CV ,  the  methodology  then  shows  an  appropriate 
interval  to  describe  the  optimal  trade-off  between  the 
aspiration  MOE  and  the  second  MOE.  The  interval  selection 
criteria  includes  both  the  regression  error  and  the  change 
in  the  contour  surface.  With  this  selection  of  contours, 
the  methodology  then  optimizes  the  second  measure  of  effec¬ 
tiveness  as  constrained  to  the  contour  lines  of  the 
aspiration  measure  of  effectiveness.  These  data  are  then 
plotted  and  displayed  in  Chapter  Four  of  the  analysis  and 
also  in  Appendix  E. 

The  analyses  of  these  data  give  the  decision  maker  a 
range  of  tools  to  gain  further  insight  into  the  optimal 
operation  of  the  system.  The  data  map  displays  the  de¬ 
cision  maker's  optimal  system  structure  given  his  goal 
aspiration  level.  As  new  demands  are  made  upon  the  system's 
measures  of  effectiveness,  this  map  shows  the  required 
changes  to  each  component  of  the  system  to  achieve  the 
higher  MOE  demands.  The  graphic  techniques  of  Chapter  Four 
offer  the  decision  maker  a  choice  of  multiple  solutions 
which  very  closely  approximate  the  optimal  trade-off  value. 
For  the  example  given,  an  aspiration/demand  of  1176  CV,  the 
decision  maker  has  a  wide  range  of  force  structures  to 
choose  from  by  simply  reducing  the  CF  value  by  1.4  from  the 
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optimal  of  855.4  This  reduction  is  less  than  two  tenths 
of  a  percent  from  the  optimal  value.  In  consideration  of 
the  accuracy  of  the  vulnerability  determination  and  single¬ 
shot  probabilities  of  kill  used  in  the  arsenal  exchange 
model,  this  small  reduction  in  the  optimization  process  is 
minimal.  However,  this  negligible  reduction  gives  the 
decision  maker  greater  flexibility  in  his  choice  of  the 
system's  operational  structure. 


Applications 

Considering  the  number  of  deterministic  optimization 
models  in  use  today,  this  methodology  should  have  wide¬ 
spread  application  when  viewing  multiple  missions  or  roles 
of  the  system. 

The  military  applications  are  several.  The  methodology 
application  to  force  structure  planning  is  well  illustrated 
in  the  research  example  and  Smith  and  Mellichamp's  work 
(30).  Beside  force  planning,  the  methodology  may  give  in¬ 
sight  as  a  measurement  and  design  tool  for  weapon  systems. 

Many  of  today's  weapon  systems  are  required  to  perform 
multiple  missions,  some  of  which  may  be  performed  by  several 
configurations  of  a  weapon  system  or  by  different  systems. 
This  methodology  may  serve  as  a  measurement  tool  to  de¬ 
termine  which  configuration  or  system  performs  best  in  the 
dual  mission.  Likewise,  the  methodology  may  be  applied  as 
a  design  tool  for  the  weapon  systems.  For  instance,  if  a 
new  bomber  was  being  designed,  it  may  have  a  strategic 
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penetration  mission  (aspiration  MOE)  and  also  a  search 
and  destroy  mission  (second  MOE).  The  operating  domain  would 
be  the  flight  characteristics  and  capability  of  the  bomber. 
This  methodology  would  choose  the  best  selection  of 
operational  capabilities  to  achieve  the  aspired  penetration 
goal  and  also  maximize  the  second  mission  of  search  and 
destroy.  It  would  also  show  how  these  capabilities  would 
have  to  be  -'hanged  to  meet  new  mission  requirements.  Thus, 
this  methodology  offers  a  planning  tool  for  new  systems  as 
well  as  total  forces. 

A  larger  scale  application  may  be  in  force  structure 
as  seen  in  an  arms  reduction  or  limitations  agreement.  In 
this  case,  if  both  parties  were  able  to  represent  their  re¬ 
spective  capabilities  in  a  response  surface,  then  each 
would  be  able  to  offer  and  counteroffer  an  optimal  selection 
of  weapon  systems  from  their  particular  arsenal  that  equated 
to  equal  levels  of  damage  expectancy  for  the  agreement. 

Further  Research 

This  methodology  was  originally  proposed  for  the  com¬ 
parison  of  non-commensurate  goals.  As  the  study  evolved, 
it  was  obvious  the  scope  would  not  allow  the  comparison  of 
more  than  two  surfaces  at  a  time.  As  such,  further  re¬ 
search  is  required  for  the  comparison  of  three  or  more 
measures  of  effectiveness. 

The  dual  measure  of  merit  case  has  the  advantage  of 
not  changing  the  optimal  solution  domain  regardless  of  the 


decision  maker's  weighting  of  the  MOE's.  However,  in  the 
case  of  three  or  more  measures  of  merit,  this  weighting  may 
determine  the  interval  contour  pattern  as  it  restricts  the 
operating  domain  of  the  maximized  MOE  under  consideration. 

Another  approach  may  be  to  use  the  dual  measure  of 
effectiveness  and  then  project  each  maximized  function  into 
a  next  higher  dimension  and  find  a  third  response  surface 
common  to  the  two  optimization  processes.  For  example, 
find  the  maximum  A  subject  to  C  and  the  maximum  B  subject 
to  C.  Then,  determine  the  maximum  C  as  a  function  of 
max  A(C)  and  max  B(C). 


APPENDIX  A 


Aggregated  Arsenal  Exchange  Model 

This  appendix  details  the  aggregated  arsenal  exchange 
model  employed  in  this  study.  The  model  was  developed 
with  the  expert  assistance  of  Captain  Robert  Bunnell.  Its 
purpose  is  to  generate  three  different  response  surfaces 
for  the  methodology  example.  As  such,  its  complexity  is 
minimal  in  comparison  to  other  arsenal  exchange  models. 

Chapter  Three's  Tables  I  and  II  are  repeated  below 
as  a  description  of  the  target  class  and  weapon  class 
characteristics . 


TABLE  I 

Weapon  System  Characteristics 


W1 

W2 

W3 

W4 

W_ 

0 

CEP 

.200 

.150 

.  240 

.600 

.200 

Warheads/Weapon 

1 

3 

10 

4 

6 

Yield 

2.0 

.17 

.05 

.35 

.20 

Reliability 

.85 

.85 

.85 

.85 

.  85 

Daily  Alert 

.85 

.85 

.85 

.33 

.33 

Max  Warheads/Type 

450 

750 

1040 

424 

360 

Type 

Number 

Vulnerability 

Number 

Diameter 

CV1 

140 

24Q0 

.56 

CV2 

215 

13P1 

.49 

CV3 

430 

3 1P6 

.00 

CV4 

520 

23Q0 

.31 

CF1 

450 

35P7 

.00 

CF2 

1000 

52P8 

.00 

CF3 

200 

39P0 

.00 

CF4 

20 

22P1 

.  36 

CF5 

150 

11P0 

.00 

CF6 

100 

10Q1 

.79 

As  mentioned  in  Chapter  Three,  the  vulnerability  number 
characterizes  a  target's  vulnerability  to  its  most  suscept¬ 
ible  threat  (i.e.,  overpressure,  dynamic  pressure,  etc.). 

The  first  two  digits  are  its  level  of  susceptibility  to  the 
threat.  The  alpha  character  is  the  most  threatening  type  of 
effect  for  the  target  in  question.  The  last  number  is  the 
target's  sensitivity  to  the  duration  of  the  particular 
weapon's  effect.  Thus,  with  five  alpha  numeric  characters, 


the  target's  worst  case  vulnerability  is  defined. 


Likewise,  the  weapon's  ability  to  create  the  effects 
used  in  the  target  VNTK  number  is  calculated  from  its 
yield  in  megatons,  circular  error  probable  (CEP),  and  the 
probability  of  reaching  the  target  (reliability).  These 
effects  are  determined  in  a  cumulative,  lognormal  damage 
function  (1:3). 

The  target  and  weapon  characteristics  are  combined 
in  accordance  with  the  physical  effects  equations  of  the 
Defense  Intelligence  Agency's  Physical  Vulnerability 
Branch  (1:34)  to  determine  a  single-shot  probability  of 
survival  (SSP^)  for  each  target  when  attacked  by  each 

o 

weapon.  Thus,  for  this  example  of  five  weapon  types,  five 
SSPg  are  calculated  for  each  target.  The  transformation 
from  the  VNTK  number  and  weapon  characteristic  to  SSPs  is 
accomplished  by  a  computer  program  used  in  the  Bunnell  and 
Takacs  study  (5). 

The  SS?s  are  converted  into  single-shot  probabilities 
of  kill  (SSP^.)  for  each  target  and  weapon  combination, 

Eq  (30). 


SSP,  =  1  -  SSP 
k  s 

Likewise,  the  single-shot  probability  of  kill  is  calculated 
for  two  weapons  of  the  same  type  against  each  target  type, 
Eq  (31). 

SSP  "  =  1  -  (SSP  )" 

rv  S 


3  b 


Thus,  each  target  now  has  five  pairs  of  SSP^ 1 s  associated 
with  it  -  one  for  single  strike  and  one  for  two  strikes 
by  the  same  weapon. 

The  arsenal  exchange  is  modeled  using  linear  pro¬ 
gramming.  The  model  uses  the  SSP  ' s  and  SSP.  "s  as  the  co- 

rv  K 

efficient  for  the  objective  function.  The  number  of 
targets  struck  is  the  variable  for  the  tai-eau.  The  L? 
algorithm  maximizes  the  sum  of  damage  possible  as  deter¬ 
mined  by  the  SSP^.  coefficients  times  the  number  of  targets 
struck.  Eq  (32)  depicts  this  example  objective  function. 

5  10 

£  £  SSP,  •  T.  +  SSP,  "  •  T.  " 

W=1  K=1  w  w 

where 

k  =  target  struck  class 
w  =  weapon  class 

In  this  tableau,  a  separate  variable  is  used  for  a  target 
struck  twice  by  a  particular  weapon  class. 

The  first  ten  rows  of  the  tableau  define  the  sum  of 
each  class  of  targets  available,  Eq  (33). 

row(i)  =  *  of  targets  in  class  (i) 

=  f=l  T i  (  2  j  —  1  )  +  Ti  (  2  j  )  " 

where 

( 2  j  —  1 )  =  odd  columns  up  to  99 
(2j)  =  even  columns  up  to  100 

=  1  through  10 
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Rows  11  through  15  define  the  sum  of  each  class  of 
weapons  available,  Eq  (34). 

row  (i)  =  #  weapons  in  class  (i) 

50 

=  Z  Ti ( 2 j-1 ) +  2  *Ti ( 2 j ) " 
j=l 

where 

(2j-l)  =  odd  columns  up  to  99 

(2j)  =  even  columns  up  to  100 

i  =11  through  15 

The  difference  between  the  CF,  CV,  and  AF  models  is 
in  the  objective  function.  The  objective  variables,  TQj' 
which  are  designated  CF  or  CV  are  given  a  coefficient  SSP^ 
of  zero  if  not  in  the  exchange  under  consideration.  Thus, 
to  run  a  CF  exchange,  all  CV  targets  are  given  an  SSP^.  of 
zero. 

Since  the  model  was  to  have  multiple  runs  for  the 
weapons  available  in  rows  11  through  15,  a  program  was 
written  which  compiled  a  tableau  for  each  MOE  (CF,  CV,  AF ) 
This  program,  CREATE 2 ,  is  an  interactive  program  which 
requests  the  user  to  enter  the  right-hand-side  values 
(RHS )  for  the  number  of  weapons  available  in  each  class 
(rows  11  through  15).  The  program  then  requests  if  any 
other  tableau  values  should  be  changed  from  the  initial 
values  already  constructed  into  the  tableau.  Upon  com¬ 
pletion  of  these  two  data  entry  steps,  the  program  sends 
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the  tableau  for  each  option  to  local  files  AF,  CF,  and  CV, 
respectively.  The  tableau  format  is  determined  by  the 
linear  programming  package  "Multi-Purpose  Optimization 
System  (MPOS)."  The  source  code  for  this  Fortran  Five 
program  follows. 


PROGRAM  CREAIE2 

DIMENSION  AI.J<0:i5,0:t00)»CFlJ<0:itf<':t00).<:  ■  >0  ! 


THIS  PROGRAM  CREATES  THE  DArA  HA.'RIV  r OR  ,v*  A  >:;•  ' a r”  A-  ; 
EXCHANGE  MODEL  Tt;  BE  SOI  '.'ED  ' HR  ,%  .  T  0!  ha-  '?*i 

TARGET  CATEGORIES  A-!.'!  F r«  F  ur Ai* ON  *  1  7'.  ’  < 

A  SINGLE  TYr  E  WEAPON  AGAINST  A  EARTH  '-.AF  -ri  r 
OBJECTIVE  (ROW  0)  MAXIMIZES  *Hi:  Mr*A%  .  rn 

TOTAL  TARGET  BASE  .  TM  MODEL  ALI.'*WS  A  -  s' S'.  Wl  :  )N  -  A 
PAIR  OF  SAID  WEAPONS  TO  5".  Af.'lCA:  !Ti  71  A  [M  ..  '  ’ -.S  V 

NOTE!  THE  RIGHT-HAND  SIDE  VAUJZS'RO  AR  .  El  .  >  *? 

THE  OBJECTIVE  COEEF.  ARE  IN  ROW  ZEn's 


INITIALIZE  MATRIX  AIJ  TO  ALL  ZEROS . 

no  200.  j=o,ioo 
no  ioo,  1=0,15 

A!  J( 1 , J) =  0.0 
CONTINUE 
CONTINUE 


BUILDING  TARGET  CONSTRAINT  GET (ROWS  1  TO  10 

DO  400,  11=0,00,20 
DO  300,  1=1,10 
AIJO,<2*I+TI>)=  1.0 
ATj(If(2?l-l*-n))=  1,0 
CONTINUE 
CONTINUE 


nnnn  onooon  n  n  onoo 


BUILD  WEAPON  CONSTRAINT  SET (ROWS  11  TO  15) 


JA=t 

JB=19 

DO  800,  1=11,15 
DO  600,  J=JA, JB,2 
AIJ< I , J)=  1.0 
AIJ< I , ( J+l) )=2.0 
600  CONTINUE 

JA=JA+?0 

)B=JB+20 

800  CONTINUE 


SET  RHS  (C0I..0  ROWS  0  THRU  10)  FOR  OBJECTIVE  AND  NUMBER  OF 
TARGETS  IN  EACH  CATEGORY. 

AIJ(0,0>=  0,0 
AIJ<1,0>=  140.0 
A1J(2,0)=  215.0 
AIJ(3,0>=  450.0 
AIJ(4,0)=  1000.0 
AIJ(5,0)=  200.0 
AIJ(6,0)=  20.0 
AIJ(7,0)=  150.0 
AIJ<3,0)=  100.0 
AIJ(9,C)=  430.0 
AlJ( 10,0)=  520.0 


SET  PK'S  IN  OBJECTIVE  FUNCTI0N<R0W  0,  COLS.  1  TO  100) 

AIJ(0,1)=  0.34797 
AIJ(0,2)=  0.97689 
AIJ(0,3)=  0.85000 
AIJ(0,4)=  0.97750 
AI J( 0,5)=  0,83249 
AIJ(0,6)--  0,97194 
AIJ(0,7)=  0,47032 
AIJ(0,8)=  0.71944 
AIJ(0,9)=  0.67217 
AIJ(0,10)=  0.89253 
A I J ( 0 , 1 1 ) =  0,34999 
AIJ(0,12)=  0,97750 
AIJ(0,13)=  0.85000 
AIJ(0,14)=  0.97750 
AIJ(0,15)=  0.85000 


90 


AIJ<0,16)= 
AIJ<0,17)= 
AIJ<0,18>= 
AIJ<0,19)= 
AIJ(0,20)= 
AIJ<0,21)= 
AIJ<0r22)= 
AIJ<0,23)= 
AIJ<0f24)= 
AIJ<  0» 25>  = 
AIJ<0,26>= 
AIJ<0f27)= 
AIJ<0,28)= 
AIJ<0,2?)= 
At J(Of 30)= 
AIJ(Of 31 )= 
AIJ<0,32>= 
AIJ(Of 33)= 
AIJ<0,34)= 
AIJ<0,35)= 
AIJ<0,36)= 
AIJ( Of 37)= 
AI J(0f 38)= 
AIJ(0,39)= 
AIJ(0,40)= 
AIJ(0,41 )= 
AX J(0,42)= 
AIJ<0,43)= 
AIJ(0,44)= 
AIJ(0,45)= 
AIJ<0,46)= 
AIJ(0,47)= 
AI J(0,48)= 
AIJ(0f49)= 
AI JtOf S0)= 
AIJ<0,51)= 
AIJ(0,52)= 
ACJ<0,53)= 
AI J(Of 54)= 
AIJ<0,55)= 
AIJ(Of 56)= 
AIJ(0,57)= 
AIJ(0,58)= 
AIJ(0f59)= 
AlJtOf 60)= 
AI  J(Of  6).  )  = 
AIJ(0,62)= 
ALJ<0,63)  = 
AI J<0,64)= 
A) J(0f65)= 
AI J(0f66)= 


0.97750 

0.84660 

0.97647 

0.84979 

0,97744 

0.76716 

0.94579 

0.35000 

0,97750 

0.61073 

0.84847 

0.15429 

0.23473 

0.37563 

0,61016 

0.83104 

0.97145 

0.85000 

0.97750 

0.85000 

0.97750 

0.71782 

0.92038 

0.81872 

0,96714 

0.44306 

0.68982 

0.84164 

0.97492 

0,14544 

0.26973 

0.02663 

0.05255 

0,08063 

0,15475 

0.51804 

0,76771 

0,84982 

0,97744 

0.84097 

0.97471 

0.22121 

0.39349 

0.52353 

0,77298 

0,36172 

0.59260 

0.32551 

0.96955 

0.11517 

0.21708 


o  o  n  o  o  n 


A1J<0,67)= 

0.02402 

A I J ( 0 , 60  >  = 

0.04747 

AIJ<0,69)= 

0 ,04995 

AI J(0,70)= 

0,09741 

AIJ<0,71 )= 

0.40469 

AI J(0,72)= 

0.64561 

AIJ<0,73)= 

0.84494 

AI J(0,74)= 

0.97596 

AIJ<  0,75)  = 

0.83678 

AIJ(0,76>= 

0,97336 

AIJ<0,77)= 

0.17125 

AI J<0,70)= 

0.31317 

AIJ(0,79)= 

0.40083 

AIJ(0,80)= 

0.64099 

AI J(0,81 )= 

0.72773 

AI J(0,82)= 

0.92587 

AIJ(0,83)= 

0.85000 

AT J(0,34)= 

0.97750 

AIJ(0,85)= 

0.48420 

AIJ(0,86)= 

0,73395 

AIJ<0,87)= 

0.10270 

AI J(0,88)= 

0.19496 

AIJ(0,89)= 

0.26052 

AI J(0,90)= 

0.45317 

AIJ<0,91)= 

0.80019 

AIJ(0,92)= 

0.96008 

AIJ(0,93)= 

0.85000 

AIJ(0,94)= 

0.97750 

AIJ(0,95)= 

0,85000 

AI J(0,9A)= 

0,97750 

AIJ(0,97)= 

0.60865 

AI J(0,98)= 

0,84685 

AIJ(0,99)= 

0.78184 

AIJ(0,100)= 

0.95241 

RHS  VALUES  REQUESTED  FOR  EACH  WEAPON  CATEGORY , 

ANY  COEFFICIENTS  MAY  BE  CHANGED  IN  THIS  SECTION.  THE  SESSION 
IS  TERMINATED  WHEN  THREE  ZEROS  ARE  ENTERED. 


900  PRINT*, 'INPUT  ROW  NUMBER, COLUMN  NUMBER, COEFFICIENT: < I ,J,CCEFF) ' 
PRINT*, 'FOR  EACH  WEAPON  RHS  AND  ANY  OTHER  COEFFICIENT  TO  BE' 
PRINT*, 'CHANGED.  ENTER  ONLY  ONE  SET  AT  A  TIME,  TO  TERMINATE' 
PRINT*, 'ENTER  THREE  ZEROS  (000)' 

READ*,  I,J,C 

I F ( I ♦ EG . 0 » AND ♦ J ♦ EQ . 0 . AND « C . EQ . 0 )  GO  TO  902 

AIJ(I,J)=C 

GOTO  900 
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902  PRINT*, 'ENTER  THE  RHS  WEAPON  ALLOCATION  VALUES' 

DO  46,1*11,15 
READ* jC 
AI. ”1,0=0 
46  CONTINUE 
C 

950  PRINT*, 'THE  FILE  CONTAINING  YOUR  ARSENAL  MODEL  MATRIX' 
PRINT*, 'IS  CALLED-MN  FOR  ALL  FORCES  AND  CV/CF  FOR  OPTIONS.' 

C 

c 

c 

C  ARSENAL  MATRIX  WRITF.N  TO  FILE  MTX .  ARSENAL  MATRIX  ALSO 

C  COPIED  FOR  CV  AND  CF  OPTIONS. 

C 

0PEN(12,FILE='MN' ) 

REWIND  12 
DO  954,  1=0,15 
DO  952,  J=0,1.00 

IF(AIJ(I,J).NE.O.O)  URITE<rj,951>  I,J,AIJO.  ,J) 

951  FORMAT (2X,I3,2X,I3,2X,FtO»5) 

CFI J( I , J)=AI J( I , J) 

CVIJCI  ,«J)=AIJ<  I ,  J ) 

C 

952  CONTINUE 
954  CONTINUE 

CL0SEU2) 

C 

C  COUNTER-FORCE  MATRIX*.  SET  PK=0.0  FOR  CV  TARGETS  IN  OBJECTIVE 

C 

OPEN' 11 ,FILE='CF' ) 

REWIND  11 
DO  13,  11=0, SO, 20 
DO  11,  1=1,4 
CFI J (0,1+11) =  0.0 

11  CONTINUE 

DO  12,  1=17,20 
CFIJIO, I+II)=  0.0 

12  CONTINUE 

13  CONTINUE 
C 

C 

C 

C  COUNTER-VALUE  MATRIX?  SET  PK=0.0  FOR  CF  TARGETS  IN  OBJECTIVE 

C 

OPEN! 13,FILE='CV' ) 

REWIND  13 
DO  16,11=0,80,20 
DO  15,1=5,16 
CVIJIO.I+Ii )=  0.0 
1.5  CONTINUE 

16  CONTINUE 

r 
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DO  21#  1=0. 15 

no  2o,j=o,ioo 

IF(CFIJ(I,J).NE»0.0)  UR I T  E  < 1 1 , 2  ?  >  I # J , CF  I  J ( I , J ) 
TF(CVIJ(I,J).NE.0,0>  WRITE<13,29>  I , J,CVIJ< I , J) 
CONTINUE 
CONTINUE 

F0RMAT(2X,I3,2X,I3,2X#F10.5) 

CL0SE<11) 

CLOSE (13) 


The  linear  programming  package  chosen  to  optimize 
the  tableaux  is  the  Multi-Purpose  Optimization  system. 

Thanks  to  the  work  of  fellow  student#  James  Cooke#  several 
flaws  in  the  MPOS  options  were  discovered  prior  to  this 
study.  Thus#  the  MIMIT  option  was  chosen  (7:38).  This 
option  solves  a  "sequence  of  primal  and/or  dual  iterations#” 
(7:38.1).  Samples  of  the  MPOS  control  cards  for  each  MOE 

and  a  sample  MPOS  output  are  given  below. 

CF  Control  cards: 

♦THIS  IS  THE  CONTROL  CARD  DECK  FOR  MPOS  AEH 
TITLE 

ARSENAL  EXHANGF.  PROBLEM 

MINIT 

VARIABLES 

XI  TO  X100 

PACKED 

MAXIMIZE 

CONSTRAINTS  15 

E+m++m++m 

FORMAT 

(2X,I3,2X#I3,2X,F10.5> 

REWIND 
READ  CF 
OPTIMIZE 


rT*vn 


v*  ‘J  t* 


CV  Control  cards: 


♦THIS  IS  THE  CONTROL  CARD  DECK  FOR  MPCI3  AFM 

TITLE 

ARSENAL  EXHANGE  PROBLEM 
MINIT 
VARIABLES 
XI  TO  X100 
PACKED 
MAXIMIZE 
CONSTRAINTS  15 
+++++++++++++++ 

FORMAT 

<2X,I3,2X,13,2X,F10.5> 

REWIND 
READ  CV 
OPTIMIZE 


AF  Control  cards: 


♦THIS  IS  THE  CONTROL  CARD  DECK  FOR  MF'OS  AEM 
TITLE  MCN 

ARSENAL  EXHANGE  PROBLEM 

MINIT 

VARIABLES 

XI  TO  XlOO 

PACKED 

MAXIMIZE 

CONSTRAINTS  15 

+++++++++++++++ 

FORMAT 

<2X,I3,2X,I3f2X,F10.5> 

REWIND 
READ  MN 
OPTIMIZE 
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I ist,f=ocf 


1  MPOS  VERSION  4.0  NORTHWESTERN  UNIVERSITY 

************************************* 

*  % 

*  MPOS  * 

*  * 

*  VERSION  4.0  t 

*  * 

*  MULTI-PURPOSE  OPTIMIZATION  SYSTEM  * 

*  * 
************************************* 


*****  PROBLEM  NUMBER  1  ***** 

*THIS  IS  THE  CONTROL  CARD  DECK  FOR  MPOS  AEM 
TITLE 

ARSENAL  EXHANGE  PROBLEM 

MINIT 

VARIABLES 

XI  TO  X100 

PACKED 

MAXIMIZE 

CONSTRAINTS  15 

+++++-H-++++++++ 

FORMAT 

<2XrI3f2Xf I3f2XfFl0.5) 

REWIND 
READ  CF 
OPTIMIZE 

1  MPOS  VERSION  4.0  NORTHWESTERN  UNIVERSITY 

********************* 

*  PROBLEM  NUMBER  1  * 

********************* 

USING  MINIT 

ARSENAL  EXHANGE  PROBLEM 
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INPUT  TRANSLATION  TIME  =  .1810  SECONDS 


DATA  STORAGE  REQUIREMENTS  ARE  004264  OCTAL 
FIELD  LENGTH  FDR  THIS  PROBLEM  050000  OCTAL 
1  MPOS  VERSION  4,0  NOPTHWESTE!  N  UN T VERS i TY 

ttMttt-Mtttrx******* 

*  PROBLEM  NUMBER  3.  * 

^  i  ^  ^  'it  ^ +  'A'  ♦  ^  4'  'k  't- 

*  •*  ^  ^  /p  ‘P  *  p  ^  ^  ^  ^  ip  ^  ^ 

USING  MINIT 

ARSENAL  EXHANGE  PROBLEM 


ENTERING  PHASE  II 


ITERATION 

1 

IN 

VAR- 

5 

OUT 

VAR- 

103 

-ZMIN* 

374 

.621 

ITERATION 

o 

IN 

VAR- 

34 

OUT 

VAR- 

3.07 

-ZMIN* 

52 1 

.2-16 

ITERATION 

3 

IN 

VAR- 

30 

OUT 

VAR- 

105 

-ZMIN* 

643 

♦  »:  /  • : 

ITERATION 

4 

IN 

VAR- 

96 

DOT 

VAR- 

103 

-ZMIN® 

741 

.028 

iteration 

b 

IN 

VAR- 

47 

OUT 

t'AR- 

113 

-ZMIN* 

764 

.  462 

ITERATION 

6 

IN 

VAR- 

32 

OUT 

VAR- 

106 

-ZMIN* 

733 

,393. 

ITERATION 

7 

IN 

VAR- 

38 

OUT 

VAR- 

115 

-ZMIN* 

799 

.  483 

ITERATION 

3 

IN 

VAR- 

56 

OUT 

VAR- 

96 

-ZMIN* 

33.3 

,37? 

ITERATION 

9 

IN 

VAR- 

68 

OUT 

VAR- 

104 

-ZMIN* 

320 

,024 

ITERATION 

10 

IN 

VAR- 

23 

OUT 

VAR- 

13.2 

-ZMIN* 

B21 

,211 

ITERATION 

1.1 

IN 

VAR- 

74 

OUT 

VAR- 

68 

-ZMIN* 

353 

,040 

ITERATION 

12 

IN 

VAR- 

52 

OUT 

VAR- 

34 

-ZMIN* 

957 

.705 

ITERATION 

13 

IN 

VAR- 

92 

OUT 

VAR- 

32 

-ZMIN* 

35S 

,093 

ITERATION 

14 

IN 

VAR- 

48 

OUT 

VAR- 

52 

-ZMIN* 

853 

.847 

ITERATION 

15 

IN 

VAR- 

27 

OUT 

VAR- 

28 

-ZMIN* 

362 

,  388 

ITERATION 

16 

IN 

VAR- 

87 

OUT 

VAR- 

47 

-ZMIN* 

384 

,347 

ITERATION 

17 

IN 

VAR- 

76 

OUT 

VAR- 

114 

-ZMIN* 

364 

,592 

ITERATION 

18 

IN 

VAR- 

54 

OUT 

VAR- 

56 

-ZMIN* 

864 

,604 

ITERATION 

19 

IN 

VAR- 

6 

OUT 

VAR- 

111 

-ZMIN* 

364 

,604 

ITERATION 

20 

IN 

VAR- 

7 

OUT 

VAR- 

6 

-ZMIN* 

364 

.604 

ITERATION 

21 

IN 

VAR- 

25 

OUT 

VAR- 

27 

-ZMIN* 

397 

con 
♦  <J  7  w 

ITERATION 

22 

IN 

VAR- 

90 

OUT 

VAR- 

5 

-ZMIN* 

904 

♦  335 

iteration 

23 

IN 

VAR- 

36 

OUT 

VAR- 

88 

-ZMIN* 

906 

,  086 

ITERATION 

24 

IN 

VAR- 

47 

OUT 

VAR- 

90 

-ZMIN* 

91  4 

.463 

ITERATION 

25 

IN 

VAR- 

26 

OUT 

VAR- 

37 

-ZMIN* 

923 

n  n  n 
♦  OOi', 

ITERATION 

26 

IN 

VAR- 

9 

OUT 

VAFi- 

25 

-ZMIN* 

93.1 

,6se 

ITERATION 

27 

IN 

VAR- 

52 

OUT 

VAR- 

92 

-ZMIN* 

932 

♦  oif  O 

1  MPOS  VERSION  4.0  NORTHWESTERN  UNIVERSITY 

*:C  ******************:!: 

*  PROBLEM  NUMBER  3.  * 

USING  MINIT 

ARSENAL  EXHANGE  PROBLEM 


.r-C  «*«.«*« 


PI 
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SUMMARY  OF  RESULTS 


:  I ABLE 

VARIABLE 

BASIS 

ACTIVITY 

OPPORTUNITY 

NO* 

NAME 

NON --BASIC 

LEVEL 

U  u  >  1 

, 

XI 

MB 

— 

.  4696100 

n 

VO 

NB 

— 

.  9392200 

3 

X3 

NB 

— 

.4696100 

4 

X4 

NB 

— 

.9392200 

5 

X5 

NB 

— 

.0779900 

6 

X6 

NB 

— 

.4081500 

7 

X7 

B 

355.0000000 

— 

8 

X8 

NB 

— 

.2204900 

9 

X9 

B 

95.0C00000 

— 

10 

X10 

NB 

— 

.2492500 

11 

Xll 

NB 

— 

.3354900 

12 

X12 

NB 

— 

.6775900 

13 

X13 

NB 

— 

,5452100 

14 

X14 

NB 

— 

.8873200 

15 

X15 

NB 

— 

.5426100 

16 

XI 6 

NB 

— 

.3347200 

17 

X17 

NB 

— 

.4696100 

IS 

X18 

NB 

— 

.9392200 

19 

X19 

NB 

— 

,4696100 

20 

X20 

NB 

— 

.9392200 

21 

X21 

NB 

— 

.2038000 

22 

X22 

NB 

— 

.4076000 

23 

X23 

NB 

— 

.2038000 

24 

X24 

NB 

•  — 

,4076000 

25 

X25 

NB 

— 

.0339400 

26 

X26 

B 

270.0000000 

— 

27 

X27 

NB 

— 

,0502200 

23 

X23 

NB 

— 

.1235300 

29 

X29 

NB 

— 

.0307300 

30 

X30 

B 

105.0000000 

— 

31 

X31 

NB 

— 

.0886300 

32 

X32 

NB 

— 

,1520200 

33 

X33 

NB 

— 

. 2794000 

34 

X34 

NB 

— 

,3557000 

35 

X35 

NB 

— 

.2763000 

36 

X36 

NB 

— 

,3531000 

37 

X37 

NB 

— 

. 2033000 

33 

X38 

NB 

— 

.4076000 

39 

X39 

NB 

— 

.2038000 

40 

X40 

NB 

— 

, 4076000 

41 

X41 

NB 

— 

,0259200 

42 

X42 

NB 

— 

,0513400 

43 

X43 

NB 

-- 

.0259200 

44 

?'44 

NB 

— 

,0513400 

45 

X  '  5 

NB 

— 

.3213300 

46 

X46 

NB 

— 

-  2229300 
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SUMMARY  0 

if  RESULTS 

MRIABi 

-E  VARIABLE 

BASIC 

ACT  I '7TTY 

Cr  PflRTiJNITY 

MfJ* 

NAME 

NON- BASIC 

'  i 

,--fV  ;T 

47 

X47 

iT' 

63C *  0000000 

..... 

48 

X48 

8 

7  *  0000000 

_ 

49 

X4? 

NB 

— 

.  1473500 

50 

X50 

NB 

— 

.0996500 

51 

X51 

NB 

— 

.2237500 

52 

X52 

B 

20.0000000 

— 

53 

X53 

NB 

— 

.1017000 

54 

X54 

B 

94.0000000 

_ 

55 

X55 

NB 

— 

.1079500 

56 

X56 

NB 

— 

.0001300 

57 

X57 

NB 

— 

. 0259200 

58 

X58 

NB 

-- 

.0518400 

59 

X59 

NR 

— 

.0259200 

60 

X60 

NB 

— 

,0513400 

61 

X61 

NB 

.... 

.0251300 

62 

X62 

NB 

.... 

.0503600 

63 

X63 

NB 

.0251300 

64 

X64 

NB 

.0503600 

65 

X63 

NB 

— 

.3503300 

66 

X66 

NB 

— 

.2741500 

67 

X67 

NB 

— 

►0013700 

68 

X68 

NB 

— 

.0036000 

69 

X69 

NB 

— 

.1777900 

70 

X70 

NB 

— 

.1555100 

71 

X71 

NB 

— 

.3363600 

72 

X72 

NB 

— 

,1206200 

73 

X73 

NB 

— 

.1053400 

74 

X74 

b 

56.0000000 

— 

75 

X75 

NB 

— 

.1114000 

76 

X76 

B 

100.0000000 

_ 

77 

X77 

NB 

— 

.0251300 

78 

X78 

NB 

, 0503600 

79 

X79 

NB 

,0251300 

80 

X80 

NB 

— 

.0503600 

31 

X31 

NB 

— 

•  1465400 

82 

X82 

NB 

— 

,2930800 

83 

X83 

NB 

— 

,1465400 

84 

X34 

NB 

— 

,2930800 

85 

X85 

NB 

.1032100 

86 

X86 

B 

130,0000000 

— 

87 

X87 

NB 

— 

.0445500 

88 

X38 

NB 

— 

,0983300 

39 

X69 

NB 

— 

. 0385300 

90 

X90 

NB 

— 

.0424700 

9: 

X91 

NB 

— 

.0622200 

92 

v  n  n 

A  7  * 

NB 

— 

-  0483700 

HPQS 

VERSION  4  * 0 

NORTH 

UHSTF.LM  UNIVERSE 

TV 

100 


JIMWLm*'.' 


\  'A  'A  .1 


WW-T «  *■-  > 


:K3:*S:i*St:t^*ric,ir*:^::!c:)t**S'4! 

$  PROBLEM  NUMBER  l  * 

*  :'<  *  %  %  *****  %  *  %  %  ******* 

USING  MINIT 

ARSENAL  EXHANGE  PROBLEM 


SUMMARY  OF  RESULTS 


VARIABLE 

VARIABLE  BASIC 

ACTIVITY 

OPPORTUNITY 

ROW 

NO. 

NAME 

NON-BASIC 

LEVEL 

COST 

NO. 

93 

X93 

N3 

— 

.2221400 

94 

X94 

NB 

— 

.2411800 

95 

X95 

NB 

— 

.2195400 

96 

X96 

NB 

— 

.2383800 

97 

X97 

NB 

— 

.1465400 

98 

X98 

NB 

— 

.2930800 

99 

X99 

NB 

— 

. 1465400 

100 

X100 

NB 

— 

,2930800 

101 

—SLACK 

B 

140.0000000 

— 

( 

1) 

102 

—SLACK 

B 

215.0000000 

— 

( 

2) 

103 

—SLACK 

NB 

— 

.4408700 

( 

3) 

104 

—SLACK 

NB 

— 

.0007100 

( 

4) 

105 

—SLACK 

NB 

— 

.2025600 

( 

5) 

106 

—SLACK 

NB 

— 

,7153700 

( 

6) 

107 

—SLACK 

NB 

— 

.9256000 

( 

7) 

108 

—SLACK 

NB 

— 

.9230000 

( 

8) 

1.09 

—SLACK 

B 

430.0000000 

. — 

( 

9) 

110 

—SLACK 

B 

520.0000000 

— 

( 

10) 

111 

—SLACK 

NB 

— 

.4696100 

( 

11) 

112 

—SLACK 

MB 

— 

,2038000 

( 

12) 

113 

—SLACK 

NB 

— 

.0259200 

( 

13) 

114 

—SLACK 

NB 

— 

.0251800 

( 

14) 

115 

—SLACK 

NB 

— 

.1465400 

\ 

15) 

MAXIMUM 

VALUE  OF 

THE  OBJECTIVE 

FUNCTION  = 

932.665560 

CALCULATION  TIME 

WAS  .0720 

SECONDS  FOR 

27  ITERATIONS. 

IT0TAI  TIME 

FOR  THIS 

PROBLEM  WAS 

.359  SECONDS 

APPENDIX  B 


Research  Design 

The  research  design  chosen  for  this  study  is  a  three- 
level,  3  ,  fractional  factorial  design  (2).  It  is  a  four- 
factor  design  which  requires  twenty-seven  samples  for  the 
model  with  variance.  For  thi$  study,  the  LP  model  has 
no  variance;  therefore,  the  repetitions  of  the  center  point 
sample  are  not  required.  Thus,  only  twenty-five  samples 
are  required. 

The  design  is  shown  with  its  respective  sample  values 
in  the  following  examples.  The  first  example  is  the  un¬ 
coded  design  and  each  of  the  respective  sampled  MOE ' s 
(CF,  CV,  AF).  This  design  was  used  for  the  revised  lower 
limits  of  300,  450,  720,  and  200  for  each  Wi  respectively. 
The  remaining  examples  are  in  coded  format  (+1,  0,  -I). 

This  format  assigns  one  for  the  upper  limit  of  the  variable, 
zero  for  the  mid-range  of  the  limit,  and  negative  one  for 
the  lower  limit.  These  coded  designs  were  used  in  the 
initial  lower  VL  limits  of  zero. 


Sampled  MOE's  and  Research  Design  for  Revised  Limits 


CF 

CV 

AF 

W1 

W2 

W3 

W4 

932.67 

1190.02 

1736.46 

*50 

750 

sso 

863.40 

1159.59 

1618.73 

450 

450 

880 

312 

862.12 

1170.66 

1652.36 

300 

750 

330 

312 

797.35 

1139.65 

1504.04 

300 

450 

880 

312 

873.84 

1165.56 

1677.36 

375 

600 

1 040 

424 

868.13 

1165.56 

1645.10 

375 

600 

1040 

200 

865.41 

1164,70 

1626.09 

375 

600 

720 

424 

859.61 

1164.31 

1535.47 

375 

600 

720 

200 

904.92 

1174,87 

1701 .47 

450 

600 

880 

424 

899.15 

1174.81 

1688 . 05 

450 

600 

880 

200 

834 , 40 

1155,39 

1602.27 

300 

600 

880 

424 

828.60 

1155.32 

1567.22 

300 

600 

880 

200 

901 . 59 

1130,83 

1717.52 

375 

750 

1040 

312 

393.14 

1179.31 

1673.61 

375 

750 

720 

31? 

837.39 

1150.23 

1594.91 

375 

450 

1040 

312 

328.37 

1143.42 

1526.33 

375 

450 

720 

312 

906.24 

1175.30 

1708.49 

450 

600 

1040 

312 

397.79 

1174.23 

1657.15 

450 

600 

720 

312 

335.76 

1155.3? 

1612.06 

300 

600 

1040 

312 

827.24 

1153.96 

1554.73 

300 

600 

720 

312 

900.26 

1130.40 

1710.36 

375 

750 

880 

424 

394.50 

1180.34 

1631.32 

375 

750 

380 

200 

836.03 

1149.35 

1585.12 

375 

450 

880 

424 

330.23 

1149.79 

1540.05 

375 

450 

880 

200 

366.77 

1165.13 

1635.83 

375 

600 

830 

312 

The  next  two  examples  are  for  coded  CF  and  CV  research 
designs  and  their  sampled  values.  These  designs  were  used 
in  the  initial  lower  limit  of  zero. 


APPENDIX  C 


Regression 

This  section  includes  a  complete  stepwise  regression 

for  CF,  and  the  final  step  and  residual  scatter  diagrams 

for  CV  and  AF  regressions.  It  concludes  with  Tables  VI, 

% 

VII,  and  VIII  which  show  the  regression  error  for  each 
surface. 

The  Statistical  Package  for  the  Social  Sciences 
(SPSS)  was  employed  with  a  stepwise  regression  option  for 
the  determination  of  the  regression  polynomials. 

Example  CF  Stepwise  Regression: 


1 

T 

84/03/31.  15.54.17.  PAGE 

ASD  COMPUTER  CENTER 
WR1GHT-PATTERS0N  AFB,0HI0 

SPSS--  STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES 
VERSION  8.3  <N0S)  —  MAY  04,  1782 


376500  CM  MAXIMUM  FIELD  LENGTH  REQUEST 


KTOiTO*ra*roiire^;  ■>  ■>  ■;  -r;  ’.^vw  ?■?■.»,?? 


RUN  NAME 

FILE  NAME 

VARIAELE  LIST 

N  OF  CASES 

INPUT  MEDIUM 

INPUT  FORMAT 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 


SECOND  ORDER  REGRESSION  FOR  4  INDEP.  VARIABLES 
REGRESSION  ANALYSIS  WITH  RGDAT 
CFtCV,AF,U1 Ty2,W3,U4 
25 

CARD 

FREEFIELD 

W11=U1**2 

W22=W2**2 

W33=W3*#2 

W44=W4**2 

W12=W1*W2 

W13=W1*W3 

W14=W1*W4 

W23-W2*M3 

W24=W2&W4 

W34=W3*W4 


CPU  TIME  REQUIRED. , 


.044  SECONDS 


REGRESSION  METHOD=STEPWISE/ 

VARIABl.£S=CF  TO  W4,W11 ,W22fU33,W44,W12,Wl3,W14,W93,W24,W34 
REGRESSIGN=CF  WITH  W.t  TO  W34/RESIDUALS/ 

REGRESSIONS  WITH  W1  TO  W34/RESIDUALS/ 

REGRESSIONS  WITH  W1  TO  W34/RESIIHJALS/ 

STATISTICS  ALL 


00056000  CM  NEEDED  FOR  REGRESSION 


1 SECOND  ORDER  REGRESSION  FOR  4  INDEP.  VAR  84/03/31.  15.56.17.  PAGE  2 
FILE  -  REGRESSI  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGDAT 

mumHHULT  IPLE  REGRESSION********* 


VARIABLE 

MEAN 

STANDARD  DEV 

CASES 

CF 

866.0364 

33.9360 

25 

CV 

1164.9860 

12.9165 

25 

AF 

1631.3460 

63,6308 

25 

W1 

375.0000 

53.0330 

25 

W2 

600.0000 

106,0660 

25 

W3 

880.0000 

113,1371 

■^5 

•u. 

W4 

312.0000 

79.1960 

25 

Wll 

143325.0000 

39373,0367 

25 

W22 

370800.0000 

127795.2464 

25 

W33 

786688.0000 

199543.6745 

25 

105 
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HN 


fiftd 


1444 

103365.1200 

49330,490 4 

25 

1412 

225000.0000 

51 143 >1314 

25 

(413 

330000,0000 

63261,3626 

25 

(414 

117000.0000 

34167,2845 

25 

1423 

523000.0000 

115327.4579 

25 

(424 

137200.0000 

53310.2049 

25 

W34 

274560.0000 

73463,7564 

25 

n? 

& 


CORRELATION  COEFFICIENTS. 

A  VALUE  OF  99.00000  IS  PRINTED 

IF  A  COEFFICIENT  CANNOT  BE  COMPUTED. 


CV 

.96255 

AF 

.94427 

.93601 

U1 

.73483 

.53064 

.55301 

142 

.66939 

.04161 

.74357 

0 

143 

.08836 

.03570 

.30704 

0 

0 

(44 

.06015 

.00292 

.19999 

0 

0 

0 

Wll 

.73340 

.53717 

.55226 

.99741 

0 

0 

(422 

.66466 

.83810 

,73667 

0 

.99596 

0 

(433 

.03858 

,03541 

.30665 

0 

0 

.99736 

(444 

.06047 

.00376 

,20068 

.00000 

,00000 

.00000 

W12 

.97779 

.98924 

,91788 

.62217 

.77771 

0 

(413 

,60135 

.42093 

,61314 

.73772 

0 

.67065 

1414 

.40811 

.26336 

.44135 

.48424 

0 

0 

1423 

.59119 

.69868 

,77587 

0 

,80584 

.53606 

U24 

.42890 

.48001 

.58204 

0 

,56753 

0 

(434 

.09315 

,01846 

,31451 

0 

0 

. 44937 

CF 

CV 

AF 

W1 

142 

W3 

Wll 

0 

(422 

0 

-.00182 

(433 

0 

-.00133 

-.00166 

W44 

.99173 

-.00260 

-.00325 

-.00237 

(412 

0 

.62056 

,77457 

0 

0 

(413 

0 

.73581 

0 

.66922 

0 

,43899 

(414 

.86916 

,48299 

0 

0 

,86197 

.30128 

(423 

0 

0 

,80258 

.58481 

0 

.62671 

1424 

.81491 

0 

,56524 

0 

.30317 

.  44133 

U34 

,38821 

0 

0 

.44891 

*  SS03t> 

0 

W4 


Wll 


W22 


W33 


W44 


!ji  ■ 
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U14 

.35724 

W23 

.39305 

0 

1424 

0 

.7022? 

.45733 

(434 

,30171 

.77200 

.26365 

♦  /  4>  V.1 

U13  U14  (423  U24 


1 SECOND  ORDER  REGRESSION  FOR  4  INDEP.  VAR  34/03/31.  15.56.17.  PAGE 


FILE  -  REGRESSI  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGDAT 

*********multiple  regression********* 


DEP.  VAR...  CF 


MEAN  RESPONSE  066,03640 

STD.  DEV. 

33,93593 

VARIABLE (S) 
1412 

ENTERED 

ON  STEP 

1 

MULTIPLE  R 

.9778 

ANOVA 

DF  SUM  SQUARES 

MEAN  SQ. 

F 

R  SQUARE 

.9561 

REGRESSION 

1,  26425.270 

26425.270 

500.493 

STD  DEV 

7.2662 

RESIDUAL. 

23.  1214.352 

52.793 

SIG.  0 

ADJ  R  SQUARE 

.9542 

COEFF  OF  VARIABILITY  ,3PCT 

VARIABLE 

B 

S.E.  B 

F  SIG. 

BETA  ELASTICITY 

1412 

.001 

.000 

500.498  0 

,97779 

.16856 

CONSTANT 

720.055 

6.685 

11601.442  0 

****** 

VARIABLE (S) 
1413 

*  *  *  * 

ENTERED 

*  *  *  *  * 

ON  STEP 

****** 

n 

Am 

*  *  *  * 

********** 

MULTIPLE  R 

.9927 

ANOVA 

DF  SUM 

SQUARES 

MEAN  SQ.  F 

R  SQUARE 

.9855 

REGRESSION 

!  2,  2 

7240.206 

13620.103  750,201 

STD  DEV 

4.2609 

RESIDUAL 

22. 

399,416 

18.155  SIG.  .000 

ADJ  R  SQUARE 

:  ,9842 

COEFF  OF  VARIABILITY 

,5PCT 

VARIABI  E 

B 

S.E,  B 

F 

SIG. 

BETA  ELASTICITY 

1412 

.001 

.000 

949.877 

0 

.88903  .1532? 

(413 

.000 

,000 

44,887 

.000 

.19327  .03951 

CONSTANT 

699.085 

5,016 

19421.403 

0 

1 SECOND  ORDER  REGRESSION  FCS  4  INDEP.  VAR  34/03/31.  15.56.1"’.  PACE  4 
FILE  -  REGRESSI  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGI'/J 


107 


rZSLSt 


W23 

-.000 

.000 

/  *  u  0  / 

.016 

-.00249 

-.000*4 

W33 

-.000 

.  000 

7,593 

.017 

-.00360 

-.00056 

W24 

-.000 

.000 

2.549 

,136 

-.00102 

-.00013 

W44 

-.000 

.000 

1.47? 

.247 

-.00031 

-.00007 

CONSTANT 

504,445 

.360 

.  196EF07 

0 

1 SECOND  ORDER  REGRESSION  FOR  4  INDEP.  VAR  84/03/31.  15.56.17.  PAGE  12 


FILE  -  REGRESSI  (CREATED  -  34/03/31)  ANALYSIS  WITH  RGDAT 

#*#*##**#MUITIP  LE  REGRESSION##***#*# 


DEP.  VAR...  CF 

VARIABLE'S)  ENTERED  ON  STEP  17 
Wli 


MULTIPLE  T. 

1.0000 

ANOVA 

DF  SUM  SQUARES 

MEAN  39, 

R  SQUARE 

1,0000 

REGRESSION 

13. 

27639. 

620 

2126.125 

,  125E+08 

STD  DEV 

,0131 

RESIDUAL 

11. 

♦ 

002 

,000 

SIG.  0 

ADJ  R  SQUARE 

1.0000 

COEFF  OF  VARIABILITY 

,OPCT 

VARIABLE 

B 

S.E.  B 

F 

SIG. 

BETA  ELASTICITY 

W13 

-.000 

.000 

7.186 

.021 

-.00272 

-.00056 

W14 

-.000 

.000 

1.320 

.275 

-.00090 

-.00012 

W2 

.296 

.000 

.357E+06 

0 

.92361 

,20473 

W1 

.472 

.001 

. 163E+06 

0 

.73693 

.20419 

W22 

-.000 

.000 

37746.079 

0 

-.25266 

-.02873 

U3 

,029 

.001 

2272,737 

0 

,09629 

>02935 

W4 

.028 

.001 

1977,203 

0 

.06502 

.01004 

W34 

-.000 

.000 

11.87? 

,005 

-.00290 

-.00040 

U23 

-.000 

.000 

7,186 

.021 

-.00249 

-.00044 

U33 

-.000 

.000 

3.523 

.  OR  7 

-.00335 

-.00052 

W24 

-.000 

.000 

2.347 

,154 

-.00102 

-.00013 

W44 

-.000 

.000 

.564 

.449 

-.00068 

-.00006 

Ull 

.000 

.000 

,046 

.834 

.00035 

.00005 

CONSTANT 

504.542 

.590 

.732E+06 

0 

F-LEVEL  OR  TOLERANCE-LEVEL  INSUFFICIENT  FOR  FURTHER  COMPUTATION 


COEFFICIENT!:  AND  CGNFIDFNCE  INTERVALS. 
VARIABLE  B  95  PCT  C.I. 


W13 

-.0000 

-.0000 

-.0000 

W14 

-.0000 

-.0000 

.0000 

142 

.2955 

.2944 

.2966 

141 

.4716 

.4690 

.4741 

1422 

-.0001 

-.0001 

-.0001 

W3 

.0289 

.0276 

.0302 

144 

.0279 

.0265 

.0292 

W34 

-.0000 

-.0000 

-.0000 

1423 

-.0000 

-.0000 

-.0000 

1433 

-.0000 

-.0000 

.0000 

W24 

-.0000 

-.0000 

.0000 

W44 

-.0000 

-.0000 

.0000 

1411 

.0000 

-.0000 

.0000 

CONSTANT 

504.5424 

503.2446 

505.3401 

15EC0ND  ORDER  REGRESSION  FOR  4  INDEP.  VAR  94/03/31.  15.56.17.  PAGE  13 
FILE  -  REGRE3SI  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGDAT 

*********MULTIPLE  REGRESSION********* 
DEP.  VAR...  CF 


VARIANCE/COVARIANCE  MATRIX  OF  THE  UNNQRMAl  I ZED  REGRESSION  COEFFICIENTS. 


W1 

.00000 

142 

,00000 

.00000 

W3 

.00000 

,00000 

,00000 

144 

.00000 

,00000 

.00000 

.00000 

Wtl 

-.00000 

-.00000 

-.00000 

-.00000 

.00000 

1422 

-.00000 

-.00000 

-.00000 

- . 00000 

,00000 

.00000 

W33 

-.00000 

-.00000 

-.00000 

-.00000 

.  00000 

,00000 

!444 

-.00000 

-.00000 

-.00000 

-.00000 

.00000 

.00000 

W13 

-.00000 

-.00000 

-.00000 

-.00000 

.00000 

.00000 

1414 

-.00000 

-.00000 

-.00000 

-.00000 

.00000 

.00000 

U23 

-.00000 

-.00000 

-.00000 

-.00000 

,00000 

.00000 

1424 

-.00000 

-.00000 

-.00000 

1 

♦ 

o 

o 

o 

o 

,00000 

.00000 

W34 

-.00000 

-.00000 

-.00000 

-.00000 

,00000 

,00000 

W33 

.00000 

W44 

,00000 

. ooooo 

WIT 

.00000 

. OOOOO 

» ooooo 

W14 

. 00000 

.00000 

'.00000 

.00000 

W23 

.00000 

.00000 

-.00000 

-.00000  .00000 

W24 

.00000 

.00000 

',00000 

.00000  -.00000 

. ooooo 

U34 

.00000 

.00000 

-.00000 

-.00000  -.00000 

.00000 

W33 

W44 

W13 

W14  U23 

W24 

W34  .00000 

W34 

1 SECOND  ORDER  REGRESSION  FOR  4  INDEP.  VAF;  84/03/31.  15.56.17.  PAGE  14 
FILE  -  REGRESSI  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGB AT 

MMUinNULTIPLE  REGRESSION******** 
DEP.  VAR...  CF 


SUMMARY  TABLE. 


STEP 

VARIABLE 

E/R 

F  MULT-R 

R-SQ  CHANGE 

R 

OVERALL  F 

S1G . 

1 

W12 

E 

500.498 

,978 

.956 

.956 

.978 

500,498 

0 

W13 

E 

44,88/ 

.993 

.986 

,02? 

,601 

750,201 

,000 

3 

W14 

E 

14.303 

.996 

,992 

.006 

.408 

318,861 

0 

4 

Wll 

E 

5.995 

.997 

,993 

,002 

.733 

761,71? 

.000 

5 

W2 

E 

53.057 

.99? 

.998 

,005 

,  669 

2206.100 

0 

6 

Ml 

E 

31,552 

1,000 

,999 

,001 

.735 

4799,383 

.000 

7 

W12 

R 

,000 

1,000 

,999 

0 

,978 

6079,851 

0 

8 

W22 

E 

54,029 

1,000 

1.000 

,000 

.665 

19216.17? 

.000 

9 

Wll 

R 

,002 

1.000 

1.000  - 

-.000 

.733 

24337.452 

0 

10 

W3 

E 

39,242 

1,000 

1,000 

.000 

.088 

61103.470 

,000 

11 

U4 

E 

3411.581 

1.000 

1.000 

.000 

.060 

.998EF07 

0 

12 

W34 

E 

6,890 

1.000 

1,000 

.000 

.093 

, 118E+08 

,000 

13 

W23 

E 

5,284 

1.000 

1.000 

.000 

.591 

,  132EF03 

0 

14 

W33 

E 

5.359 

1.000 

1.000 

,000 

,089 

.  154ET03 

.000 

15 

W24 

E 

2,459 

1,000 

1.000 

.000 

.429 

, 154E+0S 

o 

16 

W44 

E 

1,479 

1,000 

1.000 

.000 

.060 

, 147EP08 

,000 

17 

Mil 

E 

.046 

1.000 

1.000 

.000 

,733 

, 125E+08 

0 

SECOND  ORDER 

REGRE 

S3 I ON  FOR 

4  INDEP.  VAR 

34/03/31. 

15.56,17. 

PAGE 

FILE  -  REGRESS!  (CREATED  -  84/03/31)  ANALYSIS  WITH  RGDAT 


#********multiple  regression********* 


RESIDUAL  PLOT . 


Y  VALUE 

Y  EST. 

RESIDUAL 

-2SD 

0.0 

+280 

932.670 

932.655 

.015 

I 

868.400 

368.403 

-.003 

.  I 

862.120 

862.121 

-.001 

.1 

797.850 

797.370 

-.020 

♦ 

I 

873.840 

873.855 

-.015 

♦ 

I 

868.130 

363.126 

.004 

I  . 

365.410 

365.413 

-.008 

I 

859.610 

859.600 

.010 

I 

904.920 

904.912 

.008 

I 

399.150 

899.153 

-.003 

.  I 

934.400 

334.393 

.007 

I 

823.600 

828.605 

- .  005 

.  I 

901.590 

901.595 

-.005 

.  I 

393- 140 

893.148 

-.008 

+ 

I 

337.390 

837.373 

.012 

I 

828.370 

328.362 

.008 

I 

906.240 

906.247 

-.007 

» 

I 

897.790 

397.800 

-.010 

¥ 

I 

335.760 

835.749 

.011 

I 

827.240 

327.232 

.008 

I 

900.260 

900.257 

.003 

I  , 

894.500 

894.504 

-.004 

.  I 

336.030 

836.025 

.005 

I 

330.230 

830.232 

-.002 

.  I 

366.770 

866.770 

.000 

♦ 

NOTE  -  (*) 

INDICATES 

ESTIMATE 

CALCULATED 

WITH 

MEANS  SUBSTITUTED 

R 

INDICATES 

POINT  CUT 

OF  RANGE  OF 

PLOT 

NUMBER  OF  CASES  PLOTTED  25. 

NUMBER  OF  2  S.D.  OUTLIERS  0  OR  0  PERCENT  OF  THE  TOTAL 

VON  NEUMANN  RATIO  1.79113  DURBIN-WATSON  TEST  1.71754 

NUMBER  OF  POSITIVE  RESIDUALS  12. 

NUMBER  OF  NEGATIVE  RESIDUALS  13, 

NUMBER  OF  RUNS  OF  SIGNS  15, 

EXPECTED  NUMBER  OF  RUNS  OF  SIGNS 
EXPECTED  S.D.  OF  RUN  DISTRIBUTION 
UNIT  NORMAL  DEVIATE- 

Z= <  OBSERVED-EXPECTED+0 , 5 > /S . D . 

PROBABILITY  OF  OBTAINING  ,GE,  ABSCZ) 


111 


13. 

2.44328 

.82676 

.20419 


CF  Summary  Regression 


:*  *  %  X  t 

*  #  *  M 

U  1.  T  I  P 

L  E  RE  3 

R  E 

S  S  *  0  N  t  % 

%  s:  %  3*  * 

REP.  VAR , » ♦ 

CV 

VARIABLE (S) 

ENTERED 

ON  STEP  1 

W4 

MULTIPLE  R 

1.0000 

ANOVA 

DF  SUM 

SQUARES  MEAN  SQ, 

F 

R  SQUARE 

.9999 

REGRESSION 

12, 

4003, 

743  333.646 

13552.390 

STD  DEV 

,1569 

RESIDUAL 

12. 

* 

295  .025 

SIG,  .000 

ADJ  R  SQUARE  ,999? 

COEFF  OF  VARIABILITY 

.  OPCT 

VARIABLE 

B 

S.E.  B 

F 

SIG, 

BETA  E 

LAST I CITY 

W12 

-.000 

.000 

3,416 

,039 

-.05103 

-.00249 

142 

.120 

.006 

350 , 895 

,000 

,93394 

,06171 

Ml 

,  166 

,014 

133 . 4?9 

,000 

.68034 

,  05338 

143 

,031 

.007 

17,845 

.001 

.26941 

,02323 

U 13 

-.000 

.000 

/  ♦  1 6  5 

,020 

- ,08571 

-.00496 

1423 

-.000 

.000 

7,165 

.020 

■■ » 07346 

-  ,00397 

W33 

-.000 

,000 

4,195 

.063 

-.11542 

-.00504 

1444 

.000 

,000 

,523 

.484 

, 02076 

.00048 

1434 

-.000 

.000 

1.545 

.238 

-.03305 

-.00123 

W22 

-.000 

.000 

.929 

,354 

-.03958 

-.00127 

1411 

-.000 

.000 

.353 

.374 

-.04734 

-.00139 

144 

.002 

,006 

,099 

.758 

.01169 

.00051 

CONSTANT 

1027,03? 

7.121 

20799.923 

0 

F-LEVEL  OR  TOLERANCE-LEVEL  INSUFFICIENT  FOR  FURTHER  COMPUTATION  * 
COEFFICIENTS  AND  CONFIDENCE  INTERVALS ► 


VARIABLE 

B 

95  F'CT 

C.I. 

1412 

-.0000 

-.0000 

,0000 

142 

,1198 

.1060 

» 1336 

Wt 

,1653 

,1346 

,1971 

143 

,0308 

,0149 

.0466 

1413 

-.0000 

-.0000 

-.0000 

1423 

-.0000 

-.0000 

-.0000 

1433 

-.0000 

-.0000 

,0000 

1444 

.0000 

-.0000 

,0000 

W34 

-.0000 

-.0000 

,0000 

1422 

-.0000 

-.0000 

,0000 

1411 

-.0000 

-.0001 

,0000 

144 

,0019 

-.0113 

,0151 

CONSTANT 

, 10E+04 

. 10E+04 

,  IOEK'4 

112 


RESIDUAL  PLOT 


Y  VALUE 

Y  EST, 

RESIDUAL  -2SD 

0.0 

1190,020 

1190.021 

-.001 

» 

1159,590 

1159.565 

.025 

I  , 

1170,660 

1170. 633 

.027 

I  . 

1139,650 

1139,596 

,054 

I 

1165.560 

1165,618 

-.058 

» 

T 

X 

1165.560 

1165.706 

-.146 

* 

]■ 

1.164,700 

1164.501 

,1.99 

I 

1164,310 

1164,200 

,110 

I 

1174,870 

1175.004 

-.134 

* 

I 

1174,810 

1174,897 

-.087 

♦ 

T 

1155,390 

1155.325 

,065 

I 

1155.320 

1155,219 

.101 

I 

1180.830 

1180.663 

.162 

I 

1179.810 

1179,776 

.034 

r 

X  » 

1150.230 

1150.341 

-.061 

♦ 

I 

1143.420 

1148,610 

-.190 

♦ 

I 

1175.300 

1175.137 

.163 

I 

1174.230 

1.174,246 

.034 

I  . 

1155.320 

1155,879 

-.059 

♦ 

I 

1153,960 

1154.147 

-.187 

♦ 

V 

X 

1130,400 

1130.534 

-.134 

♦ 

I 

1180.340 

1180,428 

-.080 

♦ 

1 

1149.350 

1149.783 

.062 

I 

1149,790 

1149.681 

,109 

I 

1165.130 

1165.130 

.000 

* 

MOTE  -  (SO 

INDICATES 

ESTIMATE  CALCULATED  WITH  ME 

ANS  SUBSTITUTED 

R  INDICATES  POINT  OUT  OF  RANGE  OF  PLOT 


NUMBER  OF  CASES  PLOTTED  25, 

NUMBER  OF  2  S.D.  OUTLIERS  0  OR  0  PERCENT  OF  THE  TOTAL 

VON  NEUMANN  RATIO  1,71919  DURBIN- WATSON  TEST  1,650 


‘•iW  T*  .■»  TTTTTJT?  '.'■  'XV  V  ’V/V.V;?. 


VV, 

& 


AF  Summary  Regression: 


mEP,  VAR...  AF 

VARIABLE (S)  ENTERED  ON  STEP  14 
W44 


MULTIPLE  R 

.999? 

ANOVA 

DF  SUM 

SQUARES 

MEAN  SQ, 

F 

R  SQUARE 

.9999 

REGRESSION 

14,  9 

7163,036 

6940.220 

6349 ,7' 

STD  DEV 

1.0066 

RESIDUAL 

10, 

10.132 

1,013 

SIG.  ,0' 

ADJ  R  SQUARE 

.999? 

COEFF  OF  VAR 

IABILITY 

.  1PCT 

VARIABLE 

B 

S.E.  B 

F 

SlOi. 

BETA  El 

..AST IC  [1 ' 

W12 

-.001 

,000 

238,497 

.000 

-.55530 

-.09529 

W34 

-.000 

.  000 

17,245 

,002 

-.14332 

-.01943 

W3 

CTI 
«•  yj 

.04? 

125. 163 

,000 

,23190 

W2 

1.335 

,  042  1 

028 . 334 

,000 

2,22576 

,491.11 

ui 

1.357 

,094 

203,661 

,000 

X , 13109 

.31196 

W22 

-.000 

.000 

110,893 

.000 

- . 56309 

-  ,06373 

U23 

-.000 

.000 

144.143 

,000 

-.45330 

-.03149 

W4 

.435 

.048 

81.119 

,000 

.54150 

,08321 

U24 

-.000 

,000 

59,509 

,000 

-.21173 

-.02652 

W13 

-.000 

.  000 

8,853 

.014 

-.12407 

-.02524 

U33 

-.000 

,000 

7.355 

,019 

-  .20606 

-.03169 

Ull 

-,000 

,000 

3,723 

.083 

-.12373 

-.01305 

U14 

-.0000 

,  656 

,43? 

-.02605 

-.00343 

W44 

-.000 

,000 

,250 

♦  o23 

-.01370 

-.00151 

CONSTANT  3 

23.743 

46,559 

48.351 

,000 

ALL  VARIABLES  ARE  IN  THE  EQUATION, 
COEFFICIENTS  AND  CONFIDENCE  INTERVALS, 


VARIABLE 


B 


95  F’CT  C,I, 


W12 

-.0007 

-.0008 

-.0006 

W34 

-,0001 

-.0002 

-.0001 

W3 

* i.  < 

.4185 

,  6267 

W2 

1,3353 

1.2425 

1,4281 

Ul 

1.3571 

1,1478 

1.5665 

W22 

-.0003 

-.0003 

-.0002 

W23 

-.0003 

-.0003 

1 

o 

o 

o 

ro 

W4 

.4351  , 

5427 

W24 

-.0002 

-.0003 

-.0002 

W13 

-.0001 

- >0002 

-.0000 

U33 

-.0001 

-.0001 

-.0000 

Wll 

-.0002 

-.0004 

,0000 

W14 

-.0000 

-.0002 

.0001 

W44 

-.0000 

-,0001. 

,0001 

CONSTANT 

323.7475 

220,0071  4 

27.4379 

SECDND  ORDER 

REGRESSION 

FOR;  4  TNDEP 

,  VAR 

114 


r-'AGF 


TABLE  VI 


CF  Regression  Error 


continued 


Table  VI  continued 


Random  Sample 


441 

469 

803 

226 

864.41 

864.289 

.121 

328 

458 

752 

286 

808.84 

808.915 

.075 

391 

604 

923 

333 

876.80 

876.780 

.021 

300 

645 

969 

260 

841.71 

841.990 

.280 

423 

711 

853 

357 

912.46 

912.689 

.228 

350 

480 

825 

275 

826.08 

825.885 

.195 

429 

577 

726 

208 

880.68 

880.385 

.294 

403 

585 

789 

241 

872.62 

872.423 

.196 

408 

470 

947 

410 

857.71 

857.549 

.160 

489 

545 

791 

216 

904.31 

903.527 

.783 

355 

462 

894 

392 

829.02 

828.921 

.099 

Random  Sample  Mean  Error,  5  =  0.223 

Random  Sample  Standard  Deviation  of  Error  =  0.204 


TABLE  VII 


CV  Regression  Error 


Design  Points 


W1 

W2 

W3 

W4 

cv 

m 

cve 

Residual 

450 

750 

880 

312 

1190.02 

1190.021 

.001 

450 

450 

880 

312 

1159.59 

1159.565 

.025 

300 

750 

880 

312 

1170.66 

1170.633 

.027 

300 

450 

880 

312 

1139.65 

1139.596 

.054 

375 

600 

1040 

424 

1165.56 

1165.618 

.058 

375 

600 

1040 

200 

1165.56 

1165.706 

.146 

375 

600 

720 

424 

1164.70 

1164.501 

.199 

375 

600 

720 

200 

1164.31 

1164.200 

.110 

450 

600 

880 

424 

1174.87 

1175.004 

.134 

450 

600 

880 

200 

1174.81 

1174.897 

.087 

300 

600 

880 

424 

1155.39 

1155.325 

.065 

300 

600 

880 

200 

1155.32 

1155.219 

.101 

375 

750 

1040 

312 

1180.83 

1180.668 

.162 

375 

750 

720 

312 

1179.81 

1179.776 

.034 

375 

450 

1040 

312 

1150.28 

1150.341 

.061 

375 

450 

720 

312 

1148.42 

1148.610 

.190 

450 

600 

1040 

312 

1175.30 

1175.137 

.163 

450 

600 

720 

312 

1174.28 

1174.246 

.034 

300 

600 

1040 

312 

1155.82 

1155.879 

.059 

300 

600 

720 

'312 

1153.96 

1154.147 

.187 

375 

750 

880 

424 

1180.40 

1180.534 

.134 

375 

750 

880 

200 

1180.34 

1180.428 

.088 

375 

450 

880 

424 

1149.85 

1149.788 

.062 

375 

450 

880 

200 

1149.79 

1149.681 

.109 

375 

600 

880 

312 

1165.13 

1165.130 

.000 

Design  Point  standard  Deviation  of  Error  =  0.1569 


continued 


Random 

Sample 

441 

469 

803 

226 

1159.99 

1159. 

955 

.035 

328 

458 

752 

286 

1143.91 

1143. 

285 

.625 

391 

604 

923 

333 

1167.73 

1167. 

811 

.081 

300 

645 

969 

260 

1160.21 

1160. 

778 

.568 

423 

711 

853 

357 

1182.59 

1182. 

649 

.059 

350 

480 

825 

275 

1149.49 

1149. 

118 

.372 

429 

577 

726 

208 

1169.02 

1169. 

070 

.050 

403 

585 

789 

241 

1166.83 

1167. 

267 

.437 

408 

470 

947 

410 

1156.36 

1156. 

467 

.107 

489 

545 

791 

216 

1173.88 

1174. 

003 

.123 

355 

462 

894 

392 

1148.52 

1148. 

406 

.114 

Random  Sample  Mean  Error,  5  =*  0.234 

Random  Sample  Standard  Deviation  of  Error 

=  0.222 

TABLE  VIII 

AF  Regression  Error 


Design  Points 

W1 

w2 

w3 

w4 

AFm 

AFe 

Residual 

450 

750 

880 

312 

1736.46 

1737.320 

.860 

450 

450 

880 

312 

1618.72 

1619.040 

.310 

300 

750 

880 

312 

1652.86 

1653.337 

.477 

300 

450 

880 

312 

1504.04 

1503.967 

.073 

375 

600 

1040 

424 

1677.36 

1677.435 

.075 

375 

600 

1040 

200 

1645.10 

1645.622 

.522 

375 

600 

720 

424 

1626.09 

1626.355 

.265 

375 

600 

720 

200 

1585.47 

1586.182 

.712 

450 

600 

880 

424 

1701.47 

1701.778 

.308 

450 

600 

880 

200 

1668.05 

1666.600 

1.450 

300 

600 

880 

424 

1602.27 

1603.065 

.795 

300 

600 

880 

200 

1567.22 

1566.256 

.964 

375 

750 

1040 

312 

1717.52 

1716.390 

1.130 

375 

750 

720 

312 

1673.61 

1673.215 

.395 

375 

450 

1040 

312 

1594.91 

1594.650 

.260 

375 

450 

720 

312 

1526.83 

1527.305 

.475 

450 

600 

1040 

312 

1708.49 

1708.939 

.449 

450 

600 

720 

312 

1657.15 

1656.674 

.476 

300 

600 

1040 

312 

1612.06 

1612.405 

.345 

300 

600 

720. 

312 

1554.73 

1554.150 

.588 

375 

750 

880 

424 

1710.86 

1710.299 

.581 

375 

750 

880 

200 

1681.32 

1682.070 

.750 

375 

450 

880 

424 

1585.12 

1584.239 

.881 

375 

450 

880 

200 

1540.05 

1540.480 

.430 

375 

600 

880 

312 

1635.88 

1635.880 

.000 

Design  Points  Standard 

Deviation  of  Error  *  1 

.0066 

continued 
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Table  VIII  continued 


Random 

Sample 

441 

469 

803 

226 

1591.84 

1588.798 

3.042 

328 

458 

752 

286 

1494.51 

1496.747 

2.237 

391 

604 

923 

333 

1658.73 

1658.556 

.170 

300 

645 

969 

260 

1612.72 

1613.919 

.199 

423 

711 

853 

357 

1711.99 

1711.877 

.113 

350 

480 

825 

275 

1541.90 

1540.694 

1.206 

429 

577 

726 

208 

1614.18 

1614.414 

.234 

403 

585 

789 

241 

1619.45 

1619.095 

.355 

408 

470 

947 

410 

1626.35 

1628.740 

2.390 

489 

545 

791 

216 

1654.74 

1654.590 

.150 

355 

462 

894 

392 

1571.61 

1571.837 

2.770 

Random  Sample  Mean  Error <  1  =  1.233 
Random  Sample  Standard  Deviation  of 

Error  ■  1 

.154 

APPENDIX  D 


Sequential  Unconstrained  Minimization  Technique 

(SUMT) 

The  SUMT  program  was  developed  by  Mylander,  et  al  (23) 
based  on  the  nonlinear  programming  work  of  Fiacco  and 
McCormick  (12).  A  descriptive  pamphlet  of  the  program  was 
developed  by  Captain  R.  M.  Floyd  for  the  Air  Force  Institute 
of  Technology.  With  his  permission,  this  pamphlet  is  re¬ 
produced  in  this  appendix.  Following  the  pamphlet  is  an 
example  of  both  an  executive  program  that  was  used  in  the 
study  and  a  sample  output. 
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A .  Purpose 


This  program 

function  subject 

finds  the  minimum  of  a  multivariable,  nonlinear 

to  nonlinear  inequality  and  equality  constraints:. 

Minimize 

F(X1,  X2,  ....  XfJ) 

Subject  to 

k  ■  1,  2,  ...,  M 

W  X2 . V  *  °* 

k  »  Mel,  M*2,  ...,  M*MZ 

B.  Method 

The  procedure  was  developed  by  Fiacco  and  McCormick  (1).  The 
technique  uses  the  problem  constraints  and  the  original  objective  function 
to  form  an  unconstrained  objective  function  which  is  minimized  by  any 
appropriate  unconstrained,  multivariable  technique. 

The  algorithm  proceeds  as  follows: 

1)  A  modified  objective  function  formulated  consisting  of  the 
original  function  and  penalty  functions  with  the  form 

M  M-MZ 

P  .  F  -  r  E  $nG  *  Z  HT/r 
feat  K  k-M+l  K 

where  r  is  a  positive  constant.  As  the  algorithm  progresses,  r  is 
reevaluated  to  form  a  monotonically  decreasing  sequence  r^>r2>...  '>0- 
As  r  becomes  small,  under  suitable  conditions  P  approaches  F  and  the 
problem  is  solved. 

2)  Select  a  starting  point  (feasible  or  nonfeasible)  and  an  initial 
value  for  r. 

3)  Determine  the  minimum  of  the  modified  objective  function  for  the 
current  value  of  r  using  an  appropriate  technique  (several  options 
available). 
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4)  Estimate  the  optimal  solution  using  extrapolation  formulas  (1), 

5)  Select  a  new  value  for  r  (1)  and  repeat  the  procedure  until 
the  convergence  criterion  is  satisfied. 

A  logic  diagram  for  this  method  is  given  in  Figure  1. 

Use  of  the  SUNT  Library 

Execution  of  the  SUMT  routines  requires  that  the  user  submit 
a  computer  job  composed  of  the  following  elements i  1)  Job 
Control  commands,  2)  a  MAIN  program,  3)  user-supplied  sub¬ 
routines,  and  4)  NAMELIST  data.  The  layout  of  these  elements 
is  shown  below. 

JOB  CARDS 

(7/8/9) 

MAIN  PROGRAM 

USER-SUPPLIED  SUBROUTINES 

(7/8/9) 

NAMELIST  DATA 
(6/7/S/9) 

Note  that  an  end-of-record  and  an  end-of-file  is  signified 
by  a  multipunch  (7/8/9)  and  (6/7/8/9),  respectively  in  column 
1  of  a  keypunch  card  (or  by  *EOR  and  *EOF,  respectively  when 
coded  at  a  computer  terminal).  Each  element  is  described 
separately  below,  in  Section  E  an  example  is  presented  which 
combines  these  elements  for  job  execution. 

1)  Job  Cards 

The  Job  Card  sequence  required  to  execute  SUMT  is 

abc,Tt,IOk,CM100000.  PN.name 
ATTACH , lfn , SUMT , ID;T800535 . 

LIBRARY, If n. 

FTN. 

LGO. 


in  which 


3  character  identifier 

decimal  seconds  of  Central  Processor  time 

requested  for  the  job  (default  is  20) 

decimal  seconds  of  Input/Output  time  requested 

for  the  job  (default  is  70) 

user  computer  problem  number 

user  r.ar.e  or  ether  identifying  information 

2. C C 3. 1  fil,6  7*3-716  1*01*  thS  2 U 1 T  2  i *7 Z'r'. T'r 


\- v  vV 


Correct  Starting 
Point  Until 
Feasible 


Fiacco  and  McCormick  (SUMT  ALGORITHM) 
Logic  Diagram 


2)  MAIN  Program 

The  MAIN  program  identifies  three  files  and  assigns  corres¬ 
ponding  device  numbers  in  a  Program  statement.  It  communi¬ 
cates  the  device  numbers  through  the  labelled  COMMON i 

COMMON/DEVC/NI ,N0 ,NS 

The  variables  NI,  NO,  NS  are  given  values  by  MAIN  and  cor¬ 
respond  to  the  device  numbers  for  the  input  and  output  files 
and  a  save  file,  respectively.  Execution  of  the  SUMT  routines 
is  initiated  by  a  call  to  subroutine  SUMT.  The  code  given 
below  demonstrates  the  requirements  of  the  MAIN  routine* 

PROGRAM  MAIN (INPUT, OUTPUT, SA VS, TA?E5rINPUT, 

TAPE61OUTPUT . TAPE7;SAVE ) 

COMMON/DEVC/NI . NO , NS 

NIT  5 

N0=6 

NS:7 

CALL  SUMT 

STOP 

END 

3)  User-Supplied  Subroutines 

The  subroutines  described  below  are  called  by  the  given  names 
and  with  the  given  argument  lists.  Inclusion  of  the  routine 
RESTNT  is  mandatory,  while  routines  GRAD1  and  MATRIX  are 
optional  (default  routines  are  coded  in  SUMT).  Each  user- 
supplied  subroutine  must  contain  the  COMMON  cardi 

COMMON/S HARE/X ( 100) ,DEL(100 ) ,A (100,100) ,N(5) 

This  cpmmon  region  is  used  for  communication  between  the 
SUMT  routines  and  the  user-supplied  routines.  For  each 
routine,  the  argument  IN  is  an  input  argument  referring  to 
the  objective  function,  the  inequality  constraints,  or  the 
equality  constraints  as  follow* 

IN  r  0  1  objective  function  is  being 

considered 

IN  Z  1,2 , . . . ,M  1  the  INth  inequality  constraint 

is  being  considers! 

IN  I  M+l,M+2 . M+MZ  :  the  1st,  2nd,...,  Mlth  equality 

constraint  is  being  considered 


RESTNT(IN.VAL) 

This  subroutine  sets  the  argument  VAL  to  the  value  of  the 
objective  function,  an  inequality  constraint,  or  an  equality- 
constraint  according  to  the  value  of  the  argument  IN.  The 
current  values  of  the  independent  variables  are  found  in  the 
vector  X  of  the  labelled  COMMON/SHARE/. 


GRADl(IN) 

This  subroutine  computes  the  gradient  of. the  objective 
function,  an  inequality  constraint,  or  an  equality  constraint 
according  to  the  value  of  the  argument  IN.  The  gradient 
with  respect  to  the  Ith  independent  variable  is  3tored  in 
location  DEL(I).  The  current  values  of  the  independent 
variables  are  in  the  vector  X,  and  both  X  and  DEL  are  in 
labelled  COMMON/SHARE/.  The  elements  of  QSL  are  not  set  to 
zero  prior  to  entry  into  GRAD1.  The  routine  GRAD1  included 
by  default  in  SUMT  computes  a  numerical  approximation  for 
each  gradient  using  the  central  difference  technique.  If 
the  user  wishes  to  provide  some  of  the  gradients  analytically 
and  others  by  numerical  differencing,  the  central  difference 
approximation  for  any  desired  gradient  as  indexed  by  the 
argument  IN  may  be  computed  using  a  call  to  the  SUWT  library 
routine  DIFF1  as 


CALL  DIFFl(IN) 


with  selective  execution  according  to  the  value  of  IN.  In 
the  usual  case,  user's  rely  on  the  default  GRAD1  routine 
rather  than  provide  a  replacement  for  it  with  analytically 
determined  gradients. 


MATH IX ( IN. L) 


This  subroutine  computes  the  upper  triangle  and  diagonal 
elements  of  the  matrix  of  second  partial  derivatives  of  the 
objective  function,  an  inequality  constraint,  or  an  equality 
constaint  according  to  the  value  of  the  argument  IN.  The 
current  values  of  the  independent  variables  are  in  the  vector 
X,  and  the  computed  record  partials  are  stored  ir.  array  A. 


of  the  array  A  below  the  main  diagonal  must  not  be  modified 
within  MATRIX i  the  upper  triangle  of  A  is  set  to  zero  before 
entry  into  MATRIX.  The  second  argument  (L)  is  provided  so 
that  the  user  may  inform  SUMT  that  the  matrix  of  second 
partial  derivatives  of  a  function  is  zero«  L  should  be  set 
to  1  according  to  the  value  of  argument  IN  if  the  corres¬ 
ponding  matrix  of  second  partials  is  zero.  The  routine  MATRI' 
included  by  default  in  SUMT  computes  a  numeriaal  approx¬ 
imation  for  each  matrix  of  second  partial  derivatives  using 
central  differences  to  second  order.  If  the  user  wishes  to 
provide  some  of  the  matrices  analytically  and  others  by 
numerical  differencing,  the  central  difference  for  any  des¬ 
ired  matrix  of  second  partials  as  indexed  by  the  argument  IN 
may  be  computed  using  a  call  to  the  SUMT  library  routine 
DIFF2  as 

CALL  DIFF2 (IN) 


with  selective  execution  according  to  the  value  of  IN.  In 
the  usual  case,  user's  rely  on  the  default  MATRIX  routine 
rather  than  provide  a  replacement  for  it  with  analytically 
determined  second  partial  derivative  matrices. 


NAMELIST  Data 

Data  is  entered  into  SUMT  using  a  NAMELIST  read  from  the 
input  file)  the  NAMELIST  is  referred  to  as  DATA  and  includes 
the  following  variables* 


NAMELIST/DATA/EPS I , RHOIN .THETAO , RATIO , TMMAX ,M , N ,MZ , 
X ,NT ,XSP1,XE?2,NEX0P1,NEX0P2 


The  variables  X  and  NT  are  vectors  of  length  100  and  10, 
respectively.  All  of  the  variables  in  the  NAMELIST  are 
assigned  default  values  by  SUMT  (see  Section  ^b  below).  To 
assign  different  values  to  any  subset  of  these  variables  an 
assignment  is  made  to  each  appropriate  variable  name  in  the 
NAMELIST  input.  The  complete  data  set  is  delimited  on  the 
left  by  $D ATA  and  on  the  right  by  $ZND ;  individual  data 


assignments  are  separated 


ecnr.nsi  t.ne 
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over  one  or  more  80  column  fields.  No  information  should 
be  placed  in  column  1  of  any  data  field,  since  the  NAMELIST 
read  processes  information  beginning  with  column  2.  Refer 
to  an  appropriate  FORTRAN  manual  for  more  information  on 
NAMELIST  input. 

a.  Definition  of  NAMELIST  Data  Items 

EPSI  Tolerance,  €,  used  to  decide  if  an  unconstrained  minimum 
has  been  achieved  for  each  subproblao  (seeNT(7)) 

8HOIN  Possible  initial  value  of  r,  r  ,  (often  set  at  1.0)  (see  N T  (1)) 

THETAO  Tolerance,  9o,  used  to  decide  if  the  solution  to  the  problem 
has  been  approximated  (see  NT  (5)) 

RATIO  Parameter,  C,  (>l)  used  to  compute  consecutive  values  of  r; 
ri+  a  r^/C  (often  set  at  16.0) 

HIM  AX  Maximum  amount  of  time  for  solving  problem  (sec.) 

M  Number  of  inequality  constraints 

N  Number  of  independent  variables 

MZ  Number  of  equality  constraints 

X  Independent  variables  (starting  point  on  input) 

NT(1)  Option  key  for  r  values,  as  follows: 

»1  The  value  for  r  is  made'  bv  finding  an  approximation 
solution  min{Crp(X°,  r)  !>5P(X?,  r)  J-^CX0,  r)  ] ]  which 
is  a  good  approximation  only  when  X°  is  close  to  the 
boundary  (i.e.,  for  some  i,  Gj^X)  is  close  to  zero)  or 
when  V2F(X°)aO  and  when  MZaO 

>2  rt  is  given  by  formula  8.65  [Ref.  (l)p.  191  .]  (Only 

can  be  used  when  MZaO) 

-3  rt  -  RHOIN 
(normally  set  to  3) 

NT (2)  Option  key  for  constraints  as  follows: 

•  1  The  requirements  (trivial  constraints)  that  *  0 
for  i  m  1,  ...,  N  are  to  be  automatically  included 
in  the  problem 

•2  The  only  constraints  on  the  problem  are  those  inputed 
by  the  user 

NT  (3)  Option  key  for  printout  as  follows: 

■  1  Standard  printout  (this  includes  a  call  to  OUTPUT  after 
the  solution  of  overy  subproblem.  Also  the  estimates 
of  the  "Lagrange  multipliers"  and  first-  and  second- 
order  solution  estimates  are  printed) 

•2  For  additional  printout  (includes  standard  printout 
and  every  iniet.vedi.ite  point,  gradient  of  P  and  the 
vector  S) 

(normally  sot  to  :) 

;i  f  (L )  Option  'toy  or:  firm.'.  t  c  once  ns  follows: 

«1  Pinal  cor.veracr.1.0  is  ditc-r.iunod  on  the  basis  oi  irrent 
solution  to  the  .vabpreb!  ■  m 
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NT  (6) 


Final  convergence  is  determined  on  the  basis  of  the 
first  order  estimates.  The  first  order  estimate  of 
the  solution  vector  must  be  close  to  feasible. 

See  below. 


>3  Final  convergence  is  determined  on  the  basis  of  the 
second  order  estimates.  The  second  order  estimate  of 
the  solution  vector  must  be  close  to  feasible  before 
the  convergence  check  is  made.  If  X  is  a  solution 
estimate  it  is  considered  close  to  being  feasible  if 
dt(X)  .  Mo.i.l,  J,  ....  M. 

(normally  set  to  l) 


Option  key  on  final  convergence  as  follows: 

The  convergence  criterion  determining  the  problem  has  bden 
solved  (Only  use  *  1,  whenN T  (** U  l) 

....  G  -  F  Cx(rw)] 

-1  *** wh,n  siYCr  TO-xxfrn <  9o 


■2  Quit  when  |r  2  .*  in  G  ,[x(r.  )  3 1  <  8  (use  only  fbr  M  =/» 
J* .  J  k  o 


Quit  when 


first  order  estimate  of  v 


GLXlr^) .^(r^) .XCr^) J 


-  1  <  8 


(normally  set  to  l) 

Option  key  for  extrapolation  as  follows: 


No  extrapolation 


■2  Extrapolate  through  last  two  minima 


*3  Extrapolate  through  last  three  minima 
(normally  set  to  l) 


NT  (7) 


Option  key  for  subproblem  convergence  (fixed  value  of  r) 
as  follows:  r  -.2 .-i-i  . 

-1  Quit  when  ^JP^X^.r)  7«P(Xl,r)|  <  c 


follows:  r  Tfiofr  il*"1  i 

Quit  when  (X^.r)  ax'^x'l  7XP(X^r)  I 

L  i  j  J 

Quit  when  jv^tX^r)  1  7xp(xi*r)l 

<  PtX1"1)  -  PCX1) 


.3  Quit  when  |’7xP(X1,r)  !  <  € 


NT  (8) 


Option  key  on  problem  linearity  as  follows: 
al  At  least  one  nonlinear  constraint 


Linear  constraints 


■3  Linear  constraints  and  linear  objective  function 
(i.e. ,  a  linear  programming  problem) 


When  Option  -  *  3  MATRIX  Uhe  uj«r  subroutine  'implying 
the  second  cj.  nai  derivatives)  will  not  be  caiied,  and 
when  Option  d  *  2  it  will  be  called  only  to  get  the  second 
partials  of  F(X). 


v.v.v.v.v 


NT  (9)  not  used 


NT(10) 

XEP1 

XEP2 


NEXOPl 


NEXOP2 


not  used 


Finite  difference  parameter  used  in  DIFF1  if  numerical 
derivatives  used.  Usually  XEPl  >  O.OOOl  is  satisfactory. 

Iteration  improvement  limit.  When  minimizing  the  P-function 
for  a  given  value  of  r  (XHO)  the  value  of  P  must  decrease 
by  an  amount  exceeding  XEP2  for  each  iteration  after  the 
first.  If  it  does  not,  then  the  code  prints  out  the 
message  "apparently  roundoff  errors  prevent  a  more  accurate 
determination  of  the  minimum  of  this  subproblem,"  and  it 
is  assumed  a  minimum  has  been  found.  (Normally,  set  XEP2 
equal  to  0. ) 

Key  for  checking  derivatives  as  follows: 

■1  Solve  problems  without  checking  derivatives 

m2  Solve  problem  after  checking  only  first  derivatives 

•3  Do  not  solve  ■problem  after  checking  only  first 
.  derivatives 

■4  Solve  problem  after  checking  first  and  second 
derivatives 

■5  Do  not  solve  problem  after  checking  first  and  second 
derivatives 

If  the  user  wishes  to  check  his  first  or  second  partial 
derivatives  by  having  them  printed  along  with  experimental 
SUMT's  numerical  approximations,  he  can  accomplish  this 
by  appropriately  setting  NEXOPl. 

In  all  cases  the  first  values  printed  out  are  those  gotten 
from  the  user's  code,  end  the  second  values  are  those 
gotten  by  numerical  differencing. 

Key  for  choosing  unconstrained  minimization  technique  as 
follows: 

■1  The  method  for  minimizing  the  unconstrained  penalty 
function  is  to  be  the  generalized  Newton-Raphson 
method  as  modified  to  .landle  indefinite  Hessian 
matrices.  This  method  requires  function  values,  first 
and  second  derivatives. 

■2  Same  as  1,  except  that  when  an  "orthogonal  move"  is 
made  because  of  an  indefinite  Hessian  matrix  -TP  is 
added  to  the  orthogonal  move  vector. 

■3  Steepost  descent  is  used  to  minimize  P-function. 

■4  The  method  for  minimizing  the  unconstrained  penalty 
function  is  McCormick's  modification  of  the  Fletcher- 
Powoll  method  as  reported  in  the  Fiocco-McCorni ck 
book  ■'.'  This  requires  function  values  and  first 

derivatives. 


b.  Default  Values  of  NAMELIST  Data  Items 

Variables  in  $CATA  are  given  the  following  values  in  SUM? 

prior  to  execution  of  the  NAMELIST  read. 


EPS  I 

RHOIN 

THETA 0 

RATIO 

TMMAX 

M 

N 

MZ 

X(I) 

NT(1) 

NT(2) 

NT(3) 

NT{4) 

NT{5) 

NT(6) 

NT(7) 

NT  (8) 

NT(9) 

NT (10) 

XEP1 

XEP2 

NEX0P1 

HEX0P2 


l.E-9 

1. 

l.E-6 

16. 

30. 

0 

0 

0 

0.  I  zi . 100 

3 

2 

1 

1 

3 
1 
1 
l 

l 

l.E-4 

0. 

1 

4 
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The  program  consists  of  a  main  routine,  two  control  subroutines 
(BODY,  FEAS),  twenty  one  special  purpose  subroutines  (CONVRG, 

ESTIM,  EVALU,  FINAL,  GRAD,  INVERS,  OPT,  OUTPUT,  PEVALU,  SAVE, 
REJECT,  RHOCCM ,  SECORD,  STORE,  T CHECK,  TIMEC,  SET,  S-IOVE,  DIFFl, 
DIFF2,  CKCKER)and  three  user  supplied  subroutines  (RESTNT, 
GRADl,  MATRIX).  Input  is  coordinated  by  the  main  routine. 

Output  is  from  the  main  routine  and  subroutines  BODY, 
CHOKER,  CONVRG,  ESTIM,  FEAS,  INVEHS,  OPT,  OUTPUT,  PUNCH,  TCHEOC, 
TIMEC. 

2)  Subroutines  Required; 

SUBROUTINE  BODY  coordinate?  all  subroutines. 


SUBROUTINE  CHOKER  computes  first  derivatives  of  objective  function 
using  GRADl  and  DIFFl. 

SUBROUTINE  CONVRG(Nl)  checks  for  convergence. 

SUBROUTINE  DIFFl(IN)  computes  first  derivatives  by  central  difference 

SUBROUTINE  DIFF2(IN)  comj  jtas  second  derivatives  by  central 
difference. 

SUBROUTINE  ESTIM  estimates  Lagrange  multiplier  (AJ  values  and  final 
solution  extrapolation. 

SUBROUTINE  EVALU  evaluates  objective  function  and  constraints. 

SUBROUTINE  FEAS  determines  feasibility  of  starting  point;  if  not 
feasible,  a  feasible  point  is  sought;  if  no  feasible  point  possible, 
error  message  printed. 

SUBROUTINE  FINAL(N2)  checks  for  convergence. 

SUBROUTINE  GP-AD(IS)  computes  first  derivatives  of  penalty  function. 
SUBROUTINE  INVERS(NSME)  solves  linear  system  of  equations. 

SUBROUTINE  OPT  performs  one  dimensional  search  by  Colden  Section 
method  (modified  Fibonacci). 

SUBROUTINE  OUTPUT  (K)  prints  out  results  of  each  iteration. 

SUBROUTINE  PEVALU  computes  value  of  penalty  function  and  dual. 
SUBROUTINE  REJECT  returns  stored  values  to  original  locations. 

SUBROUTINE  R'riCCOM  computes  initial  value  of  r. 

SUBROUTINE  SAVE  writes  a  NAMELIST  data  description  cf  the  or.i 
ter. iiticr.3  ts  ieftte  e-er.t  tf  ?:<e  :t_;n  ifcirt  i 


SUBROUTINE  SECQRD(IS)  computes  second  derivatives  of  penalty  function. 


SUBROUTINE  SET( UMAX)  stores  the  time  at  start  of  run. 

SUBROUTINE  STORE  stores  values  of  current  point. 

SUBROUTINE  TCHECK  checks  run  time  for  90%  of  maximum  specified. 
SUBROUTINE  TIMEC  checks  elapsed  time  of  run  for  possible  termination. 

SUBROUTINE  XMOVS  determines  saarch  directions  by  one  of  several 
options  (modified  Nevton-Raphson,  Steepest  Desceiit,  Fletcher-Powell). 

SUBROUTINE  R£STNT(IN,VA1.)  specifies  objective  function  and  constraints 
(user  supplied). 

SUBROUTINE  GRADl(IN)  evaluates  analytical  first  derivatives  for 
objective  function  and  constraints  (user  supplied). 

SUBROUTINE  MATRIX (J,L)  evaluates  analytical  second  derivatives  for 
objective  function  and  constraints  (user  supplied). 

3)  DIMENSION  Requirement s; 

The  general  subroutines  in  the  code  are  currently  dimensioned  to 
handle  a  problem  with  up  to  100  variable*  (N)  and  200  constraint* 
(M-fMZ).  This  should  be  sufficient  for  most  problems  attempted  with 
this  method.  If  the  limits  need  to  be  changed,  a  detailed 
description  of  the  various  parameters  are  given  elsewhere  (2). 

4)  Output 

The  program  title  is  printed  first  followed  by  a  listing 
of  the  parameters,  option  keys,  and  starting  point  infor¬ 
mation.  This  is  followed  by  intermediate  result  printouts 
until  convergence  is  achieved.  At  convergence,  the  final 
answer  information  is  printed.  Execution  time  is  listed 
with  the  intermediate  results.  In  the  event  that  the  ex¬ 
ecution  time  limit  set  by  TMMA.X  is  attained,  the  stopping 
point  information  needed  to  recommence  computation  at  the 
same  point  achieved  is  output  to  device  NS  in  NA^iELIST  form. 
This  file  (if  generated)  may  then  be  cataloged  and  used  as 
the  incut  NA.'.IILIST  data  in  a  subsequent  execution. 


134 


v"  V 


B.  Example  Problem 

The  following  minimization  problem  illustrates  the  effective 
use  of  the  SUMT  library  routines! 

Minimize  P(X1,X2)  m  (Xl-1. )2+(X2-2. )2 
Subject  to 

Cl  »  X1+X2S-2. 

C2  i  X2-X1=1. 

C3  «  XI— 0. 

C4  t  X2^0. 

Here,  the  objective  function  is  nonlinear,  while  all  constraints 
are  linear.  Constraint  Cl  is  an  inequality,  C2  is  an  equality, 
and  C3  and  C4  are  trivial  inequality  constraints.  Imposing  the 
constraint  format  defined  in  Section  A,  constraint  Cl  i3  re¬ 
written  as 


Cl  .  2.-X1-X2S0. 
and  constraint  C2  becomes 

C2  i  X2-X1-1.=0. 

After  selecting  initial  values  for  the  X  vector,  and  choosing 
to  employ  the  generalized  Newtin-Raphson  unconstrained  optim¬ 
ization  algorithm,  the  computer  run  is  coded  as  (see  next  cage) 
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<yrs,T35,ioioo,CMicoooo.  T820163, floyd 
ATTACH , SUMT , SUMT , ID=T8005  35  - 
LIBRARY, SUM? . 

FTN. 

LGO. 

*30R 

PROGRAM  MAIN( INPUT , OUTPUT , SAVE , TAPE5*IN?UT , TAPE6=CUT?UT , 

1  TAPE7-3AVE) 

COMMON/DSVC/NI ,  HO ,  .'IS 

NI»5 

NO=6 

NS-7 

CALL  SUMT 

STOP 

END 

SUBRCUTirrE  RE3TNT( IM, VAL) 

COMMON /SHAH  E /X ( ICO) , DEL( ICO) , A( 100 , 100 ) , N( 5 ) 

IF(IN)  10,10,20 

10  VAL=(X(l)-l.)*-*2+(X(2)-2.)**2 

RETURN 

20  GO  T0(21,22) , IN 

21  VAL=2.-X(1)-X(2) 

RETURN 

22  VAL=X(2)-X(l)-1 
RETURN 

END 

*ECR 

$DATA  N*2,M«i,MZ*l,X»0. ,1. ,NT(2)*1,NT(5)*2,NT(8)«2,NEX0P2«1  5 END 
•EOF 


The  solution  was  obtained  in  about  0.25  seconds  of  Central 
Processor  time.  The  final  printout  is 


FINAL  VALUE  OF  F  =>  5  .COCCOCCUS-l 

FINAL  X  VALUES 

X(  1)  »  1; .  99099998E-I  X(  2)  *  1 . 500CCCC0E0 


G{  1) 


3-7 G(  2)  * 


5 . 5L57GcC5E-ll 


136 


.  -v-r 


s 


(1)  Fiacco.A.V.  ,  and  G.P. McCormick.  Nonlinear  Sequential 
Unconstrained  Minimisation  Techniques.  John  Wiley  and 
Sons,  Inc.,  New  York,  1968. 

(2)  Mylander.W.C. ,  R.L.  Holmes,  and  G.P. McCormick.  A  Guide 
to  SUMT-Version  4i  The  Computer  Program  Implementing  the 
Sequential  Unconstrained  Minimization  Technique  for  Non¬ 
linear  Programming.  Research  Analysis  Corporation,  McLean, 
Virginia,  1971.  (AD731  391) 


(3)  Kuester.J.L. ,  and  J.H.  Mize.  Optimization  Techniques  with 
FORTRAN.  McGraw-Hill  Book  Company,  New  York,  1973. 

(4)  Fiacco.A.V.,  and  G.P.  McCormick.  "Computational  Algorithm 
for  the  Sequential  Unconstrained  Minimization  Technique  for 
nonlinear  Programming,"  Management  Science  (10),  601-617 
(July  1964). 

(5)  Bracken, J,,  and  G.P.  McCormick,  Selected  Applications  of 
Nonlinear  Programming.  John  Wiley  and  Sons,  Inc.,  New  York 


V 


The  following  program,  STCFCV,  is  a  Fortran  Four 
executive  program  for  SUMT.  This  program  is  run  to 
calculate  the  maximum  CF  for  given  contours  of  CV.  In 
its  current  form,  the  program  will  request  a  data  entry 
in  the  form  of  a  "namelist"  for  each  CV  contour  value. 
These  contours  are  set  at  two  unit  increments  and  the 
user  enters  the  starting  value. 


PROGRAM  STCFCV  ( INPUT ,  OUTPUT ,  flSHT ,  TAPES- 1 NPUT ,  !'APS6=0UTPUT ,  TAS  T 
I  =0SMT) 

COHMON/DEVO/NI,NO,N8 

COMMON/CNTUR/CVI 


C 

c 

c 


c 

c 

c 

c 


THIS  PROGRAM  USES  SUMT  TO  SOLVE  FOR  THE  MAXIMA  OF  RESPONSE  ONE, 
CF,  CONSTRAINED  TO  A  VALUE, CVI,  OF  RESPONSE  TUG.  THE  PROBLEM 
IS  FIJTHER  RESTRICTED  BY  THE  OPERATING  DOMAIN  UNDER  OONSIDERSTIOf 


Nl=>5 

N0=6 

MS=/ 

11=0 


13 


c 

c 

c 


WRITE (NO, 13) 

FORMAT  cm  ENTER  INITIAL  CVI  VALUE*****) 
READ(NI.*)CVI 
CVI=  1122,0 
DO  IV  11=1,34 


CVI--  CVI r  2.0 


11 


1.2 

IV 


WRITE (NO, 11)  CVI 
WRITE (NS, II)  CVI 

FORMAT ( 1X-F7.2,  *=*CVI*F0R  THE  FOLLOWING  RUN  *****  **.****:::***#*  ‘ 

CALL  SUMT 

WRITE (NS, 12)  CVI 

rnriUAT  f  -l  V  f  "7  •*)  |  *»  •/  «/  «/  y  «/  •/  "/  V  V  ■/  *'  «/“■"'  •/  •/  •'  “/  */  •'  V  */  V*'  */  */  *•'  * ' v  *'  •’  *  \ 

I”  U  ft  n  n  I  ■■  1  A  }  I  /  *  •.  f  .'it  /u  ftt  ,‘m  tm  s’ 4  /•»  .•  it  .• «  /«  Sm  /m  .-m  in  ■n  fn  /•»  in  i  «  ■  ,  /* .  <:  „  ■>,  in  tn  .•>  Jm  in in  .•»  ,  j  .  t  fn  ■  -  'j 

UR  I TE ( NO .12)  CVI 
CONTINUE 


STOP 

END 
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SUBWUT  f NE  RLHTNT (IN, UAL > 

ROHKGN/SUnRE/X ( ICO)  ,DEL(  100)  .A<  100.  j  00)  .M(f. : 

C:ftf?QN/CNTUR/r.VI 

c 

c 

IF ( IN)  10,10,20 

10  VAL--504 . 542366-0 ,4715619*X<1 )-0.295Sl 07*X < 2 ) -0 , 0?38839*X ( 3 ) 

C  -0 ,0273609*X ( 4 ) -0 , 0000003* ( X ( 1 ) **2 > +0 , 0000671 * ( X ( 2 ) **2 ) 

C  10 . 0000006* ( X ( 3 ) **? ) +0 . 0000005* ( X ( 4 ) **2 ) f 0  *  00000 1 5*X ( 1 ) *X ( 3 ) 

C  +0 , 00000Q9*X (1)*X(4) +0 . 0000007*X( 2 ) *X ( 3 ) +0 , 0000006*X < 2 ) *X ( 4 ) 

C  f 0 . 00000 13* X ( 3) *X ( 4 ) 

RETURN 

C 

20  GO  70(22,23,24, 25, 26, 27, 29, 29, 30), IN 
C 


22  VAL=  450  -  X<1> 
RETURN 

73  VAL  -  X  < 1 ) -300 
RETURN 

24  VA L=  750  -  X ( 2 ) 
RETURN 

25  UAL=  X ( 2 ) -450 
RETURN 

26  VAL=  1040  -  X(3) 
RETURN 

27  UAL=  X (31-720 
RETURN 


23  VAt=  424  -  X<4) 

RETURN 

29  Unl=  X ( 4 ) -200 
RETURN 

30  VAl= ( 1027 . 039-CVI ) +0 , 1658222*X (1 ) +0 . 1198222*X ( ? ) +«) . 0307M74*X < 

C  +0 . 001 9064*X <4)40. 0000054* ( X ( 4 ) **2 ) -0 , 0000153* ( X ( 1 ) **2 ) 

C  -0 , 0000040* ( X ( 2 ) **2 ) -0 , 0000075* ( X(3)**2 ) -0  - 0000129*X ( 1  ) *X , 


3) 


-0.0000 175‘X (1 )*X( 3)- 
RETURN 


■0 , 0000087JX  <  2 )  *X  ( 3 )  -0  •  0000054:! 


END 


An  example  of  the  SUMT  output  for  the  final  iteration 
to  maximize  CF  for  a  CV  value  of  1160  is  shown  on  the 
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NOTE:  Since  the  SUMT  is  a  minimizing  algorithm,  the  minimum  of  the  negative  CF 
was  calculated.  Therefore,  the  final  CF  objective  function,  f,  is  negative 
for  this  maximization. 


APPENDIX  E 


Analysis  Data 

This  appendix  first  contains  the  program  "surface" 
which  generated  data  for  the  CF  and  CV  plots.  Following 
the  program  are  the  enlarged  contour  plots  for  each  MOE 
and  the  enlarged  overlay  plots  used  in  Chapter  Four. 
Likewise,  the  enlarged  MOE  trade-off  graphs  are  included 
Finally,  Tables  IX  through  XII  contain  the  actual 
SUMT  data  for  each  of  the  MOE  trade-offs  in  Chapter  Four 
and  also  the  AF  trade-off  cases. 
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OPEN ( 11 ?F ILt-'P!  OT ' ) 

REWIND  U 
IW2=450 
WRITE  (11, 0.1) 

FORMAT ( 3X , ' CFF ' , 3X , '  W1',3X,'U2  ',3X,'  M3 ' » 3X , '  W4d 


LOOP  CF»W1  AND, M2  THROUGH  THEIR  RESPECTIVE  RANGES  TO  SOLVE  FOR  W 

DO  23,ICF=800,920,20 
DO  22 , I W 1 =300 , 450 , 25 
DO  21,1144=200,420,20 

M3  =  ( 1  ,'Q . 0289 )  * ( 1 CP -504 , 5424-0 . 471  A* I W 1-0 , 2995* I W2-0 .02791  .HI W 
C  +0.0001*1142**2) 

IF<W3.LT. 1041,0  .AND.  W3.0E, 720.0)  THEN 
MR ITEd 1 , 02 )  I OF , IU1 , I M2 , W3 , IW4 
FORMAT (2X,I3,2X,I3,2X,I3, 2X , FA ♦ 1 , 2X ,13) 

ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 

URITEdl  ,03) 

FORMAT { 5X , '  *-*-*-*-*-*-#-#-*-*-*-*-*-*-*-*-#-*-*-*-*-*-*-*-*-* 

WRITE < 11,04) 

FORMAT  (3XV  CV  ',2k,-’  141  ',2X,'  M2  '  ,2X,'  M3  ',2*.'  U4') 

LOOP  CV,W3, AND  M2  THROUGH  THEIR  'RESPECTIVE  RANGES  TO  SOLVE  FOR  U 


DO  33, ICV=1 155, 1225, 10 
DO  32,  [141=300,450,25 
DO  31, 1144=200,420, 20 

U3= ( 1/0 . 0303 ) * ( ICV-i 027 . 039-0 , 1658* tWl-0 . 1 198* 1W2-0 . 00 19* [U4 

IF(W3»LT. 1041.0  .AND.  W3.GE.720  0)  THEN 

URITEdl, 05)  ZCV,IWJ  ,IU2,W3,IU4 

FORMAT ( 2X , 14 ,2X , 13 , 2X , 13 , 2X , F6 . 1 ,2X , 13) 

ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 

PRINT  *, 'RUN  COMPLETE-CONTOUR  MAP  DATA  IN  FILE-PLOT  ' 

CLOSE  (11.) 


(1175,400) 


(1175,410) 


(1175,420) 


(1175,430) 


Contour  Plot 


400 


Fig  14.  Contour  Plot 


JV*  4l» 
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Fig  20 
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Hounded  Regression  Coefficients 


Fig  22.  Maximum  CF  subject  to  CV 
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Maximum  CV  subject  to  CF 
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TABLE  IX 

SUMT:  Maximum  CF  given  CV 


Max  OF 

CV 

W2 

"3 

W4 

791.14 

1154 

300.90 

450.00 

727.41 

390.96 

797.32 

1156 

314.39 

450.00 

720.06 

393.68 

803.05 

1158 

326.44 

450. CO 

720.04 

395.05 

809.27 

1160 

338.04 

450.00 

724.90 

420.98 

814.42 

1162 

350.43 

450.00 

720.65 

396.45 

820.30 

1164 

362.52 

450.00 

720.06 

403.47 

825.45 

1166 

374.46 

450.00 

721.93 

384.54 

831.71 

1168 

386.64 

450.00 

720.03 

404.99 

836.54 

1170 

399.13 

450.00 

720.04 

366.74 

840.10 

1172 

412.27 

450.01 

720.03 

272.27 

848.70 

1174 

422.86 

450.00 

720.04 

401.55 

855.36 

1176 

434.94 

450.00 

720.05 

400.34 

860.47 

1178 

446.60 

450.00 

720.09 

421.44 

863.35 

1180 

450.00 

457.23 

741.20 

392.99 

866.67 

1182 

450.00 

473.03 

744.62 

393.85 

870.81 

1184 

450.00 

493.35 

729.32 

413.37 

374.32 

1186 

450.00 

511.60 

722.72 

421.92 

876.92 

1188 

450.00 

523.33 

742.77 

409.83 

879.39 

1190 

450.00 

536.63 

757.00 

393.54 

883.63 

1192 

450.00 

558.53 

734.97 

422.14 

886.68 

1194 

450.00 

575.42 

734.21 

421.97 

888.70 

1196 

450.00 

587.05 

755.41 

397.80 

892.42 

1198 

450.00 

605.34 

747.59 

423.83 

895.31 

1200 

450.00 

622.92 

744.21 

422.81 

898.31 

1202 

450.00 

642.42 

733.31 

422.44 

899.78 

1204 

450.00 

649.79 

771.48 

391.81 

903.12 

1206 

450.00 

656.37 

770.46 

415.01 

906.12 

1208 

450.00 

686.26 

757.65 

421.90 
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Table  IX  continued 


TABLE  X 


SUMT:  Maximum  CV  given  CF 


Max  CV 

CF 

W1 

W2 

w3 

W4 

1156.85 

790 

300.69 

456.56 

806.43 

225.08 

1158.72 

792 

300. 10 

462.81 

846.26 

219.77 

1160.91 

794 

300.03 

467.21 

901.65 

203.39 

1160.65 

796 

300.05 

478.98 

844.41 

249.25 

1163.77 

798 

300.03 

480.44 

942.41 

209.29 

1165.03 

800 

300.06 

492.10 

938.52 

201.72 

1165.48 

802 

300.06 

502.08 

912.54 

230.11 

1166.98 

804 

300.44 

510.39 

927.30 

222.20 

1168.28 

806 

300.38 

518.72 

937.08 

227.26 

1169.72 

808 

300.07 

536.28 

918.73 

203.67 

1170.33 

810 

300.24 

542.09 

912.69 

239.71 

1173.13 

812 

300.38 

545.69 

990.97 

203.77 

1174.13 

814 

300.03 

561.75 

963.09 

203.99 

1175.35 

816 

300.04 

573.60 

956.37 

205.12 

1176.17 

818 

300.01 

586.11 

933.72 

220.24 

1178.12 

820 

300.02 

589.45 

983.99 

218.42 

1179.31 

822 

300.39 

601.34 

974.51 

218.53 

1180.71 

824 

300.02 

616.73 

962.25 

213.37 

1182.49 

826 

300.02 

627.16 

980.34 

202.34 

1183.86 

828 

300. 29 

636.29 

987.55 

206.58 

1184.49 

830 

300.02 

648.30 

960.61 

238.85 

1186.93 

832 

300.02 

659.63 

997.97 

204.95 

1188.22 

834 

300.04 

675.27 

978.87 

205.20 

1189.88 

836 

300.01 

685.26 

994.07 

204.56 

1191.43 

838 

300.00 

697.02 

998.58 

205.44 

continued 


Table  X  continued 


1192.24 

840 

1193.80 

842 

1195.41 

844 

1197.93 

846 

1199. 14 

848 

1199.01 

850 

1199.33 

852 

1201.03 

854 

1201.72 

856 

1202.64 

858 

1202.74 

860 

1203.20 

862 

1203.58 

864 

1204.23 

866 

1205.05 

868 

1206.74 

870 

1207.41 

872, 

1208.07 

874 

1208.78 

876 

1208.60 

878 

1209.29 

880 

1210.87 

882 

1210.73 

884 

1211.67 

886 

1211.99 

888 

1212.56 

890 

1213.94 

892 

1214.94 

894 

300.01 

706 . 00 

300.24 

718.39 

300.32 

731.12 

300.02 

747.69 

300.58 

749.83 

305.98 

748.98 

310.22 

749.25 

312.26 

749.99 

317.12 

749.15 

322.60 

749.95 

325.91 

749.14 

330.81 

749.92 

335.20 

749.95 

339.66 

749.95 

343.32 

749.29 

347.41 

749. 19 

350.57 

749.83 

354.37 

749.94 

358.96 

749.95 

364.25 

749.94 

368.66 

749.95 

371.46 

749.57 

376.69 

749.95 

380.79 

749.93 

384.72 

749.37 

392.34 

747.70 

393.18 

749.95 

395.98 

749.68 

988.06 

237.35 

939.37 

236.49 

991.56 

235.77 

1012.72 

203.01 

1039.14 

226.71 

1007.94 

243.83 

993.86 

256.96 

1035.67 

247.04 

1035.57 

241.01 

1034.44 

217.00 

1020.92 

251.00 

1006.25 

251.01 

994.36 

260.70 

991.37 

260.00 

1000.6“ 

263.64 

1036.13 

230.08 

1037.53 

243.50 

1037.51 

250.42 

1036.07 

245.94 

1000.83 

264.71 

999.48 

263. 19 

1038. 10 

249.54 

1002.66 

267.59 

1011.75 

260.63 

1002.17 

278.75 

988.09 

244.99 

1018.62 

259.10 

i036 . 51 

266. 23 
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Table  X  continued 


1215.65 

896 

1216.46 

898 

1217.10 

900 

1217.32 

902 

1217.47 

904 

1219.34 

906 

1219.84 

908 

1220.04 

910 

1220.63 

912 

1221.94 

914 

1222.20 

916 

1222.20 

918 

1223.39 

920 

1223.57 

922 

1224.33 

924 

400.98 

749.94 

404.98 

749.95 

408.89 

749.95 

414. 12 

749.25 

419.60 

749.28 

422.18 

749.95 

425.93 

749. 16 

428.75 

749.05 

433.66 

749.95 

437.33 

749.93 

439.04 

749.90 

442.96 

749.07 

445.30 

749.72 

448.49 

749.94 

449.99 

749.98 

1032.46 

256.32 

1037.14 

255.45 

1036.68 

261.57 

1017.83 

267.99 

993.09 

272.62 

1038.37 

250.20 

1036.78 

264.29 

1026.34 

299.74 

1015.91 

294.55 

1039.70 

279.57 

1036.08 

326.19 

1016.19 

356.60 

1039.40 

361.37 

1025.31 

392.38 

1039.90 

423.55 

TABLE  XI 


SUMT:  Maximum  CF  given  AF 


Max  CF 

AF 

W, 

X 

W2 

W3 

W4 

796.13 

1440 

321.83 

452.41 

727.13 

200.04 

799.73 

1445 

330.88 

450.04 

720.76 

200.04 

802.53 

1450 

336.80 

450.05 

721.00 

200.02 

804.89 

1455 

341.44 

450.03 

727.01 

200.03 

808.22 

1460 

348.89 

450.11 

720.45 

200.03 

810.83 

1465 

354.23 

450.07 

723.60 

200.03 

813.63 

1470 

360.08 

450.22 

723.75 

200.21 

816.69 

1475 

366.78 

450.08 

721.63 

200.03 

819.61 

1480 

372.93 

450.28 

720.74 

200.03 

822.43 

1485 

378.83 

450.39 

721.03 

200.23 

824.85 

1490 

383.49 

450.69 

726.90 

200.03 

827.98 

1495 

390.44 

450.32 

724.42 

200.00 

831.09 

1500 

397.13 

450.56 

720.81 

200.37 

834.03 

1505 

403.47 

450.28 

721.48 

200.00 

836.54 

1510 

408.15 

451.23 

724.47 

200.64 

838.33 

1515 

409.29 

456.35 

730. 19 

202.12 

842.83 

1520 

422.09 

450.40 

721.21 

200.00 

844.41 

1525 

421.73 

457.78 

729.56 

200. 14 

847.60 

1530 

429.29 

456.07 

728.20 

200.47 

851.31 

1535 

439.78 

450.37 

726.31 

200.01 

854.77 

1540 

447.52 

450.32 

720.09 

200.03 

857.14 

1545 

449.97 

456.00 

721.49 

200.08 

858.87 

1550 

449.73 

464.28 

726.03 

201.45 

860.38 

1555 

449.98 

469.43 

735.01 

204.39 

862.92 

1560 

449.98 

482.99 

733.27 

200.12 

864.93 

1565 

449.98 

492.74 

735.68 

200.31 

867.07 

1570 

449.98 

503.75 

735.30 

200.27 

continued 
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868.61 

1575 

449.95 

509.38 

751.41 

200.02 

871.20 

1580 

450.00 

524.54 

739.16 

200.34 

873.10 

1585 

449.98 

533.53 

746.20 

200.36 

874.19 

1590 

449.98 

534.88 

772.94 

202.43 

876.57 

1595 

449.98 

547.93 

772.99 

200.31 

878.47 

1600 

449.99 

556.52 

783.66 

200.22 

880.76 

1605 

449.97 

569.31 

782.06 

200.24 

883.62 

1610 

449.98 

588.61 

760.89 

200. 19 

884.32 

1615 

449.98 

584.25 

809.42 

202.98 

885.76 

1620 

449.98 

588.42 

829.18 

207.62 

888.59 

1625 

449.98 

605.85 

825.34 

203.04 

891.68 

1630 

449.89 

628.21 

803.24 

200.03 

886.08 

1635 

449.99 

574.87 

852.24 

282.03 

895.77 

1640 

449.99 

646.63 

835.78 

200.57 

898.28 

1645 

449.96 

662.83 

830.71 

200.70 

900.67 

1650 

449.98 

676.19 

838.94 

200.09 

903.28 

1655 

449.99 

694.81 

826.61 

200.59 

904.98 

1660 

449.99 

702.35 

836.39 

209.36 

905.97 

1665 

449.98 

703.39 

847.57 

227.58 

910.13 

1670 

449198 

731.02 

866.00 

206.81 

910.13 

1675 

449.98 

731.02 

866.00 

206.81 

913.53 

1680 

449.97 

741.14 

908.44 

231.17 

914.25 

1685 

450.00 

736.25 

932.99 

257.05 

912.59 

1690 

449.99 

730.61 

854.22 

309.40 

918.45 

1695 

449.99 

745.65 

1003.14 

285.49 

917.32 

1700 

449.98 

749.40 

908.52 

323.19 

920.36 

1705 

449.97 

747.56 

1010.79 

336.56 

921.27 

1710 

449.99 

749.51 

1007.47 

361.80 

919.99 

1715 

449.98 

745.02 

949.79 

399.31 

920.62 

1720 

450.00 

749.79 

923.95 

423.48 

TABLE  XII 


SUMT:  Maximum  CV  given  AF 


Max  CV 

AF 

W1 

w2 

W3 

W4 

1157.03 

1440 

305.07 

458.14 

784.02 

201.15 

1157.90 

1445 

309.22 

465.75 

759.82 

206.86 

1158.70 

1450 

310.47 

470.11 

761.75 

211.15 

1159.54 

1455 

310.33 

468.82 

794.98 

211.38 

1160.71 

1460 

314.70 

470.24 

804.15 

207.80 

1162.22 

1465 

317.56 

483.48 

786.11 

201.47 

1163.16 

1470 

317.63 

483.01 

818.51 

201.20 

1163.59 

1475 

318.86 

485.61 

814.87 

213.69 

1165.36 

1480 

323.39 

499.68 

794.19 

202.81 

1166.34 

1485 

325.21 

497.28 

825.53 

202.18 

1167.59 

1490 

324.46 

509.41 

822.75 

201.10 

1168.66 

1495 

327.98 

512.84 

825.46 

201.35 

1169.66 

1500 

336.48 

508.32 

829.77 

201.29 

1170.85 

1505 

336.82 

516.96 

832.78 

201.06 

1171.97 

1510 

339.98 

521.66 

833.94 

201.27 

1173.39 

1515 

331.14 

545.48 

834.94 

201.13 

1174.09 

1520 

347.00 

528.89 

836.61 

204.07 

1175.85 

1525 

331.89 

563.27 

841.72 

201.11 

1176.55 

1530 

354.52 

538.08 

840.45 

201.47 

1177.26 

1535 

346.41 

553.92 

844.65 

214.32 

1179.63 

1540 

335.71 

588.97 

843.71 

201.28 

1179.22 

1545 

364. 23 

545.91 

843.35 

216.70 

1182.43 

1550 

330.94 

614.37 

861.49 

201.38 

1183.72 

1555 

334.84 

618.73 

865.49 

201.00 

1185.29 

1560 

329. 17 

636.46 

877.90 

201.44 

1185.34 

1565 

341.27 

623.89 

862.55 

217.12 

continued 
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1186.87 

1570 

347.46 

625.62 

872.48 

212.30 

1188.43 

1575 

364.29 

619.93 

855.41 

201.20 

1190.64 

1580 

344.51 

660.93 

874.08 

201.29 

1190.45 

1585 

387.94 

602.55 

861.10 

201.48 

1191.17 

1590 

396.31 

602.60 

856.67 

200.08 

1195.90 

1595 

336.58 

706.61 

910.03 

200.52 

1197.38 

1600 

342.57 

706.17 

927.26 

200.93 

1199.17 

1605 

341.86 

719.82 

936.12 

201.36 

1201.69 

1610 

339.21 

730.43 

991.08 

200.92 

1201.03 

1615 

369.62 

704.46 

906.78 

201.14 

1204.96 

1620 

347.77 

745.36 

993.11 

201.11 

1205.27 

1625 

361.87 

745.66 

926.02 

200.80 

1206.77 

1630 

373.00 

731.89 

968.28 

202.38 

1207.57 

1635 

384.88 

735.04 

918.18 

200.99 

1209.94 

1640 

390.98 

723.08 

1008.67 

202.00 

1212.01 

1645 

390.98 

743.71 

995.71 

200.74 

1213.84 

1650 

396.73 

749.46 

1001.79 

200.09 

1215.48 

1655 

405.39 

747.86 

1015.02 

200.05 

1216.79 

1660 

413.96 

749.80 

1003.54 

200.24 

1218.21 

1665 

422.19 

749.66 

1005.76 

201.92 

1220.05 

1670 

430.96 

749.81 

1017.72 

200.07 

1219.98 

1675 

438.21 

749.96 

975.28 

212.82 

1221.39 

1680 

449.53 

749.63 

961 . 53 

209 . 07 

1222.87 

1685 

445.90 

749.41 

1021.35 

235.69 

1220.88 

1690 

447.89 

749.93 

94S.86 

268.78 

1222.85 

1695 

447.86 

749.95 

1011.35 

235.78 

1223.06 

1700 

446.64 

749.92 

1023.81 

315.50 

1219.17 

1705 

449.31 

749.92 

880.35 

355.75 

1222.55 

1710 

446.30 

749.94 

1005.20 

372.56 

1224.09 

1715 

449.86 

749.96 

1035.14 

387.36 

1223.22 

1720 

449.18 

749.94 

1008.83 

417.94 
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